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The Ferromagnetic Faraday Effect 


at Microwave Frequencies 


and its Appli ‘ations 


The Microwave Gvrator 
BY C. L. HOGAN 


A new microwave circuit element dependent on the Faraday rotation of a 
polarized wave has been developed. The element violates the reciprocity 
theorem and, because it shares this property with a gyroscope and because 
it is dependent on gyromagnetic resonance absorption, it has been termed 
a microwave gyrator. It is a low-loss broadband device with many applica- 
tions. Among these are one-way transmission systems, microwave circula- 
tors, microwave switches, electrically controlled variable attenuators and 
modulators. 

The microwave gyrator has been realized by making use of the Faraday 
rotation in pieces of ferrite placed in the waveguide. Polder has previously 
shown, in his analysis of the gyromagnetic resonance phenomenon, that 
ferromagnetic substances should show appreciable Faraday rotations at 
microwave frequencies. In the present study, Polder’s analysis has been 
extended to include a wave being propagated through a ferromagnetic sub- 
stance with dielectric and magnetic loss, and data are presented which give 
experimental verification of the theory. In addition an experimental tech- 
nique is described which may be of some interest in studying the properties 
of ferrites at microwave frequencies. 
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Photograph of the experimental setup shown diagrammatically in Fig. 5 


INTRODUCTION 


; » . ry 1 ‘ 

In a recent series of articles, Tellegen’ has discussed the possible 
application < of a new circuit element which he calls a gyrator. He defines 
the ideal gyrator, in principle, as a passive four-pole element which is 


described bv: (see Fig. 1 


1 (1) 


i 


Since the coefficients above are of opposite sign, the gvrator violates the 
theorem of reciprocity. Any network composed of the usual electrical 
circuit’ elements—resistors, inductors, capacitors, and transformers 

will satisfy the theorem of reciprocity. In simple terms, this theorem 
states that if one inserts a voltage at one point in the network and 
measures the current at some other point, their ratio (called the transfer 
impedance will be the same if the positions ol voltage and current are 


ter hin ie dl In the yvrator, however, this transter Impedance tor one 


ee 


| tour pole 
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direction of propagation 1s the negative of that for the other direction of 
propagation. Essentially this moans that a 180° phase difference exists 
between the two directions of propagation. For this reason it has been 
suggested that the element could be more aptly called a directional 
inverter.” 

Network synthesis today is based upon the existence of four basic 
circuit elements: the capacitor, the resistor, the inductor, and the 
ideai transformer. It is apparent that the introduction of a fifth circuit 
element, the gyrator, would lead to considerably improved solutions for 
many network problems. In fact, Tellegen’ has shown that the synthesis 
of resistanceless four-pole networks would be much simplified by its 
introduction. In addition, McMillan” has shown that it would be possible 
to construct a one-way transmission system if a gyrator were available, 


: 4 . ° . ' 
and Miles’ has shown that it would be possible by use of a gyrator to 


construct a network which is equivalent to a Class A vacuum tube 
amplifier circuit. While the realizable power gain of these gyrator circuits 
is necessarily always not greater than unity, many other networks in- 
cluding gyrators are possible which have properties analogous to vacuum 
tube circuits and some of these may be of practical importance. Since 
this new element offers such interesting possibilities in network synthesis, 
a study has recently been made in these Laboratories of possible methods 
for realizing the gyrator. 

A gyrator was employed by Bloch® in his measurement of the mag- 
netic moment of the proton. Bloch made use of the phenomenon that 
if two crossed coils with a mutual core are adjusted so that there is zero 
mutual inductance between them and if a steady magnetic field is 
applied perpendicular to the axes of both coils, then an ac voltage ap- 
plied to one of the coils will induce a voltage in the second due to the 
gyromagnetic resonance phenomenon. This induction is_ ordinarily 
extremely small unless the magnetic field is adjusted so that the exciting 
frequency coincides with a gyromagnetic resonance frequency of the 
material which forms the mutual core of the two coils. In Bloch’s ex- 
periment, the magnetic field was held constant and the exciting fre- 
quency was adjusted until it coincided with the gyromagnetic resonance 
frequency of the proton. If they were wound over a paramagnetic or 
ferromagnetic material, the two crossed coils would form a gyrator 
when the magnetic field was adjusted so that the frequency of the 
exciting field coincided with the gyromagnetic resonance of the unpaired 
electrons. The fact that this structure constituted a gyrator was first 
recognized by Tellegen’ and has been discussed by Beljers and Snoek' 
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in a paper which gives a very satisfying physical model with which to 
interpret gyromagnetic phenomena occurring within ferrites. 

Physical analysis indicates that the properties of ferromagnetic 
materials can be explained by assuming that the electron behaves as if 
it were a negatively charged sphere which is spinning about its own 
axis With a fixed angular momentum. This rotation of charge imparts 
to the electron a magnetic moment which is a function of the electric 


‘harge on the electron, the angular velocity of the electron, and its 


size. Thus the electron behaves as if it were a spinning magnetic top, 


whose magnetic moment lies along the axis of rotation, and its behavior 
can be understood by considering a spinning gyroscope suspended in 
gimbal rings at a point not coinciding with its center of gravity. If a 
gyroscope, thus supported in a gravitational field, is lifted away from 
its position of minimum potential energy and then released, it will not 
return to the position of minimum energy but will precess about the 
vertical axis. This is illustrated in }“g. 2 where the spinning gyroscope 
makes an angle #@ with the vertical 2, “AXIS, Its equilibrium motion, In 
the absence of damping, is a precessional motion about the vertical 
axis with a velocity w». 


If the gyroscope be regarded as initially hanging vertically downward 














xX 


Precessional motion of a gyroscopic pendulum in a gravitational 
Precessional motion of a gyroscopic pendulum in a gravita- 
cillates between the directions A and B 
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as indicated in Fig. 3 and then a gravitational force is suddenly made to 


act along the y axis so that the net gravitational force acts along A, 
it is obvious that the gyroscope will begin to precess about the gravita- 
tional field direction as indicated by the small dotted circle. However, if 
after completing a half cycle, the horizontal component of the gravita 
tional field is reversed so that now the net gravitational field acts along 
the vector B, the gyroscope will begin to precess about B as indicated 
by the intermediate size dotted circle. If the horizontal component of 
the gravitational field is again reversed after the gyroscope completes 
another half cycle in its precession, the gyroscope will again begin to 
precess about the direction A and the actual path of the precessional 
motion will be along the path a-b-c-d. If this process is continued in- 
definitely, the gyroscope will precess in larger and larger circles around 
the vertical until the damping becomes large enough to contain the 
gyroscope in some equilibrium circle (assuming that the damping is 
large enough to accomplish this). 

The above model affords a classical picture which can be used quite 
readily to describe the motion of the electrons in a ferrite. If the ferrite 
is initially saturated along the z axis by a steady magnetic field, the 
electrons will come to rest with their magnetic moments lying along 
the 2 axis, as the gyroscope in Fig. 3. If now an alternating magnetic 
field is applied along the y axis, the electrons will begin to precess in 
larger and larger circles about the z axis until they finally reach some 
equilibrium position under the influence of the magnetic fields and the 
damping. Thus it is apparent in the gvromagnetic resonance experi- 
ments described above why an alternating field applied perpendicular 
to a steady magnetic field in a ferrite will give rise to a varying flux 
perpendicular to both the steady field and the alternating field. It is 
also apparent why the alternating flux along the x axis is 90° out of 
phase with the alternating flux along the 7 axis. Since precession of the 
top will always be in the same direction regardless of whether the alter- 
nating field is applied along the a or y axes, consideration of Fig. 3 
makes it apparent how the two crossed coils with ferrite at their center 
can constitute a gvrator which violates the reciprocity relation in a 
manner described by Equations (1). To the present time, however, no 
practical circuit element making use of this phenomenon has been 
constructed because the coefficient of coupling between the coils is 
always small, even in the vicinity of the resonant frequency, and also 


because the losses im the materials available are so high in the vicinity 
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of resonance that the insertion loss of such a device would be prohibitively 
large 

MeMillan in his original article,” showed that a gvrator could be 
realized by means of mechanically coupled piezo-electric and electro- 
magnetic transducers. Later, McMillan’ pointed out that a gyrator 
could be realized by means of the Hall effect in a square plate of bismuth, 
is Was also predicted by Casimir.” Another similar possibility would be 
an electrical-electrical coupling through a gy roscopic link. A gvrator has 
heen built by W. P. Mason of these Laboratories which makes use 
of the Hall effect in a erystal of germanium, This gvrator showed an 
insertion loss somewhat higher than the theoretical loss of 12.3 db. 
R. O. Grisdale of these Laboratories suggested that these losses could 
be greatly reduced if the same Hall effect principle were applied to a 
vacuum tube which contained four electrodes which could both emit 
ind collect electrons. This device is no longer passive, but such a strue- 
ture has been built and showed an insertion loss of about 7 db, only 
slightly higher than the theoretical loss which would be expected from 
this geometry 

In view of the substantial losses found to exist in the earlier forms of 
gvrator discussed above, a study of other ‘‘anti-reciprocal” phenomena 
which might lead to the realization of a relatively low loss gyvrator was 
undertaken 


It has long been known that the Faraday rotation of the plane of 


polarization in optics is anti-reciprocal. In order to observe the Faraday 


rotation, polarized electromagnetic waves must be transmitted through 
a transparent isotropic medium parallel to the direction of the lines of 
force of a magnetic field. The effeet is usually produced by placing the 
material along the axis of a solenoid. The rotation is “positive” if it is 
in the direction of the positive electric current which produces the field 
and “negative” if in the opposite direction. All optically transparent 
substances show the Faraday rotation. 

Its anti-reciprocal property distinguishes the Faraday effect from 
opti il rotations caused by birefringent crystals, or by the Cotton- 
Mouton effeet, which are reciprocal. That is, if a plane polarized light- 
wave is incident upon a birefringent erystal in such a manner that the 
plane of polarization is rotated through an angle @ in passing through 
the crystal, then this rotation will be cancelled if the wave is reflected 
back through the erystal to its source. In the Faraday rotation, however, 
the angle of rotation is doubled if the wave is reflected back along its 
path. Hence, if the length of path through the “active” material is 


} 


adjuste dso as to give a 90° original rotation, the beam on being reflected 
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will have its plane of polarization rotated a total of 180° in passing in 
both directions through the material. Thus, the Faraday rotation in 
optics affords an anti reciprocal relation quite analogous to the antti- 
reciprocal property of the gyrator postulated by Tellegen. 

Lord Rayleigh’ described a one-way transmission system in optics 
Which makes use of the Faraday rotation. Lord Rayleigh’s ‘‘one-way”’ 
system consisted of two polarizing Nicol prisms (oriented so that their 
planes of acceptance made an angle of 45° with each other), with the 
material causing the Faraday rotation placed between them. Thus, 
light which was passed by the first crystal and whose plane of polariza- 
tion was rotated 45° would be passed by the second crystal also. But, 
in the reverse direction, the rotation would be in such a sense that light 
which was admitted to the system by the second crystal would not be 
passed by the first. 

Although Rayleigh’s one-way transmission system can be actually 
realized, it is experimentally difficult since most substances show ex- 
tremely small Faraday rotations. In fact, large rotations for transparent 
substances in the optical region are of the order of one degree per cm 
path length for an applied magnetic field of 1000 oersteds. To realize 
a rotation of 45° would require maintaining a field of 1000 oersteds over 
a distance of approximately one-half meter. The Faraday effect in 
ferromagnetic substances, however, is unique in that it shows rotations 
many orders of magnitude greater than the rotations exhibited by any 


other substances. For instance, Kénig'’ reports rotations of 382,000°/em 


by passing light through thin layers of magnetized iron. These data, of 
necessity, however, were taken on extremely thin sections and the total 
rotation obtained for any specimen did not exceed 10°. In order to 
obtain appreciable rotations in a device of practical size, it is necessary 
to obtain a material which shows a rotation at least intermediate be- 
tween those reported for iron and other ordinary materials. In addition, 
in order to make effective use of these rotations, the material must be 
transparent to the radiation which is being used. 


THEORY OF THE FERROMAGNETIC FARADAY EFFECT 


Polder" has shown in his analysis of the ferromagnetic resonance 
phenomenon, that a plane electromagnetic wave at microwave fre 
quencies should show appreciable Faraday rotation when propagated 
through a ferromagnetic material which is magnetized in a direction 
parallel to the direction of propagation of the wave. Polder has neglected 
both magnetic and dielectric losses in his analysis and although for the 


ferrites which are of greatest interest, this approximation is quite 
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valid, nevertheless the more complete theory is developed below. The 
exact theory of this phenomenon should, of course, be approached 


through quantum mechanics, but since the classical theory, in this 


particulal case, gives a result as satisfactory as the quantum thebry 


and since it lends itself more aptly to a fundamental physical inter- 
pretation of the phenomenon, it is the classical theory which is developed 
here Quantum mechanically, Faraday rotations in the optical region are 
recounted for by the Zeeman splitting of the spectral lines. 

The classical model which proves quite adequate for the description 
of ferromagnetic resonance is that illustrated in Figs. 2 and 3 which 
regards the electrons of the material which contribute to the magnetism 
us being spinning magnetic tops. The angular momentum of each 
electron Is 


(2) 


Angular momentum ol electron vin cm, sec 
h Planck’s constant (6.62 «& 10 erg sec) 


Phe magnetic moment which arises due to this rotation ts: 
eh 
arm 
where: 


Magnetic moment of electron Sohn magneton ) 
Charge on electron (4.80 & 10 K.S.U.) 

Mass of electron (9.10 & 10 

Velocity of light (38 & 10° em 


The so-called gvromagnetic ratio of the electron is the ratio of these 
quantities and is given by: 


‘ Va 


2m e 


If a steady magnetic field is applied to the sample such that the elec- 
tron sees an effective field /7, then a torque will be applied to the electron 
Which tries to turn the electron so that its magnetic moment lies along 
the field direction. However, as indicated in Fig. 2, the electron will 
precess around the field direction until damping forces dissipate the 


energy of precession, The equation of motion ot the electron is: 


dJ 1 dup 


au MH 
: i dt dt 
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The equation of motion of the magnetization per unit volume can 
thus be written: 
dM 


— ‘ 
di yM XH 


where: 


M = Magnetization of medium 


H = Macroscopic internal magnetic field 


The above equation, however, does not include damping. The damping 
force, regardless of its origin, must be so introduced into the above 
equation that it tends to cause the electron’s axis of rotation to line up 
with the field direction. It has been shown by Yager, Galt, Merritt 
and Wood’* that the shape of the resonance absorption line can be 


accounted for if the damping term is introduced in the following way: 


dM 
dt 


yM XH - (Mx Of Xx TH] 

The vector MX (MX H) is simply a vector which is in the proper 
direction to act as a damping force (torque) and the coefficient is chosen 
so as to give the correct units along with the parameter, a, which must 
be determined experimentally and which gives the magnitude of the 
damping torque. 

Equation 7 then is the equation of motion of the magnetization of 
an arbitrarily shaped body under the action of an arbitrary internal 
field, H. In the appendix, it is shown that if a steady magnetic field, 
H,, is applied along the z axis and then a small alternating field is ap- 
plied in an arbitrary direction to a sample which is infinite in size, the 
equation relating the resulting alternating flux density, 6, and the ap- 
plied alternating field, h is: 


b, = ph. — jKhy 


b, = jKh, + ph, 


where 
(9) 
kK = Kh’ — jk” (10) 


Kquations which give » and K in terms of the applied magnetic field 
and fundamental atomic constants are given in the appendix. 
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I:quations (8) are easily interpreted in terms of the spinning gyroscope 
model of Fig. 3. If magnetic losses had been ignored (.e. a 0) then 
both w and A would have been real. Under this condition, it is seen that 
if an alternating field, hy , is applied along the y axis, then an alternating 
flux, b, , is created along the y axis which is in phase with A, , and an 
alternating flux, b, , is created which is 90° out of phase with A, . Reci- 
procity between the 2 and y directions would demand that both terms 
containing JA should have the same sign. Thus, Equations (8) give a 
quantitative expression for the results which were previously qualita 


‘ly deduced by means of the electronic model illustrated in Figs. 2 


veguide is filled with a ferromagnetic material such as a ferrite 

‘na steady magnetic field is applied along the axis of the 
waveguide, 1l 1S necessary in order to describe this wave to find a solu- 
tion to Maxwell’s equations which is consistent with Equations (8) 
and in which 6, A, 2 and D are all proportional to exp [ywt — Te]. This 
proble m Is not solved exactly. However, in the appendix a solution is 
obtained for an infinite plane wave. It is found that the ferromagnetic 
medium ean support only a positive or a negative* circularly polarized 
wave or a combination of both, It is also shown in the appendix that the 
propagation constants lor these two « ircularly polarized waves are dift- 


ferent and are given by the following eXpressions: 


+ A )\le] 


Propagation constant 


Angular frequeney of wave 


Velo ity of light in unbounded space (3 & 10° em/see 


Complex dielectric constant of medium 
In Equations (11) and (12) it is apparent that the effective perme- 
ability of the medium to a positive circularly polarized wave, for in- 


is used here \ re the positive component 1s the com 


thie positive electric current which 
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stance, is given by the expression (u + A), and not by the usual perme- 
ability, b,/hz = u. It is also apparent that the quantity » + A can 
vary over wide limits in the vicinity of the ferromagnetic resonance. 
For this reason, care must be taken in interpreting permeability data 
for ferromagnetic materials which now occur in the literature and which 
were obtained by means of impedance measurements at microwave 
frequencies, since the above equations indicate that this method does 
not measure the same quantity that is measured at low frequencies by 
means of a toroidal sample overwound with two coils. The low fre- 
quency measurement of permeability obviously measures the quantity 
which is designated as uw in Equation (8). 

If Equations (11) and (12) are solved for the attenuation constants, 
a,, and the phase constants, 8, , the following results are obtained: 


tan 6,,{4 tan 6g + tan 6,,(1 + tan? 6g)| + tan 


7 
1 — tan 6,, tan | 


4 4] (1 + tan 6,,[4 tan 6g + tané,(1 + tan? 64)[ + tan ~ 


+ 1+ tan 6, tan 6, 


where: 
” rep 
+ K 


. Mi 
tan 0m = ; 
w+tk 


(The + sign must be used for a positive circularly polarized wave; 
the negative sign for the negative circularly polarized wave. 


: e 
tan 64a = dielectric loss tangent 


= complex dielectric constant 
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It is almost impossible to get a feeling for what these equations mean 
with respect to a wave travelling through the medium, especially since 
mw ind A are vivel by equations which are almost “is difficult to per- 
ceive. An appreciation ol these equations can be obtained however, hy 
reference to hig t which vives qualitatively the behavior predicted by 
these CXPresslons Essentially, @ and 9 are functions of two variables. 
‘These are w, the frequency of the wave, and H, , the applied magnetic 
field. In Fig. 4, the index of refraction and attenuation of the positive 


circularly polarized component are given relative to these values for 





OF REFRACTION ABSORPTION OF 
COMPONENT POSITIVE COMPONENT 


T 


| 





(N-Nz) Mz 


| 
| 
| 


SAMPLE 
AL SYMMETRY 











APPLIED MAGNETIC FIELD (ARBITRARY UNITS) 
ndex of refraction and absorption of a positive eireularly polarized 
to the same quantities for a negative circularly polarized wave 
ited through a» magnetized medium 


the negative component. Hence, both the index of refraction and 
attenuation of the negative component are represented by the abscissa 
of the graph In Fig. 4 these quantities are plotted as a funetion of the 
applied magnetic field for a wave of a fixed frequeney. Many of the 
properties of the medium are clearly displayed in this graph. In partic- 
ular, as the field necessary for ferromagnetic resonance is approached, 
the attenuation of the positive Component becomes larger and larger. 
Kventually this component will be ubstantially completely absorbed 
and only the negative circularly polarized component will be propagated. 
Hence it should be possible to establish a circularly polarized wave in a 
waveguide simply by passing the dominant mode through a ferromag- 
netic material which is subjected to a longitudinal magnetic field of the 
proper amplitude. However, there will be an absorption of one-half 


oft the power being propagated It ig. } had been plotted us a function 
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of the frequency of the wave for a fixed magnetic field, a similar set of 
curves would have resulted. This set would indicate the frequency 
dependence of the Faraday rotation. If the frequency of the wave is far 
removed from the resonance frequency, the difference between the 
indices of refraction of the positive and negative component is not 
frequency dependent. However, near resonance, this difference is a 
very rapidly varying function of the frequency. It is to be remembered 
that these equations were derived for an infinite plane wave. However, 
it would be expected that these equations would describe quite accurately 
the propagation of the dominant mode in a waveguide. The approxima- 
tion would, of course, be better when the cut-off wavelength was much 
greater than the unbounded wavelength. This condition is met when 
the waveguide is filled with ferrite and for these cases quantitative 
agreement is obtained. 

The above analysis shows that if « dominant mode wave (plane 
polarized) is incident upon a ferromagnetic material which is magnetized 
along the length of the waveguide, the wave will split into positive and 
negative circularly polarized waves whose phase constants are given by 
Equation (14).) Since the two circular components travel with different 
velocities in the medium, they will upon emerging from it unite to form 
a plane polarized wave whose plane of polarization has been rotated with 
respect to the incident polarization. The angle of rotation of the polariza- 
tion Is given by: 


where: 
path length through ferromagnetic material (em) 


In order to evaluate Equation (15), it must be combined with Equation 
(14). However, a few approximations are valid in Equation (14) which 
make it much simpler. In particular many ferrites exist for which the 
magnetic losses are extremely small as long as the internal field within 


the body is kept small so that the frequency of the wave does not ap- 


proach the ferromagnetic resonance frequency. This field can be kept 
small if the magnetic field is not raised above the point necessary to 
saturate the ferrite. Kittel has shown that for a finite body the effective 
internal magnetic field that determines the resonant frequency is given 
by: 


H? = (H, + (N. — N,)M,|H. + (N 


4 
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The terromagnetic resonance frequency is 
28H, V1 + oc megacycles (16) 
it the following formula is approximately valid 


easily shown that 
rromagnetic body with a circular or square cross-section 


ler 


1) 
| 


ir + N(u 
H, 
\ 


H 
tor t ¥i 


ration. Where: 


true de permeability at saturation 


demagnetizing factor in 2 and y directions 
1000 is assumed for the de permeability, then 


ivernuge value a] 
© readily computed for various shapes 


i thin dise: 
tor 


H, 
1000 


‘saturates at 1500 gauss. then: 


Lo oersteds 


H, 


~ 4.2 megacycle 


lree ~~ 2.4 


ng thi pel ceil: 


HH LOOOTT,, 


For this case the body could: be saturated with a field of about 1.5 
oersteds 


H, 1500 oersteds 


1200 megacyeles 


the resonance lrequency 1s so close to the operating 
due to lerromagnetic resonance become pro- 


losses 
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hibitive, it is wise to then raise the applied field to some high value, 
so that the resonance frequency will fall well above the operating fre- 
quency. Thus, for many cases of interest it is possible by various means 
to place the ferromagnetic resonance absorption frequency sufficiently 
far from the operating frequency so that magnetic losses due to this 
phenomenon are negligible.* The data accumulated to date indicate 
that the major component of the magnetic losses at microwave fre- 
quencies is due to this phenomenon. Only in a few cases have data been 
taken which have indicated that other factors, such as domain wall 
relaxation, contribute to the magnetic loss at microwave frequencies. 

If then, the magnetic field is controlled so that the ferromagnetic 
resonance absorption is negligible, Equations (13) and (14) can be 
simplified to: 

+ K’)e’ 


{ ‘\ + tan dy (18) 


1+ tan*ég+ 1 


where uv’ and K’ are given in the appendix. 

If Equation (21) is now inserted into Equation (15) a formula for 
rotation is obtained which is valid within the limits of the above ap- 
proximations. If in addition, the frequency of the wave is sufficiently 
greater than the resonance frequency, so that: 


Z “oe 
Wree Kw 22) 


then Equation (15) takes the particularly simple form: 


toM oy toM x (23) 


) 
- WwW Gd 


Most ferrites saturate at 2,000 gauss or less. Hence, for a frequency of 


* This is not always possible, for some ferrites, in the polyerystalline state, 
exhibit extremely broad ferromagnetic resonance absorption lines and it is diffi 


cult to operate at any frequency without appreciable absorption 
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9 OO0 megacycles, 


imMiy — 2000 X 17.6 X 10 
w GO) XK 2r &K 10° 





0.622 (24) 
Hence, the following approximation will be valid to within 5 per cent. 


’ ie Miy _ 4 ta ") 


With this approximation, Hquation (23) reduces to: 


+ 
{ Deo 4 ») 
) 


Equation (25) is quite remarkable. Not only does it- predict large 
| | . 


© [40M] (25) 


rotations, but it also predicts that, within the above approximations 
the rotation will not depend upon the frequency of the incident radia- 
tion. For the assumed values, 


to 1000, 
Mquation 25 predicts rotations of, 


i 
{ 


bo°/em. 


DESCRIPTION OF EQUIPMENT AND MEASURING TECHNIQUES 


The Faraday rotation has been measured in a large number of ferrites 
in order to verify the above theory and in an effort to improve the 
characteristics of the microwave gyrator. A diagram of the experimental 
equipment is given in Fig. 5, and a diagram of the test chamber in which 
the rotations were measured is given in Fig. 6. In the test chamber, 
two rectangular waveguides are separated by a circular waveguide, 


the proper nonreflective transitions being made at each end of the 






circular section, Which is about twelve inches long. One rectangular guide 
is supported so that it can be rotated about the longitudinal axis of the 
system. The dominant 7k) mode is excited in one rectangular guide, 
and by means of the smooth transition this goes over into the dominant 
if ode ty the circular guide. The rectangular guide on the opposite 
ne oowill cecept only that component of the polarization which coincides 
wih the JI ode in that guide, the other component being reflected | 
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at the transition. Absorbing vanes, inserted in the circular section, 
absorb this reflected component. The circular guide is placed in a solenoid 
to establish an axial magnetic field along its length. 

The ferrite cylinders to be measured were placed at the mid-section 
of the circular guide. When a evlinder was used which did not fill the 
cross-section of the guide, it was supported along the axis of the guide 
by means of a hollow polystyrene cylinder which did fill the guide. 

In addition to measuring the Faraday rotation, measurements of 
insertion loss were made by determining the power transmitted under 


identical conditions with the ferrite cylinder removed, and the ellipticity 


of the transmitted wave was determined by measuring the power trans- 
mitted when the rectangular guide on the detector side was rotated to 
both positions of maximum and minimum transmission. Power trans- 
mission measurements could be repeated within 0.2 db. Measurements 
of the angle of rotation of the plane of polarization could be repeated 
within 3° except in the region close to the gyromagnetic resonance where 
rotations were large and ellipticity so great that it was difficult to decide 
the positions of maximum and minimum transmission. These errors 
increased up to the point where the transmitted wave was circularly 
polarized where it was impossible to measure the angle of rotation. 


EXPERIMENTAL RESULTS 


Equation (25) indicates that the rotation per unit path length through 
the ferromagnetic material is proportional to the’ magnetization of the 
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Fig. 5—Experimental equipment set-up used to measure Faraday rotations 
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sample and is not dependent directly on the applied magnetic field 


hig. 7 shows the dependence of rotation upon magnetization for a 
sample of manganese zine ferrite, and indicates that after the sample Is 
saturated, the rotation is sensibly independent of the appled magnetic 
field. fh vuidition, the complex dielectric Constant and the saturation 
magnetization of this sample were measured. From these the rotation 


per centimete! path can be computed from the above theorv using 
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Fig. 7 Angle of rotation versus applied magnetie field for a thin dise of mar 
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Equation (25). For a particular sample of manganese zine ferrite, 
following measurements were made: 
(Sample No. 1) 
17 
4 
bo M .., 1500 gauss 
Using this data, equation (25) predicts: 


4 
f 


121.2°/em 


It is seen in Fig. 6, that the actual measured rotation at saturation 1s 
approximately 123°/em. Hence an extremely good agreement with 
theory has been obtained for this particular sample. 


Equation (25) also indicates that the rotation per unit path length 


should be sensibly independent of frequency within the above ap 
proximations. The data are shown in Fig. 8. However, it will be noticed 


that the frequency difference between these two sets of data is relatively 


small (3 per cent), and the cumulative experimental error in measuring 
angles is such that it is difficult to state that the rotation is closet 
than 1° between the two sets of data. This represents a possible differ 
ence of 5 per cent in the rotation for a change of 3 per cent in the fre 
quency. Thus, even though these preliminary data support Equation 
(25), it cannot be accepted as conclusive evidence until more measure 
ments can be made over a wider band width. 
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Fig. S Dependence of Faraday rotation upon frequency 
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The loss characteristics of different ferrites as a function of the ap- 
plied magnetic field differed distinctly from each other. Some ferrites, 
such as manganese zine ferrite showed extremely high loss which was 
associated with the imaginary part of the dielectric constant. This loss 
was not affected by the application of a magnetic field but remained 
substantially constant as the field was applied. However, as the field 
approached that necessary for ferromagnetic resonance, the total power 
absorbed by the ferrite increased, since the positive circularly polarized 
component Was almost completely absorbed by the sample. In faet by 
measuring the ellipticity of the transmitted wave, it is possible to com- 
pute the difference between the absorption of the positive and negative 
circularly polarized components. This has been done for Sample No. 1 
and the result is indicated in Fig. 9. If the curve were continued to 
higher fields, it would represent the shape of the ferromagnetic resonance 
absorption line 

Some ferrites, such as Ferramic G, showed an almost zero dielectric 
loss but on the other hand caused an extremely large absorption at 9000 
megaeveles due to magnetic losses. The major contributions to magnetic 
loss at this frequeney should be either losses associated with a domain 
wall relaxation or ferromagnetic resonance absorption due to anisotropy 
fields. Unequivocal data can be obtained by the above techniques to 
identify which loss is predominant. If the loss were due to domain wall 


relaxation (or resonance) it would absorb both the negative and positive 


circularly polarized components equally. Thus as the magnetic field 


was applied and as the ferrite became saturated, the losses in both 


components should decrease as the domain walls disappeared. However, 
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Fig. 10—Rotation of the plane of polarization and absorption of the positive 
and negative circularly polarized components when an electromagnetic wave 
(X-band) is propagated through a cylinder of Ferramie G 


if instead the loss is associated with ferromagnetic resonance absorption, 
the absorption of the positive component should begin to increase as 
soon as more domains are lined up in a direetion where they can absorb 
the positive component. Thus even before the sample is saturated the 
absorption of the positive component should be much larger than the 
absorption of the negative component. Of course, in a polycrystalline 
sample with a large anisotropy both components can be absorbed by 
ferromagnetic resonance absorption when the sample is not completely 
saturated, since the random orientation of the domains which occurred in 
zero field has not been completely eliminated until complete saturation 
occurs. Fig. 10 illustrates the rotation per cm path length versus applied 
magnetic field for a sample of Ferramic G. Superimposed on the same 
figure are curves showing the absorption of the negative and positive 
circularly polarized components. It will be noticed that as soon as the 
sample is saturated, the sample becomes transparent to the negative 
component but almost completely absorbs the positive component. 
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Hence the transmitted wave at this point is almost completely circularly 

polarized, even though the applied magnetic field would indicate that 

the resonance absorption frequency was far removed from 9000 me. 
Table I gives the data taken on several ferrites at 9000 me. 


APPLICATIONS OF THE FERROMAGNETIC FARADAY EFFECT-~THE MICRO- 

WAVE GYRATOR 

\s pointed out in the introduction, the haraday rotation affords an 
antl reciprocal phenomenon from which a microwave gvrator can be 
constructed. Such a gvrator is illustrated in Fig. 11 along with diagrams 
vhich help explain its action. Beneath the gyrator are construction lines 
vhich indicate the plane of polarization of a wave as it travels through 
the gvrator in either direction. On each diagram is a dotted sine wave 
Which is for reference purpose only and indicates the constant plane ol 
polarization of an unrotated wave. It is noticed that for propagation 
from left to right in Fig. 11, the screw rotation introduced by the twisted 


rectangular guide adds to the 90° rotation given to the wave by the 


ferrite element making a total rotation of 180°. For a wave travelling 


in the reverse direction, these two rotations cancel each other, producing 
a net zero rotation through the complete element. The unique property 
of the Faraday rotation becomes immediately apparent from this dia- 
gram. In the case of the rotation induced by the twisted rectangular 
guide, the wave rotates in one direction in going from left to right through 
the twisted section, and rotates in the opposite direction when it trans- 
verses the section from right to left. For the case of the rotation induced 
by the ferrite element, the direction of rotation is indicated by the 
arrow in the upper figure for either direction of propagation. The im- 
portant characteristic of the element is the time phase relation between 
two points such as A and B in the upper diagram. It is seen with the help 
of the diagrams illustrating the rotating waves that the field variations 
are in phase at points A and B for propagation from left to mght, and 
they are ISO™ out ol phase for propagation from right to left. In other 
words the transmission line is an integral number of wavelengths long 
hetween A and B for propagation from left to right and is an odd integral 
number of half wavelengths long for propagation from right to left. 
Krom the above deseription of the properties of the gvrator, many 
of its appheations in microwave technology become immediately ap- 
parent. Before discussing these applications in more detail, however, 
it is advantageous to introduce standardized terminology and cireult 
symbols which apply to the gvrator and to other circuit: elements 


derivable from it 
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TABLE I 


Number 
oersteds 
(db 
Hi 


db 


Material Dimension (em 


Field 
pat h 


Applied Magnetic 
Rotation/em 
Insertion Loss 


Sample 
SWR on Input 
Line 


0.447 X 2.28 0 0 10 
length & dia 245 15.610: 
9) 33.510 
735} 58.2) 9.2 
O80 6 9 
1225 7 9.2 
1470 12 10 
1715 l 
1960 l 
2206 ] 
2450) 12: l 
2695 I: 
2940 l 
3185 
3675 


3 


0 
245 
190 
735 
980 
1225 
1715 
1960 
2450 


Ferramic A 2.5 0 
245 

190 

735 

GSO 

1225 

1715 

1960 

2450 


2695 


Ferramic G 23 0 : >30) 
245° 38 s/21. 23 
190 77 7 
735 124 : 2 
980 157 ( ] 
1225 170 0. 
1470 180 0 
3430 ¢.p. 0. 


_* Data given is the difference in db between the major and minor components 
of the elliptically polarized transmitted wave 
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The “active” element of the device, the ferrite evlinder, has been 


66 ” 


termed a ‘Faraday Plate. 

As was pointed out earlier, the fundamental property of the gyrator 
is the 180° phase difference introduced between the two directions of 
propagation through it. Thus the gyrator may be thought of as a four 
terminal circuit element having no phase shift for one direction of trans- 
mission, and having a 180° phase shift for the opposite direction of 
transmission. A convenient circuit symbol for the gvrator, which indi- 
cates this property, is shown in Fig. 12. 

If the rectangular waveguides on each side of the Faraday Plate are 


rotated about their common axis so as to make an angle of 45° with 


[ — | 


--— 7 [a 


Fig. 12—-Cireuit symbol for gyrator 
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CIRCUIT SYMBOL FOR 
POLARIZATION CIRCULATOR POLARIZATION CIRCULATOR 


Fig. 13--Schematie diagram of polarization circulator. 


each other, then a one-way transmission system can be created which 
is similar to Lord Rayleigh’s one-way transmission system of optics, but 
with the important difference that this one-way transmission system 
does not depend upon frequency but is broad band. This one-way trans- 
mission system can be used, for example, to isolate the generator o1 
detector from the waveguide in microwave systems. In this application 
it has the great advantage over the attenuators which are presently 
used for this purpose in that it can be made practically lossless for the 
direction of propagation which is desired but the reflected wave will be 
completely absorbed and hence more complete isolation can be effected. 

A more complex and more useful circuit element, than this simple 
one-way transmission property would at first indicate, is obtained by 
adding a second connection on each side of the 45° Faraday Plate. It is 


suggested that this device be called a polarization circulator. Thus, the 
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ion circulator actually has four output branches corresponding 
vo different polarizations at each end. The polarizations of the 
tout branches are indicated in Fig. 18. [t is noticed that powe1 


Into the polarization eireulator with polarization a Is turned imto 


polarization b, also 6 is turned into ¢, € Is turned into d, and d is turned 


into minus a. This property is indicated very clearly by the circuit sym 


bol suggested in ‘Fig. 13, the phase inversion between arms d and a being 


wo 
j 


“| by the minus sign between the d and a arms. 
Another one-way transmission system can be created by combining 
the gyrator with two-normal hybrids. This combination is indicated in 


Vig. 14. Sinee this device has all of the fundamental properties of the 
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Fig. 14 Schematie diagram of circulator. 


polarization circulator with the exception of the phase inversion be- 
tween arms d and a it is suggested that it be called a ‘‘cireulator’’ and 
the cireuit symbol suggested which indicates its properties is also given 
in big. 14 

Phis list of applications is obviously not complete since it includes 
only the fundamental elements from which innumerable specific applica- 
tions can be made. 

In addition to the applications discussed above, which depend upon 
the anti-reciprocal property of the element for their operation there 
are several simple applications which are based only upon the fact that 
the amount of rotation can be controlled externally by adjusting the 
magnetic field. Among these uses are electrically controlled attenuators, 


modulators, and microwave switches 
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\ PPENDIX 

The equation of motion of the magnetization of a ferromagnetic 
material is: 

OM 

ot 


=» = x B= 7 iw x (M x HD) 


where 
HT = internal magnetic field (oersteds) 
tr.JJ = magnetization of medium (gauss) 


a = parameter which measures the magnitude of the damping force 

on the precessing dipole moment of the sample 

Y = gyromagnetic ratio of the electron (7 ge/2me where g is the 

Landé g factor for the eleetron). 

If a ferromagnetic material is subjected to a steady magnetic field, 
H,, along the z axis and if then an alternating field is applied in an 
arbitrary direction, Equation (1) must be solved in order to find the 
behavior of the magnetization of the material. To solve this problem, 
the following notation is introduced: 

lr. = magnetization of medium in absence of alternating field 

externally applied steady magnetic field (oersteds) 

Rs, Res & components of applied alternating magnetic field 

m,, m,, m. = alternating components of magnetization 

A Me 
h; = hz — Nem-z 

P 
h; h, — Nym, 


H, = H, +h. VAM. + m,) 


H; = components of internal magnetic field 


demagnetizing factors of body. 
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Hence the magnetic field, 1, occurring in Equation (1) is defined by: 


Bw hk 


V mt + Me? + 


In solving Equation (1), an exponential, exp [jot], time dependence 


is assumed for the alternating magnetic field and magnetization, and 
if the following assumption is made: 


ee Me ao 


it is easily shown that the alternating Components of the magnetization 


of the medium are given by (neglecting terms of the second order in 


small quantities 
ly MHL + @&) + JyaM wih, — jyM.wh, 
¥y I] 1 + a’) w” T }[2wyall:| 
Y VHA 4 a’) rT Jy¥a VW wlh, T JyM wh, 
VHY A + a’) — wo + jwyall?] 


where: 


since 


(3) 


it is possible by means of Equations (2) and (3) to find the 


relation 
between the alternating flux density b and the internal alternating field 
h’. Uf the ferromagnetic body is considered as being infinite, the internal 
fields and applied fields are equal. Hence, for this case: 


h uh qkh 


Kh, t uh, 


4 


h, 
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[ly Hil + a’) = w'|[4e My H,(1 + 0°)] + 84M oye, 
=I+ <= — — 
[lyHi(l + a*) — w|? + 4w*y"a°Hi 
4nM yoly Hil + a’) — «| 
ly He + a’) — w' |” r lw ya HH: 
82M .y7aH 
~ [YHA + a’) — w’|? + 4o*y'a°H, 
‘. 4M yoauly'Hi(1 + a’) + w’| 
[VHG(L + a’) — w') + 40*y*a" Hi; 

In order to find the behavior of a wave being propagated in this 
medium, it is necessary to find a solution to Maxwell’s equations which 
are consistent with the above set of equations and in which, b, h, EF, 
and D are of the following form: 

b = bo exp [ jot — r(n-r)] 
h = ho exp [jot — I'(n-r)] 
E = Ey exp [jwt -- '(n-r)] 
D = Dy exp [jolt - r(n-r)] 
where Ey, and ho are complex vector functions of the coordinates and 
which satisfy the boundary conditions imposed by the waveguide. 
Further: 
n = unit vector in the direction of propagation 
I’ = propagation constant 


Maxwell’s equations are: 


1 ab 


i = 
Vx cal 


| aD 


VXh=--— 
c at 


Inserting the values given in Equations (5), these become: 
fel 


vx E —- Tt X B) = 


( 


VX hy — THX hy) = ee 


Cc 
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which can be combined to 
vx (Vv &X h rn x ho) rn x (0 K ho — In XK ho) 


g Equation (9) in component form gives: 


ah y 0 h, Oh, Oh, . Ohy ‘ Oh. 
. - I Ne — +t 
OYOr OZ02 oy Zz Oy Oz 
oh z Oh, . a Oh ) 


OY 42 LC ee Oa 


oh D Ns oh y Oh, . Oh, Oh, 
—, In, i - + — 
OYOR 2 Or Oz 


Ohy : Oh, : Oh; 
_ rin, —+n,- 
OZ Oy OY 


yhy 4 nzhz) Ih, 


Oh, Oh, 


oy” 


Oh. , Oh, 
’ T My = 
Ox OY 


) + In,(nehz + Nyhy +- nh) 


where the subscript 0 has been dropped from all components for con- 
yvermence 
[i the wave in question is an infinite plane wave being propagated 


ilong the z axis, then: 


n 0, n | 


7] 


and the components of jo are constants, and h, is zero. For this particular 


CAN Mquations 10), (11) and (12) become: 


Ivh, 
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Equations (13) and (14) are general differential equations derivable 
from Maxwell’s equations and do not yet contain the properties of any 
particular medium. In order to find the behavior of a wave travelling 
through an infinite ferromagnetic medium which is magnetized along the 
direction of propagation, it is necessary to combine these equations with 
Equations (4) which describe the relation between } and hin the medium. 


This gives: 


(uh, jKh,) ne - = I'l 


2 
eo we 2 . 
(uh, + 7Khz) — = —I“h, 16) 
a 
The only possible solution to this set of equations is a circularly 
polarized wave where: 

Ry. = t jh, 
The positive sign above represents a so-called positive circularly 
polarized wave and the negative sign a negative circularly polarized 


wave. The propagation constants for these waves is given by 


r, =" Ve(u + K) (17) 


REFERENCES 


. B.D. H. Tellegen, Philips Research Reports, 3, pp. 81-101 (1948); 3, pp. 321 
337 (1948); 4, pp. 31-37 (1949); 4, pp. 366-369 (1949) 
A. G. Fox, unpublished memoranda, Bell Telephone Laboratories 
I. M. MeMillan, J. Acous. Soc. of Am., 18, pp. 344-347 (1946 
. J. W. Miles, J. Acous. Soc of Am., 19, pp. 910-913 (1947). 
F. Bloch, Phys. Rev., 70, pp. 460-485 (1946). 
». H. G. Beljers and J. L. Snoek, ’’Gyromagnetic Phenomena Occurring Within 
Ferrites’. Philips Tech. Rev., 11, May, 1950, pp. 313-322 
. E.M. MeMillan, J. Acous. Soc. of Am., 19, p. 922 (1947) and H. B. G. Casimir, 
Nuovo Cimento, Vol. VI, Series IX (1949) 
W. P. Mason, unpublished memoranda, Bell Telephone Laboatories 
. Rayleigh, Nature, Vol. 64. 
. H. Konig, Optik, 3, pp. 101-119 (1948) 
D. Polder, Phil. Mag., 40, pp. 99-115 (1949 
2. W. A. Yager, J. K. Galt, et al. Phys. Rev., 80, pp. 744-748 (1950). 





Dialing Habits of Telephone Customers 


BY CHARLES CLOS AND ROGER I. WILKINSON 


Manuseript received October 3, 1951 


Thi paper considers the heharvior o} CuUSIOMETS waiting lo dial calls, when 
dial tone is delayed. Tests were made in a panel dial central office, from 
which were determined: relationship between load carried by a group of 
line finders and the resultant dial tone delay; measures, by classes of ser- 
vice, of the magnitude of the generalized trunking formula’s “7” factor 
describing the degree to which customers wait when dial tone is delayed; 
COMpariSONs of observed and theoretical distributions of the number of 
simultaneous calls on line finder groups; and statistical accounts of the 


} 
achior of customers when dial tone is delayed, 


Following World War IT the conversion of great quantities of manual 
telephone equipment to dial, and the addition of large numbers of new 
telephones, mostly dial, in the Bell System has directed increasing atten- 
tion to those service problems peculiar to automatic operation. These 
problems concern chiefly the provision ol adequate amounts of equip- 
ment to give satisfactory service at all times. One of the important 
factors affecting the amount of equipment needed is the action of the 
customers themselves when their calls are momentarily blocked due to 
these equipment shortages. The actions of subscribers whose calls are 
blocked due to a shortage of trunk equipment have been reported pre- 
viously.! ‘Phis paper considers the behavior of subscribers, waiting to 
dial calls, when dial tone is delayed. 

During 1949, Bell Telephone Laboratories conducted a series of tests 
at the New York Telephone Company’s Sterling-3 panel dial central 
office in Brooklyn, N. Y., with the object of increasing the knowledge 
available regarding subscribers’ actions and their effects when dial tone 
is delayed 

Che following principal results were obtained from the Sterling-3 tests: 

|. The relationship between the load carried by a group of line finders 
and the resultant dial tone delay 


2.) Measures, by classes of service, of the magnitude of the general- 


rles Clos, ‘An Aspect of the Dialing Behavior of Subseribers and Its 
the Trunk Plant.”’ Bell Sustem Tech. J., 27, July 1948 
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ized trunking formula’s “‘j” factor describing the degree to which cus- 
tomers wait when dial tone is not immediate.” 

3. Comparisons of observed and theoretical distributions of the num- 
bers of simultaneous calls on the line finder groups. 

t. Statistical accounts of the actions of subseribers when dial tone 
is delayed. 


GENERAL PLAN OF THE TESTS 


The general plan of the tests was developed around two specially 
constructed devices: 

1. A 100-pen recorder capable of recording observations continuously 
for two hours and sensitive enough to record individual dial pulses 

2. A speed of dial tone measuring set comprising a means for manu- 
ally originating test calls, and an electric timer which automatically 
stopped when dial tone was received. 

Tests were conducted on a weekly schedule, one or two line finder 
groups being studied each week. Two 100-pen recorder tapes were run 
daily. The first was run from 10 A.M. to noon for message rate individ- 
ual, flat rate individual and message rate two-party classes of service to 
include the morning busy periods for these customers. For coin customers 


the first tape was run for two hours during the noon coin busy period. 


The second tape was run in the afternoon from 2 to 4 P.M. for all 
classes of service except coin. For the coin class the second tape was run 
for two hours during the early evening coin busy period. Three message 
rate two-party tapes were run inthe early evening busy period, and two 
coin tapes were run in the afternoons when World Series baseball games 
were being played in Brooklyn. A summary of the number of line finder 
groups observed, the number of tapes taken, the number of line finders 
made available and the maximum per cent dial tone delays over three 
seconds are given in Table I. 

Except for the morning runs on the message rate individual line 
groups where an effort was made to maintain a fixed number of twenty 
line finders throughout, the number of line finders made available was 
selected by close observation of the flow of traffic. All studies were by 
half hours during which the number of line finders was held constant 
as far as possible. At the end of a half-hour period the number of line 

2 This formula for both finite and infinite sources was developed by R. I. Wil 
kinson in 1930, and appeared in the 1936 Bell Telephone Laboratories Out-of-Hour 
Course ‘‘The Theory of Probability as Applied to Telephone Trunking Prob 


lems.’’ This formula for infinite sources was also developed by Conny Palm and 
appeared in ‘Etude des delais d’attente’’ in Erickson Technics——No. 2—1937 
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finders left in service was adjusted in order to obtain a reasonable num- 
ber of dial tone delays in the next half hour without producing a severe 
reaction from the subscribers served by the line groups under study. 

The data recorded on the tapes showed continuously the busy or idle 
conditions of certain circuits associated with the line groups under 
study. In some cases the receipt of dial pulses and the operation of 
registers were also recorded. These circuits included line finders, a few 


subscriber lines, trip circuit sub-groups, the all trunks busy register, 


the peg eount registel and the speed of dial tone measuring device. In 


TABLE | 

Nissan al Maximum 
I sag nies Percent 

a. | pies One 


5 OE conds 


Message rate individual 
Morning tapes 
Afternoon t pes 

Mes ige rate two-part 
Morning tapes 
Afternoon t Lpes 
Kvening tapes 

Flat rate individual 
Morning tapes 
Afternoon tapes 

Com 
Morning tapes 
Afternoon tapes 


] vening t Lp bt 


addition the busy and idle conditions of a sample of senders was ob- 


ved in order to note the general load level on the senders. 


RELATIONSHIP BETWEEN LOAD CARRIED AND PER CENT DIAL TONE DELAY 


One of the principal objectives of these tests was to establish as far 
as possible the relationship between the average load carried by a line- 
finder group and the corresponding dial tone service when there is a 
shortage of line finders but not of senders.? The average load carried 
was obtained by making a switch count every thirty seconds of the 
number of line finders busy as indicated by the 100-pen recorder tape. 
The dial tone tests were made with the speed of dial tone measuring 
device. Forty-five dial tone tests were made each half hour for each 

The restriction of avoiding dial-tone delays due to a sender shortage was to 


eliminate a factor external to the line-finder group and to make the results of the 
ts applicable to both common-control and non-common-control type dial 


t 
tes 


SYSteMmsS 
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line group studied. Additional dial tone tests were made on all other 
line groups in the office as a check that the delays experienced on the 
line groups under study were not due to a sender shortage. The sender 
data on the tapes were also used for this purpose. Figs. l(a) to 1(d), 
and 2(a) to 2(d), inclusive, show for various amounts of load carried, 
the per cent of dial tone tests encountering delays greater than three 
and greater than ten seconds for half-hour study periods for the most 
frequent number of line finders in the tests for each class of service 

Plotted on each of these figures is a theoretical fitting dial tone teste 
delay curve computed for the indicated dial tone delays and for the 
following 7 factor values in the generalized trunking formula, determined 
in a manner to be explained later: 


Class of Service j factor 


MRI 6.6 
MR 2-party 5.8 
FRI 6.5 


Coin 2.1 


To indicate the effect of varying j, several curves have been added 
to Fig. 1(a). Selections of curves for 7 = 0, 1 and © (which correspond 
to the three commonly used infinite source congestion formulae, Erlang 
C, Poisson and Erlang B when adapted to the tester’s delay problem) 
are shown. It is clear that with the wide differences in delays which they 
give for specified loads carried, it is highly desirable to select that 7 
formula for engineering use which most nearly describes the customer 
actions in any situation being dealt with. In the field of curves shown 
on Fig. l(a), the one labelled 7 = 6.6 was derived in a logical manner 


from the data, and shows an agreeably satisfactory fit. For example, 


during a heavy load period when, say, 20 per cent of a dial tone tester’s 
calls are meeting delays greater than 3 seconds, an actual average load 
of 16.6 erlangs (as shown by the 7 = 6.6 curve) would likely be carried. 
(Load in erlangs equals average number of simultaneous calls.) The 
Erlang C (y 0) and Poisson (7 = 1) theories would indicate the 
presence of loads of 15.6 and 16.0 erlangs, respectively, figures clearly 
too small for the circumstances shown by the data of Fig. I(a). On the 
other hand, use of the Erlang B (j = ©) theory would predict a con- 
siderably larger load carried, about 17.9 erlangs, than one would prob- 
ably be justified in assuming here for engineering purposes. 

By grouping the dial tone delay data by bands of load carried, rela- 
tionships of per cent of test calls encountering varying dial tone delays 
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up to twelve seconds were obtained. These data are shown on Figs, 3 
to 6. Fig. 3 is for the message rate individual class of service with twenty 
line finders available. The data on this figure correspond to those on 
Figs. Ifa) and I(b). Fig. 4 is for the message rate two-party class of 
service with ten line finders and corresponds to bigs. l(c) and I(d). 
big. 5 is for the flat rate class of service with ten line finders and cor- 
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responds to Figs. 2(a) and 2(b). Fig. 6 is for the coin class of service 


with 34 line finders and corresponds to Figs. 2(¢) and 2(d). 


Plotted on Figs. 3 to 6 are theoretical fitting dial tone tester delay 


curves, curves A, determined by means of the following formulae: 
1. The generalized trunking formula for determining the proportion 
of calls that encounter congestion, i.e., find all line finders busy. 
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Some of the theoretical aspects of these formulae are considered in 
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the Wo sections that follow 


GENERALIZED TRUNKING FORMULA 
generalized trunking formula combines in one expression the 
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— 


by John Riordan paralleling that of Conny Palm’s 
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trunking formulae and a large field of intermediate assumptions regard- 
ing the disposition of calls which fail to obtain immediate service. These 
assumptions are given in Table IT on page 42. 

One method for developing the generalized trunking formula® is to 
consider the probabilities of the existence of certain states and to de- 
termine the number of transitions during some convenient interval of 


5 In this article only the unlimited sources trunking formula is considered 
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time from one state to anothe By equating certain of these transitions, 
a series of simultaneous equations evolves, which when solved \ ields one 
overall expression. Of interest is the development of the transition equa- 

hus for a case of ¢ line finders arranged in a simple group, let 


r) represent the probability that a where w < oc) line finders are 


occupied and fir + 1) represent the probability that 2 + 1 line finders 


are occupied, Let m be the average number of calls offered to the line 
finders during a long interval of time, 7. Let T be the unit of time and h 


be the average holding time per call measured in terms of 7. Over a 
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period T the number of transits from state x to state x + | must equal 
(or differ by no more than one) the number of transits in the reverse 
direction. That is: 
, +i, 
nf(x) = flx 1) (1) 
h 
It may be noted that nh a, where a is the average offered traffic 
load in erlangs. 
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TABLE II 


Disposa! of Calls that do not Obtain 
ne Finder Immediately 


W iting calls are cleared out at a rate } times the rate at 
which calls are terminated when served by the line 
finders 

Calls are cleared out of the svstem immediately, that is 
no calls wait GQ x 

W titing calls are cleared out at a rate equal to that with 
which calls are terminated when served by the line 
finders (j l 

Calls wait until served (J 0) 


When « 2 ¢ a new situation is encountered. “c”’ calls are engaged in 
‘conversation and a ¢ calls are waiting for service. If the waiting calls 
are forced to wait for an unduly long period of time so that in effect 
they are being denied service, it can be expected that they will wait for 
some average period, say /7, and then abandon their attempts. On this 


basis the corresponding equation is: 


nf (ar) J\z + , f(x (2) 


It has been assumed in the above equations that the distribution of 
the holding times is exponential, an assumption which is found in most 
local systems to be reasonably justified. The distribution of the waiting 
times is also taken to be exponential. By introducing a factor 7, where 
j h H, equation (2) can be written in the simpler form: 


aft a le + ; = r (a (3) 


Solving this system 


whe mw “ ( 


f(x) 
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The probability of a call encountering congestion, which is equivalent 
to the probability of a call having a delay greater than zero units of 
time is: 


r= 


P(>0) = > f(z) 


z=c 


DELAY FORMULA FOR THE DIAL TONE TESTER 


The probability that a dial tone test call encounters congestion is 
given by expression (7). Once a test call has encountered congestion it 
will experience a delay depending upon a number of variables. The 
assumptions underlying the dial tone tester formula are: 

1. A dial tone test call when encountering a delay waits until served. 

2. A dial tone test call does not add to the load offered and carried 
by the line finders. 

3. Upon encountering a delay, a dial tone test call is served in the 
order of its arrival with respect to all other waiting calls. For example, 
if the test call finds three other calls waiting, it waits fourth in line. 

Under the third assumption as calls drop out, due to conversations 
terminating on the occupied line finders or due to waiting calls aban- 
doning their requests for service, the test call advances from an initial 


position of say fourth in line to third in line, then to second, then to 
first in line, and finally is served. The overall delay distribution of the 
test calls depends therefore upon the number of calls they find waiting 
ahead of them. The delay distribution for each such number must be 


weighted by the probability of its occurrence in order to obtain the 
overall distribution. The delay distribution for a test call which finds 
zero calls waiting is: 


Po( >t) = exp (—et jH) (S ) 


The probability is f(c) that a call made at random will find all line 
finders busy with no calls waiting. Hence the weighted delay distribution 
Po(>1), 18: 


P (>t) = fle) pol >t) f(c) exp (—et jl (9 


The delay distribution for a test call which finds one call waiting 
ahead of it is: 


pil >t) — [| + c/j — (e dD exp | -f H)} exp | ct jH) (10) 


The probability is f(e + 1) that a call made at random will find all 
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line finders busy and one call waiting. Hence the weighted delay distri- 
bution, /,(> 1), is: 
P > | } + J C/7) exp ( tH } exp | ct JM) (11) 
Is given by the following formula: 
Fi4i(t) exp (t/H) (12) 


where / 1({f) IS given by ( ONNY Palm.® The over-all delay distribution 


IS then 


P(> Po(>t) + Pi(>d + PA >) + -:- (13) 


By making appropriate substitutions and summing the result, ex- 
| 


pression (14) becomes: 


e+; e+ je + 2%) 


[ , , 
exp | —cl/jH + (a/j)\l exp (—t/H)| (14) 
Iixpression (14) is equivalent to that of Riordan involving two in- 


complete gamma functions as follows: 


y[c/j, (a/j) exp (—t/H)] 


P(>0) : 
yic J» a J) 


(15) 
where the incomplete gamma function, 
y(N, x) x xz (16) 


The theoretical dial tone tester delay curves shown on Figs. I(a) to 
lid), 2a) to 2d), and 3 to 6 were computed from expression (14), 
using the following values of 7 and H for the classes of service studied, 


these values being determined in a manner explained later: 


seconds 
2 seconds 
7 seconds 


seconds 


On Bigs. bia), 1(e), 2(a), and 2(¢), which show the per cent of dial tone 
tests encountering delays greater than three seconds for various amounts 


ition 53, loc. cit 
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of load carried, it may be noted that most of the theoretical dial tone 
tester delay curves are in close agreement with the observed data, with 
a tendency perhaps to be slightly high. On Figs. 1(b), 1(d), 2(b), and 
2(d), which are for dial tone delays greater than ten seconds, it may be 
noted that the theoretical curves have a slightly stronger tendency to 
lie on the high side of the observed data. On Figs. 3 to 6 the theoretical 
dial tone tester delay, curves A, again lie in the proximity of the curves 
of the observed data, with a tendency to lie higher than these latter 
curves, especially at the ends where the dial tone delays are greatest. 
Among the factors which account for this discrepancy are: 

1. A feature is present in panel line finder circuits for momentarily 
releasing trip circuits with waiting calls to prevent the orphaning of 
calls under certain trouble conditions. The release occurs after a call 
has been waiting from 5 to 12 seconds and reoccurs every 7 seconds 
thereafter. When such a release occurs the call yields whatever waiting 
preference it may have had to a subsequently placed call which is not 
yet affected by such a release. The dial tone test calls did not wait 
beyond 12 seconds. Hence for these test calls there was only one possi- 
bility of such a release and for many of them the release occurred near 
the end of their waiting period. Hence they were more likely to gain 
preference over other calls than to lose their preference. 

2. Subscribers while waiting for dial tone frequently become impatient 
and proceed to flash (move their switchhook up and down). While 
flashing, a subscriber may lose preference to a subsequently placed test 
call (the latter of course does not flash). 

3. Many subscribers fail to observe dial tone and proceed to dial. 
During such dialing, a subscriber may lose preference to a subsequently 
placed test eall. 

t. Line finders serve a large proportion of call attempts of short 
holding time whose presence may militate against the occurrence of 
the longer delays. In connection with the measurement of the j factor, 
the following proportions of call attempts and average holding times 
were noted on which no dialing occurred or where no more than two 
digits were dialed. 


Class of Service 


Proportion of Attempts with No 
yunts 


Average Holding Time 


eg 


MRI 35.2% .4 seconds* 
MR 2-Pty 33.3% , -conds* 
FRI 25 // o. seconds 
Coin f / ‘ ‘conds 


* 


Partly estimated 
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The individual contributions of these four factors to discrepancies 
between theory and observation are not easy to assess. The first three 
explain a tendency for test calls to get ahead of calls already waiting 
for dial tone. On Figs. 3 to 6, inclusive, additional theoretical dial 
tone delay curves, curves B, for the case where a dial tone tester 
ilways gets first in line are shown. Even these curves tend to lie above 
the curves of the observed data on Figs. 4 and 5 where ten line finders 
were available; they more nearly agree with the observed data on Fig. 3 
vhere twenty line finders were available, and they lie below the observed 
data on Fig. 6 where 34 line finders were available. This is an indication 
of the fact that with higher traffic loads (which oceurred on the larger 
line finder groups) a test call will encounter more competition from 
other calls and therefore will have a lesser chance of gaining precedence 
over all of the other calls. The fourth factor indicates that the call 
attempts served on line finders consist of two distinet holding time 
universes and not just one, as was assumed in the development of the 
dial tone tester formula. The effeet of the presence of both a short and 
long holding time universe of calls would be to introduce a change of 
slope in the delay curves which may be seen in Figs. 3 to 6 to be at about 
t seconds. There is reason to believe that the same cause may have 
Deen responsible for the tendency of the observed delay curves to fall 
away from the theoretical at the lower levels of load carried. 

Due to the reasons given above and to the fact that the dial tone 
delay observations were made by the test call method, the above re- 


sults may not directly describe service from the customer’s point of 


view. Conny Palm has developed the following formula which gives a 
slightly different measure of customers’ dial tone service. It indicates 
the proportion of calls which have neither received dial tone nor have 
dropped out at time f 


vle/), (a ) exp (—t/H)) ad 
yic/. /) | -exp (—t/H) (17) 
¥(c/j, @/7) 


P(> P(>QO) 
Curves for this formula are shown plotted on Fig. I(a) and at C on 
Figs. 3 to 6. They are quite close in many cases to the observed dial 
tone tester results. It would appear that a sufficiently good estimate 
of the customer's dial tone service, whatever its precise definition, can 
be obtained by the dial tone tester method. 
Recently revised tables for the capacity of step-by-step line finders 
have been published for Bell System use based on Palm’s formula using 


5. This was selected as being slightly conservative for 
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most applications after reviewing the above Sterling-3 results and other 
line finder data collected in step-by-step offices. 


MEASUREMENT OF THE ) FACTOR BY CLASSES OF SERVICE 


As indicated previously, the data recorded on the tapes showed the 
states of being busy or idle and of changes in these states for line finders 
and the associated trip circuits. A fully equipped line finder group of 
100 lines has ten trip circuits each of which serves two sub-groups of 
twenty subscriber lines in the following manner. When a line originates 
a call its line relay is operated. This causes a ground to appear on a 
lead which is common to all twenty line relays in the sub-group and 
starts a line finder hunting for the calling subscriber’s line. As soon as 
this hunt is completed the cutoff relay associated with the calling line 
operates and disconnects the line relay, removing the ground (unless, 
of course, another line in the sub-group has originated a call in the 
meantime). During periods of overload when line finders are not im- 
mediately available, the ground due to a single subscriber will persist 
until: 

1. A line finder is obtained, or 

2. The subscriber abandons the attempt, or 

3. The subscriber receives an incoming call which operates the cut 
off relay. 

The twenty leads from the trip circuit sub-groups were brought out 
to the pen recorder and a record taken of the grounds that occurred on , 
each lead. Except for the possibility that more than one subscriber is 
waiting for service at the same time on a given trip circuit sub-group, 
the record of the occurrences of the grounds gives a substantially ac- 
curate’ record of the demands for service and of the number of calls 
waiting for service. Hence an analysis of the events occurring on the 
trip circuit sub-groups and on the line finders as recorded on the tapes 
gives a means for determining H. The quantity H was introduced in 


equation (2) in the term 


(xr + 1 — c) al 
—— f(r + | (18) 
rr = 


For convenience in the ensuing discussion this term will be replaced by 
7 To obtain absolute accuracy would require the use of a pen recorder with one 


pen for each of the 400 subscribers served on a line finder group plus one for each 
line finder. 
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the equi lent expression : 
(19) 


The average number of waiting calls that drop out per unit 
of time during the state (z + 

The number of waiting calls 

The number of line finders occupied with calls. 

\ measure of the rate at which calls tend to abandon 
Walting 


The proportion of time that the state (z +4 y) exists. 


(on the tapes we Cah measure f 2 + y) and count Nw Hence H can 
be determined. The result is a statistical quantity subject to many 
chance factors. In the actual analysis of a tape, the composite average 
value of HZ was determined for all possible observed states where calls 
were waiting. By an analogous process, the composite average value 
of A for all possible observed states where calls were being served by 
line finders was determined. Also as a side computation, a composite 
iverage value of A’ for calls that were served by line finders but for 
Which no peg counts were scored was determined. This value of h’ is 
included in A on the basis that data for engineering line finders consist 


of estimated calls based on peg counts and of holding times which in- 


clude an allowance for these short holding time calls. The average 


values of /7, A and of the ) factor for the four classes of service studied 
are given in ‘Table IT] 

The results for /7 by individual half hours and by various percentages 
of dial tone delays greater than three seconds are shown on Figs. 7(a) 
to 7(d) respectively for the tour classes of service. On some of these 
figures an upward bulge may be noted in the center. This is not con- 
sidered to be characteristic of the habits of the subseribers but is the 
overall effeet resulting from a number of arbitrary rules followed in 


making the analysis in order to simplify the work and to offset par- 


TaBie III 


te individual 





DIALING HABITS OF TELEPHONE CUSTOMERS 





MESSAGE RATE- 
NDOIVIOUAL 


. AVERAGE . 
| ay een On eee a ai ‘ AVERAGE 
; doogpepocodcccegecs 4 














MESSAGE RATE- 
TWO-PARTY 


AVERAGE | 


(b) 


30 40 50 60 70 80 90 0 10 20 30 40 
PER CENT DIAL TONE DELAYS GREATER THAN 3 SECONDS 














AVERAGE CUSTOMER WAITING TIME 
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tially the effect of occasionally having two or more calls waiting on one 
trip circuit sub-group. The rules and the reasons for them will be de- 
scribed with the aid of Fig. 8, which shows a hypothetical section of 


one of the tapes. The rules were as follows: 


1. Initial Overlap 


Referring to Fig. 8, at 4; a subseriber has initiated a request for 


service. At fo a line finder rises to serve the subscriber. At ts the sub- 
seriber recelves service. This case is typical ol a subseriber receiving 
prompt dial tone service. 

The span from f; to fg was difficult to measure accurately because, 
for the usual case, it was about the same as the maximum error due to 
misalignment of the recorder pens. It was not measured unless the com- 
bined span from f; to ¢; exceeded one second. 

The span from f, to ¢; involves an overlap, it represents a period 
when a line finder is busy hunting for the terminal of the subseriber who 
originated the request for service. It also represents a period when a 
subscriber is waiting for service. In the analysis this span was treated 
as a case where a line finder was busy with a call and not as a call wait- 
ing for service. 

If the span from f) to ¢; and all similar cases had been treated as 
calls waiting for service and if in addition all spans from ¢; to ft which 
were not measured had also been treated as calls waiting for service, 
the average values for /1 would have increased slightly for each class 


of service. 
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4 1 hree Second Rul Jor the Bridging of Calls 


Referring to Fig. 8, again, at 4; a request for service is originated on 
trip circuit 5 and at ft this request is withdrawn. At to apparently a 
new request for service is initiated which is then withdrawn at ft). From 
manual service observations it is found that subscribers often flash when 
dial tone is slow. A few pens were used to observe individual subseribers, 
and big + shows a case where a subseriber made several flashes when 


his tone was slow. When a subseriber flashes it appears as though he 


ee er ET Pe rae Pe arte ee ene ae 






































t7 te to ty ti 
> 
hig. S—Section of a hypothetical tape showing activities on trip circuit sub- 


groups and on line finders 


were making several bids for service. Actually he is making only one 


real bid. On the trip circuit sub-group pens it was generally impossible 


to distinguish between flashes and requests for service by two or more 
subseribers. To resolve this problem many observations of subseriber 
lines recorded on the tapes were examined from which it was concluded 
that no great error would result if a break in the demand for service 
on a trip circuit sub-group of less than three seconds were considered 
asa flash and was to be bridged, and a break greater than three seconds 
was to be considered as the termination of one call attempt and the 
start of another 
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S. Treatment of Case s Whe re Two or Vore Calls We re Found lo he 
Waiting on One Trip Circuit Sub-Group 


The occurrence of several calls waiting on one trip cireuit was oecasion- 
ally noted in the analysis. Referring to Fig. 8, a case is shown on trip 
cireult sub-group 12. At fy a line finder is seized. Trip circuit sub-group 
IS shows that a subscriber is dialing before tone. The appearance of 
dial pulses on this trip circuit indicates that only one subseriber is 


demanding service otherwise the dialing would not show. Trip circuit 


sub-group 12 however appears to have two or more requests for service. 
One of these requests for service began at t;. The start of the second 
request occurred somewhere between ¢, and fg, perhaps half-way be- 
tween. At fy, one of the requests was served by a line finder. To sim- 
plify the handling of such cases, the assumption was made that the 
first attempt started at ¢; and ended at ty and the second request started 


SUBSCRIBER FLASHED 


, 








4 
LINE RELAY OPERATED LINE RELAY RELEASED 
(SUBSCRIBER TOOK (SUBSCRIBER DEPRESSED 
RECEIVER OFF HOOK) SWITCHHOOK) 


21.6 SECONDS 


Fig. 9—Example of a customer flashing for dial tone (Tape made October 18). 


at ty. The effect of this is to understate by an indeterminate amount 
the average value of H for each class of service. This understatement 
should be noticeable for the higher degrees of overload because the 
occurrence of several calls on one trip circuit sub-group is then most 
likely to occur. 

The effeet of several calls simultaneously waiting in a trip circuit 
sub-group and of one or more calls dropping out is to overstate the 
magnitude of H. For instance if two simultaneous call attempts of five 
seconds each overlap for one second and both attempts are abandoned, 
the apparent average waiting period is nine seconds, whereas it should 
be five seconds. It is believed that the above three rules tend to create 
understatements which roughly balance this type of overstatement. 


DISTRIBUTION CURVES OF SIMULTANEOUS CALLS 


The detained analysis of the tapes provided distributions of simul- 
taneous calls. For each class of service studied these distributions can 
be compared with theoretical distributions derived from the generalized 
trunking formula using the j factors developed in the analysis. Several 
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such comparisons are shown on Figs. 10 and Il. The agreement is quite 


good in most cases 


SUBSCRIBER DIALING HABITS AS OBSERVED WITH A MONITORING CIRCUIT 


ON A SENDER WITH INDUCED DIAL TONE DELAYS® 


\s a separate study a series of tests was made by means of a moni- 
toring circuit on one of the senders serving in common the subscribers 
in the Sterling-3 and Main 2 central offices, for the purpose of obtaining 
further information on subseriber dialing characteristics under overload 
conditions. A large amount of data was collected on the time intervals 
from the seizure of the sender to the first action taken by subscribers 
when encountering dial tone delays, the latter being introduced under 
the control of the observer. 

The monitoring circuit was wired to a particular sender in a group 
of 100 serving all classes of subscribers. When the circuit was in use, the 
only irregularity introduced was that the dial tone could be delayed 
even though the sender was actually available to the subscriber. The 
delay did not affect the sender in its functions if the subscriber elected 
to dial before tone 

The sender monitoring cireuit provided the following four features: 

1. A receiver was bridged across the tip and ring leads in the sender 
so that an observer could hear certain actions taken by a subseriber 
connected to the sender. The sender was of course disconnected before 
conversation. 

zZ The observer was able to preselect one of several intervals by 
which dial tone was delaved on successive calls served by the sender. 
‘This was accomplished with a capacitance resistance-vacuum tube cir- 
euit 

3. By means of a timer which started when the sender was seized, 
the observer was enabled to note elapsed time intervals to the oecurrence 
of the various actions of the subscribers. The reading of the time of 
the first action of a subseriber had to be made when the second hand 
Was in motion, which introduced certain errors later to be discussed. 

1. By means of colored lamps the observer was able to classify all 
ealls observed as being message rate, flat rate or coin. 


During the sender dial tone delay tests, observations were made only 


during the afternoons when the flow of traffic was light and the prob- 


ability of a subscriber obtaining a delay before reaching the sender was 


& minimum 
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DIALING HABITS OF TELEPHONE CUSTOMERS 


The observer was provided with a means for introducing either no 
delay or one of four values of delay 2, 5, 10 or 15 seconds into the sender 
dial tone circuit. The observer took 50 observations using a particular 
value of dial tone delay and then shifted to another so that no particu- 
lar value of delay would become evident to the customers during an 
afternoon’s test. Each group of 50 observations comprised a mixture of 
message rate, coin and flat rate calls in the approximate proportions of 
13 to 6 to 1, representing the respective volumes of traffic from these 
classes of service during the afternoon periods. It was not possible to 
distinguish PBX lines or two-party lines from the bulk of the message 
rate data nor PBX lines in the flat rate data, although to a limited 
extent the observer could identify PBX dialing by the generally faster 
pulsing. The coin data represent both public and semi-public customers. 

Fig. 12 is a diagram for explaining the results shown on Figs. 13, 14 
and 15 for the message rate, flat rate and coin classes of service, respec- 
tively as obtained with the sender monitoring circuit. Fig. 12 was ob- 
tained by the application of fitting curves to those message rate data of 
Fig. 13 for which a dial tone delay of five seconds was introduced by the 
observer. In the interval from ¢t = 0 to t = 5 seconds, three curves A, B 
and C represent the per cent of subscribers still waiting at time ¢ for dial 
tone. Curve A and its extension beyond ¢ 5 seconds represents the 
action of subscribers who would dial their calls before tone if dial tone 


were delayed indefinitely. Curve B and its extension beyond ¢ 5 
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Fig. 12—Explanatory chart for sender monitoring observations; dial tone at 
’ = 5 seconds. 
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seconds represents actions of subscribers who would disconnect if dial 
tone were delayed indefinitely. Curve C and its extension is the sum of 
the other two. Curve D in the region beyond t = 5 seconds represents the 
actions of subscribers who disconnect after tone, curve E represents the 
actions of subscribers who will dial their calls after tone and curve F 
represents the sum of the lower curves. Of interest is the fact that for an 
interval of about two seconds following dial tone (at ¢ 5 seconds), 


the observed total survivor curve F lies above the extended portion of 
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the hypothetical survivor curve C for infinitely delayed dial tone. This 
indicates that most of the subseribers who would have abandoned their 


attempts during this interval abruptly changed their minds and then 


consumed a noticeable interval of time after hearing tone before start- 
ing to dial. Thus, as might be anticipated, the subscribers exhibit a 
reaction time. 

Fig. 13 shows the results in terms of survivor curves that were ob- 
served for the message rate class of service. Five sets of curves are 
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”) 


presented, namely for 0, 2, 5, 10 and 15 seconds of dial tone delays from 
the instant of sender seizure. It should be noted that the general con- 
tour of the various curves up to the receipt of dial tone and when ex- 
tended beyond gives an estimate of the survivor curve for dial tone 
delayed indefinitely. Fig. 14 shows the results for the flat rate class of 
service (the data here are relatively meager), and Fig. 15 shows the 
results for coin customers 


Fig. 16 indicates for the three classes of service the progressive changes, 
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tesults of sender monitoring observations 


with increasing dial tone delay, of the per cent sender seizures resulting 
in dialing before tone, in dialing after tone and in disconnections for the 
five dial tone delay intervals studied. 

Some general comments concerning Figs. 13 to 16 may be made. 

1. There is a striking contrast between the message rate and coin 
classes of service. This may be due to the immediate financial stake 
that a coin customer has in his call. He is reluctant to disconnect before 
dial tone. 


2. The results for the message rate and flat rate classes of service 


appear to be similar at the shorter dial tone delays; at the longer delays 
a higher proportion of flat rate subscribers have already dialed before 


dial tone. This apparent discrepancy may be due to the relatively small 


number of flat rate observations. Flat rate service in Sterling-3 and Main- 
2 was principally for professional people, such as doctors and nurses, 
who were thought to be more demanding than ordinary subscribers for 
prompt dial tone service. These results therefore should not be con- 
sidered characteristic of flat rate customers generally. 

3. Detailed analysis of the data (results not presented in this article) 
indicated that the distributions of time to first pulse for subscribers not 
observing tone, and for those waiting for tone, are quite similar for the 
three classes of service. 

$1. Only the unsmoothed raw data have been shown on Figs. 13 to 
15 since certain inadequacies were detected in the observations. These 
were due to observer reactions and reading errors discovered as a result 
of comparing preliminary practice sender monitor test results with a 
simultaneous record obtained by the 100 pen tape recorder. The overall 
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is that the data obtained by the observer are generally displaced 
vard along the time axis by about 0.8 second. 


Phe message rate data were for individual, PBX and two-party 
subseribers and the flat rate data were for individual and PBX sub- 


eribers. Furthermore, certain of the flat rate subscribers had auxiliary 


message rate service. [It seems likely that different characteristics would 
be obtained for the individual, PBX and two-party subscribers since 
there appear to be reasons for expecting significant differences in their 
dialing habits. The PBX operator is in a position to “shop” for tele 
phone service. If she fails to get dial tone on one outgoing line, she can 
try anv other free line. This can also be done by subseribers with multi- 
line service. This “shopping” for service tends to produce a large volume 
of disconnections when dial tone is slow. The individual and the two- 
party stubse ribers cannot do this and hence they can be expected to 
show lewel disconnections. 

6. The results are in terms of intervals of time from the instant the 
sender is seized. It would, of course, be preferable to have these results 
in terms of time from receiver off hook. Since on the average the sender 
is seized in a time interval of about the same magnitude as that of the 
reaction time of the observer, Figs. 13 to 15 can be read approximately 
correctly when the abscissas are redesignated ‘time in seconds from 
receiver off hook.’ The foregoing results have been presented to furnish 
in jnereased understanding of subseriber dialing habits. In the next 
section additional results based on individual line records taken on the 


Tiuap ire presented, 


SUBSCRIBER DIALING HABITS OBSERVED BY INDIVIDUAL LINE RECORDS 


As indicated in the previous section, the results obtained by means 
of the sender monitor tests were subject to certain shortcomings, hence 
data taken on individual lines with the 100-pen recorder have been 
analyzed to augment the information concerning the dialing habits of 
subscribe rs 

\s noted heretofore, several of the pens on the L00-pen recorder were 
available for taking observations on subseribers lines. Two pens were 
used per subseriber line, one recorded the operation of the subseriber’s 

vhile the other pen marked whenever the subscriber’s line 
On an originating eall both pens started marking when the 
initiated a call. When a line finder was obtained, the line 
ceased marking and it was presumed that the subscriber 


lial tone at that instant. Dialing, hang-up and flashing by a 
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subscriber after receipt of dial tone are noted by breaks in the markings 
of the line busy pen. Dialing, hang-up or flashing before dial tone are 
noted by simultaneous breaks in the markings of both pens. Various 
intervals can be measured and the call attempts classified accordingly. 
Observations were obtained in the foregoing manner during the course 
of the Sterling-3 line finder tests on the following numbers of subscriber 
lines: 
Message rate—residential 
business 
PBX 
two party 
Flat rate individual 
Coin 


IS2 


The observed data were classified for each of the above six types of 
subscribers in terms of the following categories: 
1. Time to subscriber action before receipt of dial tone. 
a. Time from receiver off hook® to first digit dialed by subscriber. 
b. Time from receiver off hook to disconnecting action by the 
subscriber. 
2. Time from receiver off hook to receipt of dial tone. 
3. Time to. subscriber action after receipt of dial tone. 
a. Time from receipt of dial tone to first digit dialed by subscriber. 
b. Time from receipt of dial tone to disconnecting action by the 
subscriber. 
Because the data developed in this section are compared with both 


the 7 factor analysis and the sender monitor test results and because 


the treatment of subscriber dialing habits before dial tone for each ot 
these items is different, categories la and 1b are analyzed in two ways. 
In the 7 factor analysis, all actions of a subscriber prior to dial tone, 
except a disconnect, were considered to be one continuing demand for 
service. Thus for the first analysis, cases of dialing, flashing and short 
disconnections before tone lasting less than three seconds were ignored. 
In the sender monitor tests the only items considered were the time to 
the first digit dialed by a subscriber and the time, if no dialing occurred- 
to the release of the sender by the subscriber. Thus for the second analy- 
cases of dialing before tone and flashing or disconnections before tone 


SIS, 


® When a customer initiates a call, the line relay operates. For individual line 
and two-party subscribers this occurs when the subscriber takes his receiver off 
the hook. For coin customers this occurs when the customer has taken his re 
ceiver off the hook and made a proper deposit. For a PBX line this occurs when 
the PBX attendant has established a connection to an outside line. All of this 
is collectively termed ‘“‘receiver off hook.”’ 
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Fig. 17 —Subseriber dialing habits based on individual line reeords when dia 
tone is delayed indefinitely 


that caused the sender to release were each counted as separate attempts. 


The results arrived at are shown as survivor curves on Figs. 17(a) to 
19(b) for each of the six types of subseribers studied. 

he survivor curves were developed by considering events during 
0.1 second intervals. The number of cases of dialing before tone and 
the number of cases of disconnection before tone occurring during a 
0.1 second interval were divided by the number of cases waiting for dial 
tone at the start of the interval. This ratio was considered to be a 
retirement rate. The complement of this rate gave a survival rate. By 
a progressive multiplication of survival rates from time, ¢ 0, the 


resulting survivor curves were obtained. Those cases receiving dial tone 
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during a particular 0.1 second interval are omitted from the number of 
cases waiting for dial tone at the start of the next interval. 

In the development of the generalized trunking formula, the assump 
tion was made that the waiting times of calls infinitely delayed have an 
exponential distribution. By assuming that the plots for the survivor 
curves in the development comparable to the j factor analysis are ex 
ponential distributions, it is possible by reading the value of ¢ correspond- 
ing to 36.8 per cent of the subscribers still waiting for dial tone to obtain 
estimates of the values of H for the six types of subseribers. These esti- 
mates, most of which were obtained by extrapolation, are compared in 
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Fig. 18—Subscriber dialing habits based on individual line records when dial 
tone is delayed indefinitely. 
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Subseriber dialing habits based on individual line records when dial 


ed indefinitel 
IV with the values developed earlier in connection with the 

} factor analysis. 
Phese results are not considered to be inconsistent since the tails of 
the survivor curves were constructed from very meager data. Conclu- 
sions based on Figs. 17(a) to 19(b) should therefore be regarded as 


having qualitative value only. The results principally indicate that the 
wiiting-for-dial-tone characteristics of subseribers clearly vary with the 


different classes of service. 
(‘omparisons between the lower survivor curves on Figs. 17(a) to 
19(b) and the curves on Figs. 13 to 15 of the sender monitor tests are 


indicated by the percentages given in Table V of subscribers waiting for 


dial tone 5, LO and 15 seconds from the time they requested service. 
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These results agree reasonably well when it is recalled that parts of 
the individual line data were scanty and that the sender monitor tests 
inclided the effeets of observer reactions. 

Onee dial tone is received, it appears that all types of subscribers 

nd to follow a uniform dialing pattern. Figs. 20 to 23 show for a class 
interval of 0.5 second the distributions of the per cent of subscribers who 
dial at time ¢ for the six types of subscribers studied. Figs. 20, 21 and 22 
show the distributions when dial tone is received from 0.0 to 0.5, 0.6 to 


0.40 and 1.0 to 1.9 seconds after dial tone. respectively. These curves 
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Fig. 23-——Distributions of the start-to-dial times of subseribers who dial 
the receipt of dial tone 


indicate that a strong similarity exists among the six types of subscribers 
with regard to their dialing patterns once dial tone is received. 

Fig. 23 shows composite dialing distributions of the six types of sub 
sceribers for four dial tone delay intervals. For dial tone delays less than 
3 seconds the dialing patterns are seen to be similar. Beyond 3 seconds 
delay the start-to-dial curve shifts outward, although the data are too 


scattered to indicate closely where the movement begins. 


CONCLUSION 


The foregoing report of the results of the tests conducted at the 
Sterling-3 central office indicates that the dialing habits of subscribers 
Waiting for dial tone can be observed and analyzed to develop so-called 
) factors for use in a more general trunking formula than has been 
employed until recently in the Bell System. The report also presents 
descriptive data regarding the patterns subscribers follow when waiting 
for dial tone. 


ACKNOWLEDGEMENT 


The writers gratefully acknowledge the help given by the personnel 
of the Long Island Area of the New York Telephone Company, by 


their associates and others in the various phases of the study and pal 


ticularly by C. F. Bischoff who had a major part in conducting the tests 





Selective Fading of Microwaves 


BY A.B. CRAWFORD AND W)C JAKES, JR 


Manuscript ived Octo 


resulls of an extended survey of microwave propagation over two line- 
ght paths in AN ew Jerse / are dese thed, Angli of-arrwal measurements 
at 1.25-cm wavelength and selective fading observations in a 450-me fre . 
quency band centered at 8950-me show that the seve re fading can be explained 
in terms of multiple-path transmission, A computer of the analogue type 


was built to simulate the more complicated selective fading patterns. 


INTRODUCTION 


During the past few years, studies of microwave propagation have 
been made by the Radio Research group at the Holmdel Laboratory 
over two paths located in eastern New Jersey. Both of these are line- 
ol sight paths which might be considered to be typical links Ih & Cross- 
country microwave radio relay emreult. 

In conducting these studies, the usual continuous recordings of signal 
levels were made but the greater interest was centered in special experi- 
ments designed to reveal more ol the processes which can Cause fading. 
Che most relevant information has been obtained by exploring the inci- 
dent wavefronts with a narrow-beam scanning antenna (angle-of-arrival 
studies) and, more recently, by observing the transmission characteris- 
tics of the paths by means of a frequency-sweep technique and also by 
the use of very short pulses 

Some results of angle of-arrival observations have been reported pre- 
viously’ and a companion paper describes the transmission tests con- 
ducted with very short pulses.? The present paper describes some of the 
observed mechanisms associated with fading, presents typical data ob- 
tained with the narrow-beam scanning antenna and gives examples of 


the frequency-sweep observations, illustrating the frequency selective 


W. M. Sharpless Measurement of the Angle of Arrival of Microwaves,”’ 
Proc. I.RLE., 34, Nov. 1946, pp. 887-845. A. B. Crawford and W. M. Sharpless, 
Further Observations of the Angle ot Arrival ot Microwaves,”’ Proce I.R E ; 
34, Nov. 1946, pp. 845-848. H. T. Friis, ‘Microwave Repeater Research,’’ Bell 
Sustem Tech. J., 27, Part I, “Propagation Studies’ by A. B. Crawford, Apr 
1948, pp. 183-246 
-O bk. DeLange, “Propagation Studies at Microwave Frequencies by Means 


Very SI Pulses,’’ Bell System Tech. J., 31, Jan 1952, pp 91-1938 


GOS 
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nature of the fading. Some data derived from the continuous recordings 
of signal levels are presented in an appendix. 

The angle-of-arrival observations were made at a frequency of 24,000 
megacycles. The frequency-sweep experiment and the recordings of signal 
levels were made in the 3700 to 4200 megacycle frequency band as were 


the short pulse observations described in the companion paper 


GENERAL DISCUSSION OF PROPAGATION PHENOMENA 


The map of Fig. | shows the location of the experimental transmission 
paths. The path between Crawford Hill and Southard Hill is 17 miles 
long and clears the intervening terrain by 65 feet, approximately one 
Fresnel zone at a frequency of 4000 megacycles. The other path, between 
Crawford Hill and a 100-foot tower on the Murray Hill Laboratory 
property, is 22.8 miles long and has clearance of 280 feet. Fig. 2 shows 
the profiles of these two paths. 

The general characteristics of over-land microwave transmission are 
well known and need be reviewed only briefly. During the daytime 
hours, when the lower atmosphere is thoroughly mixed by rising cori 
vection currents and winds, the signals are normally stable and are near 
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Fig. 1—Map showing location of the experimental transmission path 
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the free space levels. Also during the winter months, when the humidity 

ontent of the atmosphere is low, signal variations are usually very small. 

on clear summer nights with little or no wind, non-uniform 

lons Of temperature and humidity can create steep dielectric 

nt gradients in the lower atmosphere, thus causing anomalous 
propagation and fading 

When fading occurred on out experimental transmission paths, an 

tlarm circuit connected into the continuously recording equipment was 

irranged to operate when the signal level dropped below a predetermined 


alue. This enabled observers to be present during severe fading periods 


SOUTHARD 








Profiles of the transmission paths 
ans of the special experiments, to determine the causes 


though it has not been possible to provide satisfactory explanations 


forall of the observed fading phenomena, much of the fading (occasions 


Vhen the signals are depressed to levels Lo to ) decibels or more below 


the normal daytime value) can be explained qualitatively in terms of 


imple rav pietures. Fig. 3 is intended to illustrate some of the observed 


fading phenomena. The case of multiple path transmission, the most 
common cause of fading on either transmission path, is shown in Fig. 3(a) 
Pwo, three and sometimes more signal components are found to arrive 
us angles in the vertical plane, usually above the line of sight 
interference among these components produces fading, the severity 

ch depends upon the relative amplitudes and delays of the com 

\t these times, different frequencies fade differently and the 

: received on two vertically spaced antennas also fade differently. 

of either frequency or space diversity would be effective in this 

latively rare type of fading, observed only on the Murray Hill 


believed to be caused by the mechanism illustrated in Fig. 3(b 
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Here a reflecting layer is situated between the heights of the transmitter 


and receiver. The signal then suffers attenuation due to reflection of 


part of the energy from the direct path. Widely separated trequencies 


are affected in like fashion and the outputs of antennas spaced fon 
diversity reception tend to be in agreement although the fine structure 
fading is usually different. 

On neither of the experimental transmission paths is there a rgeulat 
ground-reflected component of any consequence. Due to the roughness 
of the ground and the presence of vegetation, the effective reflection coet 
ficient is of the order of 0.2 for either path. Ground reflections thus play 














(d) 

Fig 3 Possible ray paths involved in severe fading i Multiple path trans 
mission. (b) Attenuation by reflection from an elevated laver. (¢) Abnormal 
water reflection on the Murray Hill-Crawford Hill path. (d) Substandard con 
ditions on the Southard Hill-Crawfard Hill path 
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no significant part in the fading picture with the xceptlon ol the situa 
tion illustrated in Fig. 3(c¢). Occasionally on the Murray Hill path, 
conditions of atmospheric retraction are such that a strong signal com 
ponent is received by virtue of reflection from the water surface of 
Raritan Bay. Under normal conditions, the geometry of the path does 
not permit suc ha reflection 

Normally the dielectric constant of the atmosphere decreases with 
height above ground so that the ray path usually has a curvature in 
the same direction as the earth curvature However, it is possible for the 
dielectric constant of the atmosphere to increase with height above 
ground (sub-standard conditions) so that the ray path has a curvature 
Opposite that of the earth. This results in the condition illustrated in 
hig. 3(d) where the limiting or tangent ray does not reach the receiver 
and only a weak signal is received by virtue of diffraction. Widely sepa 
rated frequencies and vertically spaced antennas are affected alike as 


regards the average signal level but not the fine structure fading. This 


effect has been observed only on the Southard Hill-Crawford Hill path 


Which has small clearance to begin with. It has been observed on several 
nights in late summer or early autumn after a radiation type ground 
fog has formed in the late evening and usually persists until the fog is 
dispelled by winds or by the morning sun 

Phere are, of course, times when the transmission conditions are con- 
siderably more complicated than those described above. Some of these 
apparently are due to a combination of the situations illustrated in big. 3 
while others may be the result of an atmospheric focussing or trapping 
phenomenon. In addition to the various phenomena just described, which, 
fortunately, occur rather infrequently, there are numberous occasions 
when the signal varies plus and minus a few decibels relative to the free 
space level. It has not been possible actually to demonstrate the mecha 
nism responsible but it seems most likely that these smaller variations 
are due to non-linear dielectric constant gradients which give the atmos 
phere the properties of a convergent or divergent lens. 

\n important result of the observations made to date is the con- 
vietion that the severe fades, signal exeursions to levels 30 decibels or 
more below the tree space field, were all caused by wave interference. It 
appears that, as the average signal level is depressed by any mechanism, 
it becomes more and more vulnerable to the effects of extra signal 
components of small amplitude that often may be present but go un- 
noticed when the signal is near normal levels. Thus, while the average 
signal level during the conditions illustrated in Figs. 3(b) and 3(d) may 


be no more than 15 to 20 decibels below the normal daytime level, there 
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is usually superimposed a fine structure fading in which short duration 
fades to levels as much as 45 decibels below free space have been ob- 
served. For this reason it is desirable to avoid paths having small clear- 
ance over intervening terrain and also paths which have a permanent 
ground reflection of sufficient magnitude to depress the signal to critical 
levels when, due to variable atmospheric refraction, the direct and re 
flected components are in phase opposition. 

The following sections describe the angle-of-arrival and frequency- 
sweep experiments on which much of the preceding discussion was 
based. 


ANGLE-OF-ARRIVAL OBSERVATIONS 


A photograph of the Crawford Hill receiving site is shown in Fig. 4. 
The building housing the receiving equipment and the associated an- 
tennas are mounted on a framework which can be rotated on a concrete 
track, permitting investigation of the transmission characteristics of 
either path. The parabolic antennas on the tower are used for continuous 


recording of 4195 megacycle signals. The long object at the left of the 
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picture is the metal-lens antenna used for making the angle-of-arrival 
observations. Its half-power beamwidth is 0.12 degree at the operating 
frequency of 24,000 me. The focal length of the lens is 48 feet and its 
feed is located in the little cupola on top of the building. The feed is 
held fixed, while the lens is moved vertically by a motor-driven mecha- 
nism; thus the antenna beam also moves vertically. The antenna scans a 
total angle of two degrees in ten seconds. It is fed by a 24,000-me radar 
set which is gated to receive only the pulses reflected from a corner 


reflector located at the distant terminal of the transmission path. The 
pot on the radar cathode ray tube moves vertically in’ synchronism 

ith the scanning antenna, and the horizontal deflection is proportional 
to the amplitude of the pulse received from the corner reflector. The 
display thus shows amplitude of the various incoming signal components 


as a tunetion ol their angles of arrival 


The antenna installation on Southard Hill is shown in hig. 5. At the 


left is the transmitting paraboloid for the 4195-me continuous wave 
transmitter, the radar corner reflector is in the center, and on the right 
is the horn-reflector antenna used in the lrequency-sweep experiments 
described below. Similar equipment is located at the Murray Hill termi- 
nus. The corner reflector is 5.5 feet on a side, and at 24,000-me has 
sufficient gain to override reflections from other nearby objects, and 

thus heeomes easily identifiable on the radar screen 
Phe radar oscilloscope for typical propagation conditions is shown in 
Fig. 6. These pictures were obtained by leaving the camera shutter open 
during the ten-second interval required for the antenna beam to scan 
through the angular range of 2°. All of these representative photographs 
vere taken on the Murray Hill-Crawford Hill path although similar 
results were obtained on the Southard Hill-Crawford Hill path with the 
exception of Fig. 64). The normal daytime transmission is shown in Fig. 
H(a) to consist of a single path arriving at an angle of —0.2° with respect 
to n fixed reference angle. The horizontal lines represent intervals of 
so that changes of 0.05° can be estimated. The other pictures in 

6 were all taken during fading conditions. 

igs. O(b) and 6(¢) are good examples of the multiple-path condition 
shown in Fig. 8(a) in which the individual components are almost equal 
in amplitude and well separated in angle. In Fig. 6(b) there are two 
components arriving at angles of O.1° and 0.6° above the normal line-of- 
sight while in 6(e) there are three components with angles of 0.05°, 0.35 
thove the normal angle. The position and amplitude of the 
omponents may change radically in a matter of minutes, and 


re is ho component that can be identified as the “normal” one. 
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Fig. 6—Representative photographs of angle-of-arrival observations on the 
Murray Hill-Crawford Hill path. (a) Normal day. (b) Two elevated paths. Sept 
S, 1950; 9:23 p.m. (c) Three elevated paths. Aug. 27, 1950; 1:11 a.m. (d) Multiple 
paths. August 26, 1950; 11:00 p.m. (e Wide angle ‘‘fill-in’’. Aug. 26, 1950; 11:04 
p.m. (f) Abnormal water reflection, Sept. 8, 1950; 11:28 p.m 


) 


During these multiple path conditions, the recordings of the 4195-me 
transmission generally show the broad maxima and sharp minima charae- 
teristic of wave interference. 


Figure 6(d) shows a case in which the various paths are not completely 


separated while Fig. 6(e) (taken four minutes later) shows that energy 


is being received almost without variation over a vertical angle of 0.4°. 
This may represent a number of ray paths which would be separable by 
a narrower-beam antenna, or it may indicate a focussing or trapping 
phenomenon. Often when the type of transmission illustrated by 6(e) 
is present, the recorded 4195-me signal may be as much as 12 to 15 
decibels above the free space levels, 
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Fig. 6(f) illustrates the case of abnormal reflection from the water of 
Raritan Bay on the Murray Hill path as indicated in Fig. 3(c). Here the 
“normal” signal component is arriving at 0.1° above the line-ot-sight 
while another component, almost equal in amplitude, is arriving at the 
very bottom of the scan, about 0.8° below the line-of-sight. It is quite 
probable that there have been times when this component was present 
but was outside the range of the scanning antenna. 

The mechanisms discussed in connection with Fig. 3(b) and 3(d) can- 
not be demonstrated by photographs such as those just presented al- 
though the angle-of-arrival radar was instrumental in furnishing the 
clues to the phenomena. Due to the two-way attenuation of the radar- 
corner reflector technique, the signal at these times rapidly falls below 
the noise level of the receiver. For the same reason, it is not possible to 
detect the extra signal components of small amplitude which were postu 
lated to account for the very deep fades sometimes observed under these 
transmission conditions. 


FREQUENCY-SWEEP OBSERVATIONS 


Since most of the fading is due to interference between waves which 
travel over different paths of, presumably, different lengths it was 
realized that the fading was likely to be frequency selective. Just how 
selective would depend on the relative lengths of the individual trans 
mission paths. The usual methods for determining path length differ- 
ences are to use short pulses, or to sweep the frequency. Since it was 
likely that the path-length differences would be measured in feet rather 
than yards, very short pulses or a wide frequency-sweep were required 
An oscillator? was available whose frequency could be swept over the 
licensed band of 500 me between 3700 me and 4200 me. The frequency 


sweep experiment was set up on the Murray Hill-Crawford Hill path 


for the summer of 1949. The following summer, the milli-microsecond 
pulse transmission tests described in the companion paper were con 
ducted over the same path. As might be expected, simultaneous observa- 
tions showed good agreement between the two methods. 

The frequency of the transmitter, located at Murray Hill, is swept 
over a 450-me band centered at 3950 me at a 60-cycle rate. At the 
receiver, a similar oscillator is used for the beating oscillator except that 
its frequency is swept linearly through the same frequeney band in one 

This oscillator was developed by M. EE. Hines and is deseribed in his paper 
published in the Bell System Technical Journal, Vol. 29, Oct. 1950. It uses a 416A 
close-spaced triode in a wave-guide cavity. The frequency is changed by means 


of a plunger which is capacitv-coupled to the plate of the tube and which is aetu- 
ated by a modified loud speaker unit 
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econd, Since the intermediate frequency amplifier of the receiver is only 

350 ke wide, (centered at 600 ke) narrow pulses are generated each time 
the frequeney of the transmitter crosses the frequency to Which the 
receive! tuned. These intermediate frequency pulses are displayed 
ertically on a eathode ray tube. The horizontal trace is synchronized 
vith the one-second sweep rate ol the beating oscillator. 


The normal daytime frequency-sweep pattern is shown in Fig. 7(a 


Phe vertical seale is linear in amplitude and the horizontal scale is almost 


linear in frequeney, with frequency decreasing from left to right. Visible 
at the extreme left is the signal used for continuous recording. Since 
there is only one transmission path involved, the amplitude of the re 

ved signal is nearly constant over the 450-me band. If another signal 
were present Which had travelled over a path of different length, the two 


difference is an even multiple of half-wavelengths and subtract when 


nals would add when the frequency is such that the path length 


the path length difference is an odd multiple of half-wavelengths. Simple 
ealeulation shows that if the path length difference is one foot, the 
frequencies at which the signals add and subtract are separated about 
500 me. Thus the limit of resolution for the frequency-sweep experiment 
isa little more than one foot 

Photographs taken on a night when the angle-of-arrival radar indicated 
two almost equal components separated about 0.4 degrees in angle are 
hown in Fig. 7(b). The time interval between the two pictures is 30 
seconds, during which the minimum had shifted about 150 me. The 
pietures can be interpreted as simple two path transmission with an 
ndieated path difference of about two feet and an amplitude ratio of 
O.7 to 1. On this night the minimum shifted back and forth across the 
frequency band—sometimes slowly and sometimes rapidly. At times the 
position of the minimum might remain fixed but its depth would change 

Photographs taken on a might when there were abnormal reflections 
from the water of Raritan Bay are shown in Fig. 7(¢). There are evidently 
two main components with path difference of about six feet, with a small 
third component causing the slight decrease in amplitude of the peaks 
from left to right. These pictures were taken 9 minutes apart, but this 


type if i nh Was obser. ed over fn period ol about three hours On this 


he tre queneyv sweep patterns are considerably more compli- 
the ones shown so far. Fig. 7(d) shows two photographs which 
nal components and perhaps more were 


ist three =] 


u 


al between the two pictures was about 30 
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Fig. 7—Representative frequeney-sweep patterns observed on the Murray 
Hill-Crawford Hill path. (Summer 1949 t) Normal day. (b) Two components 
with a path difference of two feet. (¢) Two main components with a path differ 
ence of about six feet plus a small third component. (d) Multiple component 
pattern. 
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OF FREQUENCY-SWEEP) PATTERNS 


Po aidcin the Interpretation of the complicated frequency sweep pat- 
i computor of the analogue type was built. This apparatus com- 
bines four signal components, three of which are variable in delay and 

and presents the result on a cathode ray tube in the same 


terns 


amplitude, 
form as the actual frequency sweep patterns. Thus a particular pattern 
can be synthesized on the computor and the number of components, 
together with their path differences and relative amplitudes, read directly 
from the computor dials. This is accomplished by generating four 600-ke 
signals, three of which are phase modulated at 60 cycles per second. The 
total phase deviation and relative signal amplitude are variable. The 
four signals are then summed and displayed in vertical deflection on a 
cathode ray tube hav ing a OO cycle horizontal sweep. 

The synthesis of the patterns of Figures 7(b) and 7(c) are shown in 
Fig. 8. The upper synthesized pattern is simply a combination of two 
components with relative amplitudes of 0.7 and 1 and a path difference 
of two feet. The lower pattern consists of the reference component with 


unity amplitude, a second component with an amplitude of 0.5 and a 


path difference of 5.7 feet, and a third component with an amplitude of 


0.2 and path difference of O.S8 feet. The similarity between the actual 


and synthesized patterns is obvious. 


Fig. S Synthesis of the frequency sweep patterns of Figs. 7(b) and 7(e 





Fig. 9—Synthesis of complicated frequency-sweep pattefns using four com- 
ponents. See Table I for values of the relative amplitudes and path differences 


Sl 
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examples of atte mpts to sy nthesize some more complicated frequency 

ns. taken on the mght of August 2, 1949, are shown in 

bout Components were required in each CuSse, with the path 

ind delays being summarized in Table I. Although the pie- 

all appear very different, in general major changes were required 

n component No. 2 to go from one pattern from the next, as Table I 

\ll the remaining components had to be very carefully trimmed 

in both amplitude and delay to get good synthesis (especially in the case 
f big. Od), but these changes were relatively small. 


Ol 


CONCLUDING REMARKS 


Lhe ~pecial experiments just described have led to the conclusion, 
ssed earlier, that the severe fading observed on the two test paths 

sult of multiple-path transmission in which several components 

volved. These components may arrive at the receiver at various 

iigies up to three quarters Of a degree above the normal daytime angle 


of-arrival and, in the case of abnormal water reflection on the Murray 


Hill path, as much as OLS degree below the normal angle. The path 


differences among these components may vary from a fraction of a foot 
to about ten feet. The long-delay components are usually small in ampli 
tude 
fn all cases where observations were made during periods ot excep 
tionally high signal levels, sav 10 to 15 decibels above free space level, 
thie rrequeney-sweep patterns were substantially flit. suggesting a foeus- 
or trapping phenomenon. The frequencey-sweep patterns were also 
m those nights when the signal excursions were only a few decibels 


ind below the normal davtime level. However, the severe fades 


TABLE | 


1s AM 12:24 AM 
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Fig. 10--A sequence of angle-of-arrival and frequeney-sweep patterns taker 
at ten-second intervals 
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of short duration sometimes observed when the average signal level was 
ce pressed by the mechanisms of Figs. 3 b) or 3(d) were found to he fre- 
queney se lective. 


Some of the studies described in this paper were made with vertically 


polarized waves and some with horizontally polarized waves; at times, 
15° polarization was used. In so far as it was possible to determine, the 
propagation characteristics of both paths were independent of the polari- 
zation used. 

No meteorological soundings were made in connection with this work. 
Considering the rapid changes usually observed with the angle-of-arrival 
and frequency sweep apparatus, it is doubtful that meteorological meas- 
urements made in the usual manner would show much correlation with 
the radio observations except, perhaps, in a general way. The sequence 
of pictures in Fig. 10 is included to show how the angle-of-arrival and 




















tribution curves of the signal levels observed on the Murray 
th. Data of 1947, 1948, 1949 and 1950 
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frequency-sweep patterns change with time. These pictures were taken 
at LO-second intervals. On this oecasion there was good correlation be 
tween the angle-of-arrival and frequency-sweep data. Such was not 


always the case, however, and considering the wide difference in operat 


ing frequencies, 24,000 me and 4000 me, instantaneous correlation should 


hot necessarily be expected. 

Although all the studies described in this paper were made on the two 
local paths, the results are compatible with propagation measurements 
made by another group in the Laboratories during a survey for the 


transcontinental radio relay system. 
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APPENDIN 


This appendix is included to illustrate some of the characteristics of 
the propagation as shown by the recordings of 4195-me signal levels. 
hig. Jibs a reproduction of some typical records obtained during severe 
fading periods. Fig. 1 PGi) is an example of transmission during the sub 
tandard conditions illustrated by the ray diagram of Fig. 3(d). Fig. 
l1(b) is typical of multiple-path type fading in which the signal com- 
ponents arrive from elevated angles as shown in Fig. 3(a), while Fig. 11 (¢ 
was recorded on a night when, fora time, there were abnormal reflections 
from the water of Raritan Bay on the Murray Hill path, see Figs. 3(¢) 
H(t) and Tle). The records of Fig. 11(d) show how the outputs of two 
similar antennas, spaced vertically about 30 feet, differ in regard to the 
deep fades of short duration. 

Phe chart of Fig. 12 shows how the fading varies with the time of year. 
this chart, the vertical lines represent the extremes in signal level 
observed during the twenty lout hour period from noon to noon. The 
large signal variations are concentrated mainly in the summer months. 
Phe time distribution of the signal levels recorded on the Murray 
Hill-Crawford Hill path are shown in Fig. 13. Each of the curves is for 
four-month period: the period of least fading, December, January, 
February and March; the period of most fading, June, July, August 
uid September; and the in-between period consisting of April, May, 
October and November. Data obtained in the vears 1947, 1948, 1949 and 
1950 are included. Fig. 14 shows similar data for the Southard Hill- 
Crawtord Hill path. Phe hump in the time distribution curve for the 
months of \pril, Mav, October and November is due to substandard 
conditions, illustrated by the ray diagram of Fig. 3(d) and the typical 
of Pig. Tita), which affected transmission on this path during 


several mights in October, particularly in the vears 1947 and 1950. When 


it occurred, this type of transmission usually persisted for a period of 


il 


several hours 





Propagation Studies at Microwave 


Frequencies by Means of Very 
Short Pulses 


BY O. E. DE LANGE 


Manuscript received March 27, 1951 


Microwave pulses with a duration of about 0.003 microseconds were trans 
mitted over a 22-mile path from Murray Hill, N. J.. to Holmdel, N.J.. 
in order to determine the effects of the transmission medium upon such 
pulses. During “fading” periods multi-path transmission effects with path 
diffe rencees as qreal as 7 feet were observed, as well as some other effects 


A microwave frequency of 4000 megacycles was employed. 


INTRODUCTION 


This experiment was set up with two main purposes in view: First, 
asa means of studying microwave propagation, especially with regard to 
multi-path transmission effects and second, to determine the effect. of 
a transmission path upon the shapes of very short pulses, particularly 
to learn what restrictions might be imposed upon minimum pulse length 
or spacing between pulses by distortions produced in the transmission 
medium. 

In regard to multi-path transmission the pulse method seems to be 
the most straightforward way of studying such effects. For example, if 
there is transmission by more than one path, and if the pulses are suffi 
ciently short in comparison to the path length differences involved, then 
there will be received a separate pulse for each path. Under these condi 
tions the number of paths involved, path length differences and other 
information become directly evident. If pulse duration is too great with 
respect to the path differences involved, the pulses received via the various 
paths will overlap in time and the resultant multi-path effect will be 
pulse distortion rather than reception of individual pulses. This situa 
tion is much more difficult to analyze. 


TRANSMISSION PATH 
The transmission path is the same as that used by A. B. Crawford for 
microwave propagation studies by means of the frequency sweep method, 


9] 
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‘path from Murray Hhll, N. J., to ¢ rawford’s Hill (near Holmdel), 

Phe path length is approximately 22 miles, and is partly over 
vater and partly over rough land terrain. The frequency sweep studies 
had indicated that the path differences involved in multi-path fading 
were of the order of one or two to about seven feet. In terms of delay 
times this means differences of about 1 to 7 muillimicroseconds. In order 
to resolve the paths when the path differences were only one or two 
leet We should have liked to have pulses of about ] millimicrosecond 


duration. Because of the difficulties involved in generating, amplifying 


4000 MC 
PULSES 
we same 2 


ANTENNA \ |/ 


| } 
“a y 


RAVELING 
WAVE 
AMPLIFIER 








BASE-BAND | 
t SF 


GENERATOR 


MONITOR 


ind detecting such short pulses, we accepted pulses which, when dis 
plaved on our final indicating equipment, had a length of 3 millimicro 
seconds at half amplitude. \bout 6 millimicroseconds at the base.) In 


free space this pulse would be Just about 6 feet long at the base 


PRANSMEPTING EQUIPMENT 


Phe transmitter was mounted on top of a 100-foot tower at Murray 


ball \. can be seen from Fig. 1, it consisted of a c-Ww reflex oscillator 


opel iting at 4000 megacycles, a baseband pulse generator, a modulator, 


t 


or gate, Tor modulating these pulses on the microwave carrier, a single 
stage traveling-wave amplifier and finally a horn antenna. Approximately 
it of power Was obtained from the transmitter at the peaks of 

s. The antenna area was 25 square feet and its gain 32 db above 

i dipole. A pulse repetition frequency of 10 me was employed. 


C. Jakes, Jr. ‘Selective Fading of Microwaves,”’ 
1952. pp. 68-90 
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RECEIVING EQUIPMENT 


The receiving antenna, a large horn, was mounted between two poles 


guved for support. It had an aperture of about 90 square feet and a gain 


of approximately 38 db over a dipole. The receiver circuit is shown in Fig 


2. About 60 db of gain at 4000 me was provided by either two or three 
stages of traveling-wave tube amplifier depending upon the gain of the 
particular tubes used. It was necessary to provide very good shielding 
and also careful filtering of all power leads to eliminate the tendency fo 
this amplifier to sing. The amplifier fed two crystal detectors through a 
hybrid tee junction, Each detector emploved a silicon crystal of the 
IN23B type 

Two indicator circuits are shown in Fig. 2. These circuits are very 
similar except that one employed a vertical amplifier coupled to a Dumont 
5SXP2 CRO tube, whereas in the second the baseband output of the 
crystal was fed directly onto the deflection system of a traveling-wave 
type of CRO tube. The latter CRO tube, which has been deseribed by 
J. R. Pierce in the November, 1949, issue of Electronics, has a very high 
deflection SENSITIVITY and is used with a microscope to enlarge its trace; 
hence, no amplification was required between it and its driving crystal 
The deflection system of this tube has a bandwidth of 500 to 1000 me. 
The micro-oscilloscope Wis provided primarily lor photographing pulses 
by means of a 35-mm camera attached to the microscope. (Exposure 
time was 5 to 15 seconds. The time recorded for each picture corresponds 
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Fig. 2—Receiving equipment 
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hat at the end of exposure.) A second microscope made it possible 
to view and to photograph the screen of the tube simultaneously. The 
general procedure was to observe continuously during periods of dis- 
turbed transmission, taking pictures at regular intervals of 5 to 10 min- 
utes. When conditions were seen to be changing rapidly, pictures were 
taken much more frequently. The large oscilloscope with its vertical 
amplifier had a bandwidth of about 150 me and hence caused some 
deterioration of the pulse. It, however, was less tiring than the small 
scope, especially for long periods of observation and was watched 
follow the general trend of events. It was capable of resolving the pulses 
resulting from two-path transmission when the path differences were 
large 
The sweep circuits Tor the two indicator oscilloscopes were practically 
identical. The horizontal sweep voltage for each consisted of the linear 
portion of a sine wave which was generated by a c-w oscillator operating 
at one third of the pulse repetition frequency of 10 me. Each oscillator 
vas synchronized with the incoming pulses by means of a 10-me voltage 
derived by amplifving the pulse energy through a narrow band amplifier. 
This cireuit provided very satisfactory synchronization even during the 
times when signal amplitude Was SO low us to produce a very poor signal- 
to-noise ratio. ‘Timing markers were provided on each roll of film = by 
periodically photographing a series of pulses spaced by an interval of 9 


millimicroseconds. 


RESULTS OF THE EXPERIMENT 


Phe picture at the left of Fig. 3 shows the transmitted pulse. The 
right-hand picture shows the received pulse under what were considered 
to be normal transmission conditions. It is seen that, except for the addi- 
tion of noise and widening of the pulse due to passage through the ampli- 
fiers and other equipment, the pulse shape is unaffected. The time cali- 
bration on this and the following photographs are in millimicroseconds, 
each mark representing one millimicrosecond (Q.001 ys). 

During the summer of 1950, when this experiment was in progress, 


— = 
ee RTT 


transmitted pulse (right) received pulse—normal transmission. 
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. 4—Received pulses—disturbed transmission 


there was comparatively little fading over the path in question at micro 
wave frequencies. There were, however, a few nights of considerable 
activity and some interesting results were obtained. 


The series of pictures on Fig. 4 show one good example of multi-path 


transmission where the path length difference was great enough to pro 
duce complete resolution of pulses. At 11:05-20 there are two pulses, 
each 7 millimicroseconds wide at the base and with their peaks just 7 
millimicroseconds apart; in other words, the path difference was just 
sufficient to produce two pulses with no overlap. The pulse at the left 
is presumably coming by the main path and that at the right from some 
second path resulting from bending of the rays caused by atmospheric 
effects. At 11:07-20 the second path appears to have shortened, resulting 
in a path difference of only about 5 millimicroseconds. This may actually 
have been due to a change of length of the second path or it may have 
been due to distortion of the second pulse by energy coming by way of a 
third path. The pictures taken at 11:08-20 and 11:10-50 show evidence 
of transmission by a third path. In the first of these, for example, the 
width of the disturbance at the base line indicates the presence of the two 
original pulses spaced 7 millimicroseconds apart but the midpoint of the 
two no longer falls to the base line as was the case in the first picture. 
This could be accounted for by the presence of a third pulse coming over 
a path whose length was somewhere between that of the other two. Con- 
ditions obtaining at 11: 10-50 could also be accounted for by the presence 
of pulses from three paths, that is, energy coming by way of a third path 
might cancel part of one pulse and at the same time add to the other. 
This could account for the fact that the spurious pulse is larger than the 
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Fig. 5 Received pulses—disturbed transmission 


normal one. It is also possible that more than three paths were involved. 
Qn a number of other occasions the pulse coming by way of the second 
path appeared to be of greater amplitude than the one coming by the 
Thin path This same effect has been observed by Mr. Crawford and his 
colleagues on the angle of arrival equipment. 

Information obtained from the above set of pictures shows that for 
a time-division multiplex system using the length of pulse used here 
7 millimicroseconds at the base) and operating over this path, pulses 
would have to be spaced a minimum of about 14 millimicroseconds apart 
if it were desired to avoid distortion at all times. If very much shorter 
pulses were used the spacing might be reduced to 9 or 10 millimicro- 
seconds. However, the 7-foot path difference indicated by these pictures 
is about the maximum ever observed and occurs rather infrequently so 
that if somewhat closer spacings were employed troubles would result 
only a small percentage of the time. 

he next series of pictures, Fig. 5, taken July 8, show an example 
of a more common type of multi-path transmission. Here the path dif 
ference is apparently less than for the last series. At 11:15 there are two 
distinct pulses with an apparent path difference of about six feet (6 milli- 
microseconds) if judged from the spacing between the peaks of the pulses. 
However, from the length of the disturbance at the base line, which we 
consider a better criterion, the path difference was more nearly four feet 
\t 11:22 distortion of the trailing edge of the pulse was the only indica- 
tion of a second path. For the pictures taken at 11:29 and 11:44 the 
path difference is sufficiently small that there is almost complete can- 


cellation of pulses, only the leading portion of each pulse being present. 
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On the 11:44 picture there is Just a trace of a second pulse. The next 
set of pictures (Fig. 6) were taken a little over an hour later on the same 
night and show about the same conditions, that is, pulse amplitude and 
shape change and other evidence of the presence of a sec ond pulse delayed 
about 2 to 3 millimicroseconds. 

On the night of October 2, fading, which was apparently due to trans 
mission by way of two paths with little path difference, was observed 
Some of the results are shown on Fig. 7. At 7:49 two distinct pulses are 
evident, there being 6 millimicroseconds between their peaks (ne might 
conclude from this that there was a second path about 6 feet longer than 
the main path but the total length of the disturbance along the base lin 


and the shapes of the pulses indicate that the actual path difference was 


about 2 to 3 feet. 

Apparently we had here two pulses of r-f energy overlapping in tim 
and involving a large number of frequencies These pulses are capable 
of interfering with each other in a rather complicated manner, it being 
possible for some frequencies to add and others to cancel at the same time 
depending upon their relative phases. Phase relationships ot course de 
pend upon frequency and path length differences. As a result pulses may 
be distorted and have their peaks shifted about by a considerable amount 
We must, therefore, realize that the first picture of Fig. 7 does not really 
represent two distinct pulses as appears to be the case, but actually shows 
the resultant interference pattern of two overlapping pulses. Since 
a change of path difference of only about one and one-half inches 
enough to produce a LS8O° change in relative phase at 4000 me, it is not 














feceived pulses disturbed transmission 
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Fig. 7 —Received pulses—disturbed transmission 


Surprising threat pulse shapes and amplitudes change very rapidly 


i 


Looking again at the photograph, Fig. 7, we see that at 7:54-05 there 


Vas a complete fade as far as our system is concerned, To produce this 
degree of cancellation the path difference must have been very small 
though still sufficient to give a relative phase angle of TSO° at the radio 


frequencies involved. At 8:00-15 and 8:08-00 pulse distortion is the most 
noticeable effeet of the “fading,” the pulses being considerably shorter 
than their normal value. Pietures, not shown here, taken between 7:54 
7, 


nd S:O8 show definite evidence of two-path transmission with a path 


difference ot 2 to 3 


feet: therefore the pulses of 8:00-15 and 8:08-00 
ie probably also the result of two-path effects. 

Phe first two pictures ol Fig. S show another form of pulse distortion 
observed on a number of occasions. Here the pulse is flattened out on 
top probably due to energy coming in over a second path differing in 
length by only one or two feet from the main path. Each time this type 
of pulse was observed a check was made to be sure that the flattening 
is not due to overload in our equipment. The pictures presented up to 


now have all shown comparatively slow changes of conditions. Very 


rapid changes were, however, quite common. In many cases pulse shape 


ol 


umphtude changed considerably during the 5 to 15 second exposure 
time ordinarily used. The picture taken at 2:20-45 A.M. on August 27 


sone example of such a rapid change, there being two definite sets of 


mditions shown on the one photograph. The remaining picture on Fig. 
the pulses used for obtaining time calibration of the system. 
vere spaced 9 millimicroseconds apart and by adjust- 


Lists 
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ing sweep expansion so that succeeding pulses fell on proper parts of 


the scale and by keeping this expansion Constant, it Was possible to ol 


tain a calibration 


TWO-PATH SIMULATOR 


As an aid to interpreting the results obtained from the above experi 
ment, particularly when the two pulses overlap and interfere, a circuit 
Was set up in the laboratory to simulate two-path transmission, The 
equipment, as shown on Fig. 9, consisted of a wave guide hybrid junction 
with the r-f pulse energy being fed into the E plane arm. To each side 
arm Was connected a variable attenuator in series with a few feet of wave 
guide fitted with a short circuiting plunger. Waves reflected from these 
two plungers recombine in the H plane arm where the detector is located 
There are two separate paths through the hybrid as follows: (1) Input, 
side arm A, reflecting plunger A, side arm A to output. (2) Input, side 
arm B, reflecting plunger B, side arm B to output. By adjusting the 
attenuator in either branch the amplitude of the signal transmitted by 
way of that branch could be adjusted. In the same way by adjusting 
the position of the reflecting plunger in either branch the distance 
traveled by a signal in traversing that branch could be varied. 

If the path lengths were made the same and the amplitudes adjusted 
to be equal there would be perfect cancellation due to a phase turn-over 
in the hybrid junction and hence no output from the detector. If one 
plunger were now left fixed in the above position and the other moved by 


au quarter wavelength (to produce a total shift of half wavelength or SO 


AUG 271950 2120 — 45 AM. 


Fig. 8 teceived pulses and calibrating pulses 
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tput signal would be at maximum amplitude due to addition of 
oming from the two paths. This amplitude is, of course, twice 

of the signal from one branch alone. In the experiment the plunger 

one branch was left fixed and the attenuator in that branch left set 

ut zero. The path through this branch then represented the normal 


insmission path for an actual svstem. The path through the other 


branch could be made to correspond to spurious paths having different 
nts of delay and attenuation simply by adjusting the position ol the 
eting plunger and the setting of the attenuator. A series of photo 
iphs were taken of pulses resulting from these different amounts of 
delay and attentuation. 
Phe first three pictures of Fig. 10 were taken with the path lengths 
etly equal. When the amplitudes were also equal there was complete 
ellation. As the signal in one branch was attenuated the amplitude 
the resultant pulse increased until it beeame equal to that of the 
rinal pulse as shown in the third picture. Increasing one path by one- 
lt wavelength brought the signals from the two branches into phase 
d they added up to double amplitude as seen in the fourth picture. It 
hould be pointed out that although in our experiment we changed delay 
by O86 millimicroseconds in going trom the first minimum to the first 


maXimum, in free space a change @f delay of only 0.125 millimicroseconds 














two-path transmissior 
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Fig. 10—Simulated two path transmission 


would he required The discrepancy lies Ith the large ratio between the 
phase velocity ana group velocity in the wave guide used whereas in 
free space this ratio is, of course, equal to unity. In free space the amount 
of delay required to go from a maximum to a minimum signal corresponds 
to a change of path difference of onlv about one and one-half inches 
With only this slight shift required to change conditions from those shown 
by the first picture of Fig. 10 to those shown by the last, it is not at all 
surprising that the received signal is verv unstable during time of multi 


path transmission. 


hig. 1] shows the effect ot changing relative phisise sin steps v1 


keeping the amplitudes equal to each other. [It is seen that even with 


carriers in direct Opposition cancellation is tar from complete due to 
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Fig. 12--Simulated two-path transmission 


relative delay between the two component pulses. Flat topped pulses 
seem to be characteristic of conditions where the two carriers are in phase 
quadrature and about equal in amplitude. 

hig. 12 shows a set of conditions with a constant delay difference of 2 
millimicroseconds (corresponding to a path difference in free space ot 
about 2 feet). For the pulses shown on Fig. 13 there was a constant delay 
difference of 7.34 millimicroseconds, enough to provide complete separa- 
tion of the pulses. The carriers were in phase opposition but with this 


amount of separation there is no overlap of pulses and the results would 


have been the same if the phase had been aiding. Any increase of delay 














Simulated two-path transmission. 
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beyond this point results only in moving the pulses farther apart and has 
no effect upon pulse shape or amplitude. 

The experimental set up just described proved to be somewhat un- 
satisfactory since it was not possible with it to produce phase opposition 
between the two carriers without having zero delay difference between the 
two paths or a difference of at least 0.66 millimicroseconds. For the length 
of pulse used this latter amount of delay difference is sufficient to prevent 
anything like complete cancellation of pulses. In fact the amplitudes of 
the two resultant peaks are only about 12 db below the peak amplitude 
of the original pulse. From this we know that for the natural path any 
fade which appeared to be complete must have resulted from path dlif- 
ferences of less than 0.66 feet, in fact from differences of less than about 
one-half foot. 


SUMMARY 

The pulse experiment results indicate that over one particular path 
at least there is, at times, transmission of microwaves by at least two, and 
probably more than two, paths. Path differences involved are from a 
fraction of a foot up to about seven feet, differences of less than about 
three feet being the most common. These results agree with those ob- 
tained by other methods. These multi-path effects result in bad distortion 
of very short pulses and even in the presence of entirely separate spurious 
pulses. These effects put a definite lower limit on pulselength and spacing 
between pulses in a pulse transmission system. The limit depends upon 
the amount of distortion which can be tolerated and also upon the per- 
centage of time such distortion can be accepted. No statistical data were 
recorded, 

With the laboratory equipment for simulating transmission over two 
paths, many of the waveforms obtained over the natural path could be 


duplicated. There were times, however, when the waveforms received by 


way of the natural path were too complicated to be explained by trans- 
mission by as few at two paths. 
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Properties of lonic Bombarded Silicon 


BY RUSSELL s. OHL 


Manuscript received August 23, 1951 


T} paper de al with anew and Very ante re sting technique bi which the 
properties ot silicon surtaces are altered very materially by bombardment 
with wns of sia h QaSCS As h jdroge nN, hie lium, nifroge nN and argon. The change 
me orectitying prope rties has hee Hl of Sper ral rnterest hut the re have been Con- 

dered also changes in the structural features of the material itself. The 
f Te fx af hombardment on the rectifying prope ries are illustrated by a series 
of characteristic curves systematically arranged to br ing out the effects of 
hie eral variables of ¢ rperiment such, for example, as ion velocity, in- 
tensity of bombarding current, length of time of hombardment, hind of Jas, 
and the t m perature ol the Specimen during bombardment. The effect of 
hombardment on materials contaminated with mmpurities ’s also illustrated. 
It of part cular pra heal importance that silicon contaminated with boron 
lo the porn where it shows relatively littl ectificalion can be modified by 


] / 


hombardment to make wt even better than most unbombarded materials. 


Some venrs ago, the writer discovered that the electrical properties ot 
silicon surfaces could be greatly modified by bombardment with positive 
10 Phe ions in question were generated in a low pressure discharge in 
ome gas, like hydrogen, helium or nitrogen, and after passing through a 
pertorate | cathode were aecelerated toa suitable velocity before imping- 
ng on the surface to be treated. This scheme may be contrasted with 
other methods subsequently reported for treating germanium! in which 
high-velocity lons were derived from radioactive sources. Preliminary 
results of the present research were described in a paper entitled Silicon 
Fransistors, by W. J. Pietenpol and the writer, presented at an Elec- 
tronics Conference held at the University of Michigan, June 22, 1950. 
Since that time exploration has continued with a view both to learning 


vhout basi prin iples sand about possible practical applications. 


/ Vote -Since the resurgence of interest in pomnt-contact rectifiers, 
research has been carried on into the ch irneteristics of silicon and 
Z u Phe author of this paper was a pioneer in this new field of study, as 
rt ed. for example, | Patent No. 2.378.944, applied for on July 26, 1939, and 
ent No. 2.402.839, applied for on Mareh 27, 1941. More recent work has been 
‘ rge number of text books and technical papers such as Electrons 
i aT Cond <n) ny Willi im Shockley, D Van Nostrand, 1950, and 
pers by Lark-Horowitz published mostly in Phvsieal Review. The 
es edin the secompanving paper is a continuation of this long research 

Ir {1 Pearson, 2?) he 80, Dec. 1950 
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PROPERTIES OF IONIC BOMBARDED SILICON LOD 


The present paper gives the results of some more recent experiments 
made with improved equipment. Also described briefly are some related 
experiments in which silicon is bombarded with alpha particles derived 
from radioactive polonium. The overall results of this work indicate 
rather clearly that with suitable variations of bombarding voltage, target 
temperature and time of exposure as well as impurity content in the base 
material, it is possible to prepare to specification silicon surfaces having 
a wide range of properties. From the materials so treated it has been 
possible to construct improved forms of signal rectifiers, harmonic gener- 
ators, transistors, modulators, gating devices and also photo-electric 
cells. It is particularly significant that the voltage range over which these 
newer devices can be operated has been greatly extended, thus making 
them useful in places not previously regarded as possible. Since these 
new surfaces appear to be readily reproducible in large numbers and 
since they are electrically tough, chemically stable and show no unsatis 
factory temperature or aging effects, it would appear that bombarding 
techniques should have considerable practical value. 

This paper is concerned mainly with the practical aspect of ion bom 
bardment of silicon, namely its effect on the voltage current character 
istics at low trequencies. Equally important, perhaps, are its theoretical 
aspects, particularly with regard to the interpretation of the rather 


yronounced changes in the properties in light of presently-accepted views 
> : 


of solid-state physics. These aspects are not covered in this paper 


METHOD 


The bombardment process referred to above consists of exposing the 
silicon surface to ions that have previously been accelerated to energies 
in the range trom about 100 electron volts to about 30 kilo-electron-volts. 
A recent setup is illustrated in Fig. 1. The electrons from a tungsten 
cathode are accelerated toward a grid which is at a positive potential with 
respect to the cathode. Many of the electrons pass a short distance 
beyond the grid and return for ultimate capture. Ionization due mainly 
to the impacts of electrons with gas molecules takes place in this turn- 
around region, producing amongst other things positive particles. Elec 
trodes are so proportioned that this ionization is fairly uniform over the 
grid area. 

The silicon specimen to be bombarded is made negative with respect 
to the filament. This accelerates the positively charged particles toward 
the target. The latter rests on a graphite plate heated by a coil below, 
carrying high-frequency currents. A thermocouple with suitable con 


nections to the exterior makes possible an adequate measurement of 
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temperature. The apparatus will accommodate circular surfaces as large 
as 1} inches diameter. The gases from which ions are derived are ad- 
mitted through the gas inlet. Thus far experiments have been made 
with hydrogen, helium, nitrogen and argon. The bombarding voltages 
have as already noted, been varied from 100 to 30,000 volts and the 
surface temperatures have been varied from about 20°C to 400°C. The 
effects of these several variables will be discussed more fully below. 


SAMPLE PREPARATION 


The material to be bombarded is usually prepared in batches of about 
300 grams in fused siliea crucibles roughly evlindrical in shape. After 
solidifving, the cast is ground to a cylinder approximately 1} inches 
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diameter. This has the effect of removing some of the contaminating 
impurities derived from the crucible as well as providing samples of 
convenient size. This 13-inch cylinder is then sliced transversely into 


thin wafers which subsequently are polished on one side. Except as other- 


wise noted the material covered by this paper had an impurity content 
of about 0.1 per cent. The exception will be found in the data of column 
(a) of Fig. 8. 


BOMBARDMENT PROCEDURE 


The wafer, as prepared above, is placed in the bombarding apparatus 
with the rough face contacting the graphite support. The vacuum cham- 
ber is sealed by placing the ion generator in position and the whole 
assembly is evacuated. The sample is then heated to the proper tempera- 
ture and the desired kind of gas is admitted, the pressure being estimated 
from the ion current. When stable conditions prevail, the accelerating 
voltage is applied to the target and the bombardment is carried out for 
the proper length of time. A convenient current density is 5 micro- 
amperes per square centimeter of target area. The target area of our 
present apparatus, including the silicon and a portion of the graphite 
support, being 20 square centimeters, the ion current is generally around 
100 microamperes. The dosage is sometimes specified in microcoulombs. 

After bombardment, the sample is removed from the apparatus and 
the rough surface is covered with a thin layer of evaporated rhodium 
For most of the tests outlined below the 14-inch diameter wafers were 
cut into £-inch squares, a size convenient for testing. 


GRAPHICAL REPRESENTATION 


In considering the merits of non-linear materials such as silicon, per- 
haps the simplest and most useful characteristic is the voltage-current 
relation. If this is plotted to a linear scale, it results in a smooth curve 
of the general form shown in Fig. 2a. Specifie curves obtained in practice 
may depart widely from that shown but in general, all may be regarded 
as made up of two semiparabolas, one in the first quadrant and one in 
the third, joined by a nearly horizontal straight line. For present put 
poses, we shall further simplify this idealized characteristic by consider- 
ing it as made up of three straight lines. The first, AB, is associated with 
the reverse voltage current characteristic. The third, CD, is associated 
with the forward voltage current characteristic. These two character 
istics are joined by the nearly horizontal line, BC. The slopes of these 


three lines correspond to resistances. The section BC for example, corre 
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sponds to a region in which resistance is very high. The points B and C 
ire particularly important for they represent points of inflection where 
the resistance undergoes rapid change and the material is departing 
most markedly from Ohm’s Law. Ideally they should be sharp but in 
practice there is usually considerable curvature. Though either inflection 
point could presumably be used in detection processes, the point to the 
right of the origin is for practical reasons, usually preferred. Point B 
defines a voltage Ey at which substantial backward currents flow. It is 
referred to simply as the reverse voltage. In a similar way, point C defines 
a forward voltage E, The distance between B and C (hp | Kp) will be 
referred to as the ¢nflection interval. The ditference in these quantities 
BE, Ie) is also of interest. One-half of this voltage difference is re- 
ferred to as the self-brasing voltage. It is a significant quantity readily 
measured in practice by noting the d-c voltage across a large condenser 


placed in series with the crystal and a supply of 60 eycles AC. For de- 


tectors, point © should preferably be close to the origin and /y should 


/ 





LOG I 


lealized characteristic curves 
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be small. For certain kinds of voltage limiters, Ey should be large. In 
either case the inflection interval should be large. 

In an alternate graphical representation, see Fig. 2b, voltage-current 
data are plotted to a log-log scale. This form of representation is of 
particular value when large ranges of data are to be shown, It is also of 
value in determining the resistance (#,) at small voltages. Corresponding 
points on the two curves shown in Fig. 2 are identifiable by the letters A 
B, C, and D. Curves of both kinds are used interchangeably to show the 
effects of the several variables of the experiment. 

EFFECT OF CONTACT PRESSURE 

In point contact rectifiers,’ pressure is of considerable importance. 
Usually the best pressure is a compromise between good electrical charac- 
teristics, usually obtainable only with light pressures, and good stability 
usually obtainable with higher pressures. Experiments have been per- 
formed with a range of contact pressures both on bombarded and unbom- 
barded materials. In general, the results are highly variable, particularly 
in the case of unbombarded material. From this wide range of data, 
however, two characteristics have been selected that may be regarded 
as typical for 10-gram and 60-gram pressure. They are shown in Fig. 3 
for silicon taken from nearby portions of the same sample. Significant 
points on these several curves may be compared with their idealized 
counterparts shown in Fig. 2. Although the samples chosen show some- 
What more than the usual intrinsic resistance typical of p-type silicon, 
the effects of contact pressure are nevertheless regarded as representative. 
As indicated in Figs. 3a and 3b, the effect of increased contact pressure,' 
particularly in the case of unbombarded material, is of reducing the low 
voltage resistance, ?,, see Fig. 2b. The more desirable higher resistance 
is obtainable only with light contact, a condition unfavorable for high 
mechanical stability. In the case of bombarded material, the effect of 
contact pressure is less important. Thus it is possible in this case to 
incorporate in the design higher contact pressures and obtain thereby 
higher stabilities. For purposes of this paper a contact force of 10 grams 
has been accepted as standard. 

In addition to showing the effect of contact pressure, Fig. 3 shows some 
overall effects of bombardment. It will be noted, for example, that the 
effect of bombardment, see Fig. 3b, has been that of shifting the plots 
of Fig. 3a to the left by several orders of magnitude. Thus the resistance 
(R,,) is increased by a factor of more than 10,000. It is to be noted also 


3 Scaff and Ohl, Bell System Tech. J., 26, Jan. 1947. 


4 Realy contact force 
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Characteristic curves. (a) and (c¢) unbombarded silicon, (b) and (d) 
mibarded with 30-kv helium ions 
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from a comparison of Fig. 3a with Fig. 3b that at the one volt level, the 
ratio of forward to reverse currents for the unbombarded case is about 
twenty, whereas that for the bombarded case, is more than 10,000. At 
other levels the difference is even greater. Referring particularly to Figs. 
3e and 3d, it will be seen that one effect of bombardment is that of 
separating the two significant points of inflection B and C. That is, the 
inflection interval has been notably increased. This increase is the result 
of a small increase in the forward voltage and a very substantial increase 


in the reverse voltage. 
EFFECT OF TYPE OF GAS 


Four high purity gases were tested as lon sources, namely, hydrogen, 


helium, nitrogen and argon, having atomic weights respectively of 1, 4, 
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14 and 40. While all four gases worked very well, helium was the easiest 
to handle. In the course of the tests, identically prepared samples each 
hoinch square, taken from a high-purity silicon melt, were bombarded 
with ions formed in the particular gas under test. Particular conditions 
known to be good for producing good rectifier units were adopted as 
standard for these tests. They corresponded to a total bombarding charge 
of 600: microcoulombs per sq. em., a surface temperature of 395°C a 
contact force of about LO grams and a bombarding potential of 30 kv 

Phat the effeet of bombardment varies with different gases is seen at 

‘lance from the characteristic current voltage eurves shown in Fig. }, 
Figs. 4b to de in particular indicate that as compared with an untreated 
sumple, Fig. da, the effect of bombardment is in general that already 
noted of separating the two significant points of inflection, B and C. A 
rather substantial increase in the forward voltage appears in the case 
of argon as compared with hydrogen, helium and nitrogen. In contrast 
with a small inerease in the forward voltage resulting from the bom- 
bardment of helium, there is a very substantial increase in the backward 
voltage. Though substantial for all four gases, the effect of bombardment 
is largest for helium with progressively smaller effects noted respectively 
for argon, hydrogen and nitrogen. A particular characteristic of helium 
bombardment, as compared with argon, not readily appreciated from a 
linear seale, is shown in Fig. 4f. It will be noted that at the one volt level, 
the ratio of forward to reverse current for helium is about 130 whereas 
for argon it is about 25. At other levels the difference is even greater. 
\t the moment helium is regarded as a preferred source of ions. 

The log-log current-voltage curves show as before that the lowest 
voltage at which substantial forward currents flow occurs in helium, while 
the highest forward voltages occur for argon. In a similar way the 
voltages for substantial reverse currents are highest for helium and 
lowest for argon. The sharp break in the reverse current characteristic, 
evident in these cases, has been observed so generally that it is now 


accepted as typical of bombarded surfaces. 


EFFECTS OF TEMPERATURI 


Investigations have been made of the properties of silicon surtaces 


as atfected by the temperature at which bombardment was carried out. 


Phis has been done not only for surfaces used as rectifiers but surfaces 


used as transistors and as photo-electric cells as well. In the case of recti- 
fiers, a procedure was adopted similar to that used in the previous tests. 


Mi istirements were mi le at five different temperatures ranging trom 
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25°C to 395°C each at accelerating voltages of 1 kv. 3 kv. 10 kv and 30 


kv using helium gas. The data so obtained were useful not only for study 


ing the effect of temperature but useful in the studies of the effect ol 
ion velocity as well. The latter will be discussed in the following section 

The results of the above measurements are plotted in Fig. 5. They 
are further summarized in Fig. 6a. The latter figure, in particular, indi 
cates that as rectifiers, there is little choice of surface temperature he 
tween about 250°C and 400°C. It has been found, however, that for 
temperatures below about 250°C) the point contact seems to be more 
vulnerable to electrical shock ' 
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EFFECT OF ION VELOCITY 


Phe effect of ion velocity (bombarding voltage) has been investigated 


for several types of silicon. The effects vary with the different types. 


Typical results are those given in Fig. 5 already referred to. It will be 
noted from a comparison of the data for a particular temperature, say 
300°C, that the principal effect of increased ion velocity is that of increas- 


ing the reverse voltage. Values of these reverse voltages Ey and also the 














« 











3 3.5 40 45 5.0 


MBARDING VOLTAGE) 


a) effect of temperature and bombarding 
effect of bombarding voltage on reverse 





PROPERTIES OF IONIC BOMBARDED SILICON 
TaBLeE I—Errect or BOMBARDMENT VOLTAGE ON Ep AND 


3 KV 


Surface Temp 


395 
300 


or 


on 


140 
Mean 


Mean E's 


corresponding forward voltages Ey have been scaled from the above 
drawings and have been tabulated below. Since they vary only slightly 
with surface temperature, only their mean values are regarded as signifi- 
cant. Mean values of £’, are plotted in Fig. 6b. 

In order to isolate further the effects of bombardment we have sub- 
tracted from the mean values of /, value of Ey for untreated silicon. 


Thus the curve marked Fy represents the improvement in backward 


voltage that has accrued from bombardment alone. This is also tabu- 
lated as the mean F/’, in Table I. 


EFFECT OF TOTAL CHARGE 


Tests have been made to determine the effect of time of bombard- 
ment on the rectifying properties of silicon surfaces. In these tests, 
specimens taken from neighboring regions of the same melt were ex- 
posed for progressively longer periods all at the same bombarding po- 
tential of 30 kv and the same rate and density of application, 5 micro- 
amperes per square centimeter. Representative current-voltage charac- 
teristics are plotted in Figs. 7a and 7b for two neighboring regions. 
The results are summarized in Fig. 7¢. The latter show a rather rapid 
improvement of back voltage Ly with total charge up to perhaps 50 
microcoulombs per square centimeter.’ Thereafter the improvement is 
small. For purposes of comparison, there is plotted as a vertical line a 
value of bombarding charge that would account theoretically for one 
positive ion in each unit crystallographic cell on the surface layer. This 
suggests that when all surface cells have been penetrated by a single 
ion, no marked increase in back voltage can be effected. 


6 Specifying results in terms of microcoulombs implies that bombarding ef 


fects are independent of the rate of application. This is known to be true only 
between factor limits of } and 2 of the bombarding current 








L1G THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 








volt ive 








50 
26 


k:tfect of time 


‘ 


u 








PROPERTIES OF TONIC BOMBARDED SILICON 


EFFECT OF MATERIAL COMPOSITION 


Thus far discussions have centered around a single type of high purity 
material that was regarded as representative. [It is of interest to examine 


the effect on other materials particularly those in which impurities have 


been added. For this purpose tests were made on comparable samples 


from four sources all bombarded for two minutes with 5 microamperes 
of current and each with five representative bombarding voltages. The 
results are illustrated by the curves shown in Fig. 8. The four columns 
correspond to progressively higher percentages of impurities beginning 
with (a) on the left as a material having an impurity content believed 
to be less than O.OL per cent. The impurity content of (b) is not known 
accurately except that it lies between (a) and (¢). The material repre 
sented by column (c) was produced by adding 0.02 per cent boron® to 
« material illustrated in column (a). The last column (d) was produced 
by adding 0.1 per cent boron to the material illustrated in column (b 

It is to be noted from Fig. 8 that marked changes in the voltage 
current characteristic may be effected by bombardment for all degrees 
of the impurity content shown. It is especially interesting that in columns 
(¢) and (d) corresponding to materials contaminated with boron to the 
point where nonlinearity is almost absent, rectification can not only 
be restored but indeed the product may be made better than the best 
unbombarded material 


EFFECTS OF ALPHA-PARTICLE BOMBARDMENT 


The close relationship between helium ions such as generated above 
and alpha particles such as emanate from radioactive materials suggests 
that the latter may be used for the bombardment of silicon surfaces 
\ few experiments of this kind have been made with results that are 
not only interesting but possibly useful. For these tests, four sources of 
alpha particles were obtained. They consisted of Py inch square pieces 
of nickel on which had been plated a thin coating of polonium followed 
by a covering of gold. The initial strength was 4 millicuries per square 
centimeter. The half-life of polonium is 140 days 

The process of bombardment consisted simply ot placing the polished 
surface of a standard silicon square against the laver of gold and examin 
ing the same periodically. Tests of four samples were carried out simul 
taneously. The results are given in Table IT. The data for Sample No. | 
departs so markedly trom the mean that it mav be disregarded. Since 


6 Boron is a particulariv active agent in eflecting changes in the 


of silicon 


prope rties 
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the data for self-bias are the result of a direct measurement while those 
for forward and reverse voltage are transcribed from a cathode-ray 
plot, they are perhaps the most significant. Samples | to 3 represent 
specimens of increasing degrees of purity. 

These alpha-particle bombardment experiments indicate rather defi- 
nitely that results may be obtained similar to those obtained from 
bombardment with gaseous ions and, like the ion bombardment, they 


tend to produce high resistance surfaces. 


TaB.e II 


Sample Number 


Days of Exposure 


Self Bias-Volts 
teverse Voltage 
Forward Voltage 
Self Bias-Volts 
Reverse Voltage 
Forward Voltage 
Self Bias-Volts 
teverse Voltage 
Forward Voltage 
Self Bias-Volts 
Reverse Voltage 
Forward Voltage 
Self Bias-Volts 
Reverse Voltage 
Forward Voltage 


MECHANICAL EFFECTS OF BOMBARDMENT 


The marked changes in the electrical properties of silicon imposed by 
bombardment strongly suggest that bombardment may also impose a 


corresponding change in the lattice structure and that this might be 


detected by suitable optical methods. Attempts were made at an early 
date to detect such changes. To this end a mask of nichrome ribbon 5 
mils wide and | mil thick was laid over a sample of silicon during bom- 
bardment. An optical examination of the surface showed that after 
bombardment in the case of helium the surface on either side of the 
mask was elevated whereas in the case of argon it was depressed. This 
result has since been confirmed by one of the authors’s colleagues, 
Dr. F. W. Reynolds, who has found that in cases of prolonged bom- 
bardment by helium the adjacent surfaces may be elevated by as much 
as 225 Angstroms’ while in the cases of prolonged bombardment by 


7One Angstrom is 107% em 





120 PHE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


gon the adjacent surfaces may be depressed by as much as 130 Ang 


stroms. Further investigation of this phenomenon is under way. 


PABILITY OF BOMBARDED SURFACES 


No extended test has yel been made of the stability of bombarded 
urfaces but results extending over more than two years are encouraging 
Similarly, rectifiers for the millimeter wavelength range, mounted with- 
out the usual protective impregnation, show little or no change at the 
nd of a veu 

In few stances bombarded surlaces have been subjected to rather 
evere tests with results that suggest that under normal conditions they 
may be even more stable than surfaces activated by more conventional 
methods. For example, surfaces contaminated while cutting or while 
cementing them to their mountings have subsequently been cleaned 
vith solvents such as aleohol and are substantially the same before and 
ifter treatment. In other cases, they have been heated in a flame to 
soldering temperatures with no appreciable effects. Even in the very 
evere case where the bombarded piece was heated to a cherry red and 
the superficially oxidized layer was removed with hydrofluoric acid the 
effects of bombardment were still evident. There was, however, con- 
siderable reduction in the tolerable reverse voltage. There is nothing 
nh our experience to date to suggest that bombarded surfaces treated 
in accordance with the simple straightforward methods outlined above, 


ire in any Wise temporary in characte! 


CONCLUSIONS 


Phe experiments reported above have shown that rather pronounced 
inges in the electrical properties of silicon may be produced merely 
bombarding the polished surface with positive ions. The ratio ot 
id to reverse currents, for example, which for the usual untreated 
>seldom more than a few hundred, can be made more than 10,000. 
experiments show that the effect depends to some extent on the type 
ion gas used, helium being a preferred medium. The effect depends 
soon the volocity of the bombarding particles, the total bombarding 
harge and to a lesser extent on the temperature of the specimen during 
bombardment. Crood results are obtained from bombarding potentials 


Wsaok vith eurrent densities oto microamperes per square centimeter 


for periods of one or two minutes. The temperature should preferable 


ibout 300 %¢ 


Ordinarily the properties of silicon are materially affected by impurity 
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content. In the case of bombarded silicon the effect is much less. More 
particularly it is possible to contaminate silicon with impurities such 


as boron to the point where its rectifying properties are almost com 


pletely lost and by bombardment it is possible to convert the crystal 


into a very useful rectifier. It is possible to produce results similar to 
the above by exposing the crystal to radioactive polonium. Bombarded 
materials appear to be relatively stable 


The writer wishes to express his appreciation of the encouragement 
and help of Dr. G. C. Southworth in the preparation of this paper, to 
A. J. Mohr, Jr., for his able assistance in the experimental work, and 
to numerous associates in Bell Telephone Laboratories for their assist 
ance in preparing materials and in making special tests for which the 
author was not adequately equipped. 
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Since the mechanical prope rties of solid polyme r materials are large ly dle- 
pendent on the motions that segments of the polymer chains can undergo, 
to understand these properties one must use measuring techniques which can 
determine these motions. One of the most promising methods is to measure 
the reaction of polymer materials to longitudinal and shear waves over a 
frequency spectrum wide enough to determine the relaxation frequencies duc 
to thermal motions of the principle elements of the chain. The presence of 
relaxations is indicated by a dispersion in the velocity and attenuation 
constants of the material, or a dispersion in the characteristic impedance 
of the material if the attenuation ts too high to allow velocity measurements. 
1 number of different types of measuring methods are described in this 
paper which make possible propagation and impedance measurements not 
only in solid polymers, but also in liquid polymers and in solutions of 
polymer molecules int jpical solvents. 

When these techniques are applied to long chain polymers in dilute solu- 
tions, the three relaxations observed correspond to motions occurring in isolated 
molecules since as the dilution increases, the molecules seldom touch. The 
lowest relaxation corresponds to a configurational relaxation of the molecule 
as a whole, the highest relaxation corresponds to the twisting of the shortest 
segment containing about 40 repeating units--while the intermediate re- 
laxation corresponds to a transient entanglement of chain segments, All 
three types of relaxations are present in pure polymer liquids but are spread 
out over a frequency range due to the pe rturbing effect of near ne tghbors of 
adjacent chains. The high frequency shortest chain relaxation can be traced 
in solid polymers of the linear chain type such as polyethylene and nylon 


and produces rubber-like response to mechanical shocks of very short duration. 


I. INTRODUCTION 

Phe mechanical properties of solid polymer materials are largely deter- 
mined by what motions, parts or segments of the polymer chains can 
undergo. Toughness, mechanical impact strength and ultimate elongation 
depend on the facility with which the polymer molecule can be displaced. 


122 
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If only a small motion of the polymer chain can occur within the time of 
the measurement, the material has high elastic stiffness coefficients and 
acts similar to a rigid solid. On the other hand, if significant segments 
of the polymer chain can move at the frequency of measurement, the 
elastic stiffness is much lower and rubber-like behavior results. An inter- 


mediate case, which occurs when the significant motion of the polymer 


molecule is near the relaxation time at the frequency of measurement, is 
that of a damping material such as butyl rubber. Even “long time” 
qualities of plastics such as creep, stress relaxation and recovery depend 
on the integrated displacements of rapidly oscillating segments of the 
chain. 

One of the most promising methods for investigating these motions is 
to determine the reaction of mechanical waves on the polymer materials 
over a wide spectrum of wavelengths, eventually going to frequencies 
comparable with those of thermal vibrations of significant groups or 
segments in the macromolecules. 

If one wishes to understand the origins of these motions it is necessary 
to measure the molecules in the form of liquids or solutions since then 
the segments of the molecule are less restrained by their neighbors and 
can perform all the possible vibrations. Polymer liquids are also interest- 
ing in themselves as sources of damping material. To apply these results 
to rubbers and solid materials, one then has to measure the modifications 
of the polymer chain motion caused by the close approach of near neigh- 
bors, by measuring the mechanical properties of these materials. 

By using different types of techniques, these processes can be applied 
to molecules in solution, to liquid polymers and to solid polymers. The 
principal types of methods used for liquids are the torsional crystal, the 
torsional wave propagation system and the shear wave reflectance method, 
all of which are described in Section IT. For solids an optical method and 
an ultrasonic method are described in Section V. All of these methods 
involve displacements of 10~* em or less so that non-linear effects are 
negligible. 

All of these methods depend on setting up shear or longitudinal waves 
in the medium and observing either the velocity and attenuation of the 
wave, or the reaction of the medium back on the properties of the 
transducer. If the attenuation of a wave in the medium under consid- 
eration is low enough to permit the wave parameters, i.e., the velocity 
and attenuation per wavelength to be determined, the relaxation of 
some significant part of the polymer molecule is determined by the dis- 
persion of the wave properties which occur, as shown by Fig. 1A, in the 


form of an increase in velocity and a maximum in the attenuation per 





RELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 





| ATTENUATION | 


+ 
| 
| 
| 


| 


HALF WIDTH 


| 


LOG jg OF SOUND FREQUENCY, f 


ATTENUATION IN NEPERS PER WAVELENGTH => 











Fig. 1A—Velocity and attenuation for a medium with one relaxation frequency 


wavelength curve. If the variation in this relaxation mechanism is 
studied as a function of temperature and chain length, the type of seg- 
ment may be determined. If, however, the attenuation of the medium 
is so high that its wave properties cannot be determined, some informa 
tion can still be obtained by determining the loading, or mechanical 
impedance, that such a wave exerts on the driving crystal or trans- 
ducer. If all the relaxations occur in the stress-strain relation, it can be 
shown that there is a reciprocal relation between the propagation con- 
stant I A + jB, and the characteristic impedance per square centi- 
meter Z) given by the equation 


ZV = (R + jX)(A + JB) = jop (1) 


where A is the attenuation and B the phase shift per centimeter, R the 
mechanical resistance and X the mechanical reactance per square centi- 
meter, w is 27 times the frequency and p the density of the medium. A 
typical two relaxation mechanism’ is shown by the curves of Fig. 1B. 
By assuming values for the stiffness and dissipation factors and fitting 
a theoretical curve to the measured values, the relaxation frequency or 
frequencies can be determined. 

All the relaxation mechanisms discussed in this paper are represented in 
terms of equivalent parallel electric circuits in which the resistance terms repre 
sent viscosities and the inverse of capacities represent shear elastic stiffnesses 


In mechanical terms these correspond to a series of Maxwell models as discussed 
in a paper by Baker and Heiss to be published in the next issue 
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The most information about the motions of isolated polymer chains 
can be obtained by investigating the properties of polymer solutions. 
This follows from the fact that in pure polymer liquids, and in solids, 
the mechanical properties are mainly determined by interactions be- 
tween chains on account of the close packing of the chains. If, however, 
one dissolves the polymer molecules in a solvent, the inter-chain and 
intra-chain reactions can be separated as the dilution increases. When 
the polymer is in the order of one percent of the solvent, the chains on 
the average touch very seldom and the mechanical properties of the 
solution are determined by the properties of single molecules. As dis 
cussed in Section III, three types of chain segment motion have been 
isolated, (l) a configurational relaxation of the chain as a whole, (2) a 
position change of the shortest segment and (3) twisting of the shortest 
chain segment. Above the frequency of relaxation of this chain segment 
the joints of the polymer molecule become frozen and the chain be 
comes very stiff. These shortest chain relaxations occur also in pure 
polymer liquids, in rubbers and in non rigid solids with linear chain 


segments such as polyethylene. In pure liquids a lower frequency quasi- 
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configurational relaxation also occurs for chain lengths greater than 60 
elements but for chain lengths less than 40 elements this type of relaxa- 
tion disappears. From the difference between the high frequency shear 
elasticities measured for polyethylene and nylon 6-6 and those measured 
for static pulls, it appears that there may be lower frequency relaxations 
in these materials as well. 


il. METHODS OF MEASUREMENT FOR SOLUTIONS AND PURE POLYMER 
LIQUIDS 


To measure the mechanical properties of such dilute solutions, shear 
waves have been used since for longitudinal waves the added stiffness 
caused by the dissolved polymers is very small compared to the stiffness 
of the solvent alone. The velocity and attenuation of a longitudinal 


wave are given by the equations 


= [x + 2m] (2) 
pv 
where \ and uw are the Lamé elastic constants, f the frequency, p the 
density, » the sound velocity, x the compressional viscosity and 7 the 
shear viscosity. Since for a one percent solution of polyisobutylene in 
evelohexane the shear elasticity does not exceed 90,000 dynes/cm’, 
whereas the value of \ is in the order of 2 X 10" dynes cm’, it is obvious 
that the longitudinal velocity would have to be measured to an accuracy 
of | part in 100,000 before the presence of polymer molecules could be 
ascertained, Attenuation measurements give some information on the 
added viscosity due to the chain molecules but since longitudinal attenu- 
ations are not easily measured below | megacycle, the most interesting 
Frequency range Is missed, 

\ pure shear wave in a viscous liquid is propagated according to the 


equation 


(3) 


where v is the transverse particle velocity, p the density, f the frequency, 
n the shear viscosity, / Vv —1 and z the distance. For typical sol- 
vents, the attenuation is so high that wave motion cannot be measured. 
However the viscous wave produces an impedance loading on a crystal 
generating such a wave which can be measured by the change in the 


resonant frequency and the change in the resistance at resonance. The 


mechanical impedance per square centimeter caused by such a viscous 
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Wave is equal to” 


Lo = % thnp (1 , 2 R uw T JX v7 (4 


lhis causes a change in resistance, and a change in frequency in a 


crystal generating a shear wave in the liquid equal to 
AR, = KiRy; Af = —K.Xy (9) 


where K, and Ky are constants of the crystal which can be obtained 
approximately from the dimensions and piezoelectric constants of the 
crystal but which are more accurately obtained by calibration in known 
liquids. The constants A, and Ay vary slightly with temperature and 
should be calibrated over a temperature range. 

The first instrument to use a vibrational method for measuring vis 
cosity was the vibrating wire method of Phillipoff.” In this method a 
wire Was vibrated in a liquid and the damping rate was used as a measure 
of the viscosity. Another method also applicable in the low frequency 
range is the transducer method of Ferry.’ In this method wires are 
vibrated by electromagnetic transducers and the resistance and reac 
tance drag on the wires are measured by the change in the electrical 
resistance and reactance of the transducer. From the constants of the 
transducer, the equivalent viscosity and stiffness of the liquid can be 
measured. 

In the medium frequency range a torsional crystal’ method was 
devised by one of the writers which has been applied in the frequency 
range from 10 to 150 ke. The torsional erystal is shown by Fig. 2. For 
these types of measurements the crystal usually is made of quartz with 
four electrodes of gold evaporated on the surface. Four wires are soldered 
on the surface and serve as supports as well as electrodes. The motion 
is all tangential to the surface and tests at Bell Laboratories and at the 
Franklin Institute,” where a precision study of the torsional crystal has 
been made, have shown no observable longitudinal waves from the 
crystal surface. The process of measurement consists In measuring the 


2W. P. Mason, Piezoelectric Crystals and Their Application to Ultrasonics 
D. Van Nostrand, 1950, p. 340. 

3W. Phillipoff, Physik. Zeits, 35, 1934, pp. 884-900 

‘T. L. Smith, J. D. Ferry and F. W. Schemp, ‘‘Measurement of the Mechanical 
Properties of Polymer Solutions by Electromagnetic Transducers,’’ J. App. 
Phys., 20, No. 2, Feb. 1949, pp. 144-1538. 

5 W. P. Mason, ‘‘Measurement of the Viscosity and Shear Elasticity of Liquids 
by Means of a Torsionally Vibrating Crystal,’”’ A.S.M.E., 69, May 1947, pp. 359 
367. 

éP. E. Rouse, Jr., E. D. Bailey, and J. A. Minkin, ‘‘Factors Affecting the 
Precision of Viscosity Measurements with the Torsional Crystal,’’ Laboratories 
of the Franklin Inst., Report 2048, presented to Am. Petroleum Inst., May 4, 
1950. 





128 CHE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


resonant Trequency and resonant resistance of the crystal In a& Vacuum, 
then introducing the solution to be measured, the change in the resonant 
resistance AR, and the change in resonant frequeney Af are determined 


by an electrical bridge. Several short cuts are possible if the mechanical 


impedance is not too high. By measuring the capacity at a frequency 


considerably higher than the crystal frequency, the resistance at reso- 
nance and Af can be obtained by changing the frequency and resistance 
until a balance is obtained leaving the capacity unchanged. This method 
has been used to measure viscosity, and a recent precision study at the 
Franklin Institute’ has shown that it agrees with other methods to an 
iccuracy Ol well under a per cent. 

The torsional quartz crystal has been successfully used to measure 
liquids having a viscosity up to 10 poise, but above this viscosity the 
electrical resistance gets so high that it is hard to measure it since it is 
shunted by the much smaller reactance of the static capacitance of the 


erystal. A erystal of higher electromechanical coupling such as am- 


ind SO-ke ervstal for shear viscosity and elasticity measurements 





MECHANICAL PROPERTIES OF POLYMERS 


i); 


Fig. 3—-Photograph of torsional erystal and rod 
I ; 


monium dihydrogen phosphate (ADP) will cause the electrical resistance 
component to be smaller in comparison to the reactance of the static 
capacitance and hence can be used to measure higher viscosities. How- 
ever, since wires cannot be soldered to the surface but must be glued, 
the crystal is much more fragile than quartz and its use has been aban- 
doned in favor of another method which makes use of the phase and 
attenuation change in a torsional wave in a rod caused by the surround- 
ing liquid whose properties are to be measured. 

This method, devised by one of the writers,’ consists in sending a short 
train of torsional waves, periodically repeated, down a glass or metal 
rod. As shown by Fig. 3, the torsional wave is generated by a torsional 
quartz crystal soldered or glued to the end of the rod. These waves 
travel to the free end of the rod and are reflected back to the crystal 
where they are detected, amplified, and displayed on a cathode ray 


7This method is described by H. J. MeSkimin, in a paper before the Ac- 
oustical Soc. of Am. in October, 1951. 
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oscilloscope. Echoes due to end to end reflections also appear, being 
attenuated by normal acoustic losses until they are undetectable by the 
time the next pulse is applied. 


With only air surrounding the rod, a phase reference and amplitude 


reference are obtained for the first received wave (or subsequent echoes 
if greater sensitivity is desired). The rod is then immersed a definite 
length in the liquid to be measured, as shown in Fig. 4, with a resulting 
phase retardation and amplitude reduction. These are measured by em- 
ploying the experimental circuit shown by Fig. 5. In order to allow the 








ig. 4 Photograph of complete torsional wave measuring equipment. 
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use of one crystal for both receiving and transmitting, the crystal is put 
in the bridge circuit of Fig. 5 where a resistance and capacity are used 
to balance out the transmitted pulse so that it will not overload the 
amplifier. The relatively weak voltages generated by the incoming 
acoustic Waves pass through directly. The gate circuit provides pulses 
of radio frequency voltage at repetition rates in the range of 20 to 100 
per second with a synchronizing voltage supplied to the oscilloscope for 
the horizontal sweep. The frequency range of the device is from 20 to 
200 ke. Both glass and nickel-iron rods were used, the latter having a 
very low frequency-temperature coefficient. With a 100-ke quartz 
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Fig. 5—Kxperimental pulsing circuit for measuring torsional impedance of 
liquids. 
torsional crystal, a rod length of 21 inches and diameter of 0.2 inch were 
used. The entire crystal-rod assembly is placed inside a glass tempera- 
ture control unit, as shown by Fig. 4, through which water can be cir- 
culated to provide temperatures in the range O°C to 80°C. The test 
liquid is placed either directly into the inner bore of this water jacket, 
or in another tube which can be inserted from the bottom to surround 
the rod up to a fixed mark. 

In use, both phase and attenuator settings were adjusted to balance 
the first received pulse against the continuous wave component passing 
through the attenuator. Cancellation for the duration of the pulse was 
visually indicated on the oscilloscope. A plot of balance phase and level 
is made as a function of the temperature. When the liquid is introduced 
an attenuation change AA and a phase change AB are required to 
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re-balance the circuit. These are measured by the amount of attenuation 
in nepers (1 neper = 8.68 db) and the number of radians phase shift 
required to re-establish a balance. An alternate method of measuring 
phase shift is to measure the change in frequency required to re-estab- 
lish balance. If this method is used the phase shift change of the overall 
circuit: with Trequency has to be calibrated for the uncovered rod by 





Fig. 6--High frequency shear reflection method for measuring shear imped 
inces of liquids 


noting the frequencies for which 360° phase shifts (as measured by 
balance) occur in the circuit. 

It is shown in the appendix that the torsional impedance of the 
liquid per square centimeter is given by 


Z (a) 5 [AA + jAB] (6) 


where p is the density, and vp the sound velocity in the rod, a is the 
radius and ¢ the covered length of the rod and AA and AB are respec- 
tively the change in attenuation in nepers and the change in phase shift 
in radians to re-establish balance. If a very viscous liquid is used it may 
be necessary to correct for the fact that the torsional impedance may 
differ from the plane wave impedance as discussed in the appendix. 

This device can measure liquids having dynamic viscosities from 10 
poise to 1,000 poise with an accuracy of the order of 10 per cent. The 
frequency range covered may be from 20 to 200 ke depending on the 
size of the erystal used to drive the rod. Hence it supplements the 
torsional erystal method for very viscous liquids. 

\t frequencies above 500 ke, the torsional crystal becomes too small 
to be used practically and recourse is had to a high frequency pulsing 
method.’ As shown by Fig. 6, shear waves are set up in a fused quartz 

*W.P. Mason, W. O. Baker, H. J. MeSkimin and J. H. Heiss, ‘‘Measurements 

Shear Elasticity and Viscosity of Liquids by Means of Ultrasonic Shear 
‘hys. Rev., 75, No. 6, March 15, 1949, pp. 936-946. See also H. T. O'Neil, 


ind Reflection of Plane Shear Waves in Viscoelastic Media,’’ Phys 
o. 6, March 15, 1949. pp. 928 936 
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rod by means of a Y-cut or AT cut crystals soldered to a silver paste 
layer baked on the fused quartz surface. The particle motion of the 
shear wave is parallel to the large reflecting surface and hence only 
shear waves are reflected from this surface. These impinge on a second 
shear crystal which is connected to an amplifier and oscillograph. Since 
the attenuation in fused quartz is so low, a long series of reflected pulses 
appear on the oscillograph. When a liquid, whose shear properties are 
to be measured, is placed on the fused quartz surface, this causes a 
change in the amplitude and phase of the reflected wave. By using the 
balance method shown by Fig. 7, in which two identical fused quartz 
rods are used, one of which has a liquid layer and the other does not, 
and by using a phase shifting network and an attenuator to balance out 
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Fig. 7—Method for obtaining resistance and reactance terms for high fre 
quency shear reflection method 


pulses, the shear impedance of the liquid can be determined. If R is the 


loss per reflection expressed as a current ratio, @ the change in phase 


angle required to rebalance the circuit and ¢ the angle between the wave 
normal and the reflecting surface, it can be shown® that the shear im- 
pedance of the liquid is 


‘ — ; 1 — R° + 2jR sin 6 
Zu = Ry + jx u ie COS ¢@ | “8 °° u pa 
L1 + R® + 2R cos 0 
where Z, is the impedance pv for shear waves in the quartz. This is 
equal to 


ZL, = 2:20 X 3.76: X 10 8.27 X 10° mechanical ohms 


Since this impedance is much larger than that of the liquids that are to 
be measured, the sensitivity is increased by making ¢ large. In practice 
g Was taken as 80°. This method is applicable from 3 me up to 100 me 
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and complements the other methods. Fig. 8 shows a photograph of the 


eq upment 
lil. MEASUREMENTS OF POLYMERS IN SOLUTION 


When such methods are applied to a polymer solution, it is found 
that the resistance and reactance components are no longer equal but 
the resistance is invariably larger than the reactance. This indicates the 
presence of a shear elasticity in the solution. If the molecules have a 
single relaxation frequency, it has been found that the shear properties 
of the liquid can be represented by a stress-strain equation of the type 


Os | 

Ol | | 
a , 
Os Te) 


Na 
(3) 
18 ot 
where 7’ is the shearing stress, S the shearing strain, NA the solvent 
Viscosity, nx a molecular viscosity of some particular motion of the 
chain which disappears when the reactance of the chain stiffness uw, of 
this motion is low enough so that the motion can follow the applied 
shearing stress at the frequency of the measurement. When this type of 
mechanism is present in the liquid, it has been shown® that the impedance 
the liquid presents to the erystals is 
pusns + j| woninn + * (ne + me) (9 

Z Rau + jXm “7 = _ 
Ne T Mp/W 
Fig. & shows a plot of the resistance and reactance components of an 
assumed solution having a single relaxation frequency, and a viscosity 
30 times the solvent viscosity. At very low frequencies, the resistance 
and reactance follow that of a solution, but for frequencies comparable 
with the relaxation frequency, the resistance becomes larger than the 
reactance while for very high frequencies the two come together on a 
line determined by the solvent viscosity. If there is more than one 
relaxation frequency, the resistance and reactance may coalesce for 
several intermediate stages. A continuous distribution would give a 
definite relation between the Trequeney dependence of resistance and 
reactance 

The torsional crystal and the shear wave reflection method have 


been applied to long chains of polyisobutylene dissolved in various sol- 


~W. PRP. Mason. Piezoelectric Crystals and Theti {pplication to Ultrasonics, 


1D). van Nostrand, 1950, p. 358 
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Fig. 8—Dual reflecting block assembly for measuring high frequency shear 
impedances of liquids. 


vents with concentrations ranging from zero per cent to 10 per cent. 
Polyisobutylene is a polymer molecule having the chemical formula 


CH; H 
C C 


CH, H 


Non-planar zigzag segments can be expected in the liquid state. Fig. 10 


shows measured curves for 20 ke of the resistance and reactance for 
solutions of viscosity average molecular weight of 3,930,000 in eyclo- 
hexane. Four values of concentration were used and tWo temperatures 


were measured. For pure cyclohexane, the resistance and reactance 
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components are equal but as the percentage of polyisobutylene 
creased, the resistance increases more rapidly than the reactance 
By solving equation (9) for n4 , 7p and wes in terms of R and XY measured 
at one frequency and n, + ne» the solution viscosity, we find 
2RX (R° — X°)°/wp ; ) 
"4 = oDY ? 1B (na + mB) — 1 
wp wolna + ne) — 2RX . 
' A (10) 
(R° " A) WB 
Me = ae 
wpl(na + ne) — 2RX 
Applying these formulae to the measured results, the curves of Fig. 11 
result. The shear elasticity is directly proportional to the concentration, 
the viscosity 4 is only slightly larger than the solvent viscosity while 
the main part of the measured viscosity resides in n»g the viscosity asso- 
ciated with chain motion. Fig. 12 shows these three quantities for a one 
per cent solution measured as a function of temperature. The apparent 
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stiffness decreases with increase in temperature for a single frequency 
measurement. However, when measurements were made at 20, 40, 80 
and 150 ke it was found that the elasticity was a function of the fre- 
quency, which indicates the presence of more than one relaxation and 
complicates the determinafion of the temperature relationships. Fig. 13 
shows measurements of the shear elasticity over a frequency range from 


2.5 ke’ up to 14 megacyeles for 25°C. There is a gradual rise up to 


about 300 kilocycles after which there is a sharp break to a stiffness of 
about 90,000 dynes em fora | per cent solution of the 3,930,000 molecu- 
lar weight polymer in cyclohexane. If one analyzes the frequency varia- 
tion of the elasticity he finds that it can be fitted by three relaxation 
frequencies, one having a frequency of 230 cycles, one around 66,000 
cycles and one around 4 megacycles. A possibility exists for a fourth 
relaxation. The lowest relaxation is thought to be a configurational re- 
laxation of all the elements of the chain. The highest one appears to be 
a relaxation of the twisting motion of the smallest segment of the chain 
such as a Kuhn segment. The intermediate relaxation appears to be due 
to the motion of the ends of the smallest chain segment from one posi- 
tion of entanglement to an adjacent position. This interpretation is 
based partly on the fact that the associated viscosity of the motion is 
very similar to that for the relaxation of the twisting motion of the 
smallest chain segment and partly from data presented in the next 
section on pure polymer liquids which shows a lower frequency relaxa- 
tion agreeing in frequency assymptotically with this one, which involves 
chain motions of approximately 30 to 40 chain elements. Temperature 
variations of these elastic components show that the lowest relaxation 
mechanism has a stiffness that increases slightly with temperature in 
agreement with the kinetic theory of elasticity. The corresponding 
viscosity (q2) which comprises most of the viscosity for a solution, when 
plotted against the reciprocal of the temperature, as shown by Fig. 14, 
indicates an activation energy of 3.9 kilocalories per mole which is 
slightly higher than that of the solvent cyclohexane alone, which is 
about 3.2 kilocalories per mole. This difference of 0.7 kilocalories pre- 
sumably represents the added energy required to bend the chain in its 
configurational motion. Measurements with another chain length of 
1.18 & 10° molecular weight showed that the stiffness of the lowest (con- 
figurational) relaxation decreased from 310 to 160 indicating that the 
stiffness of this motion is approximately proportional to the square root 

1° The lowest frequency, 2.5 ke, was measured by means of a quartz crystal 


tuning fork which will be described in another paper. This instrument makes 
possible the direct measurement of configurational elasticities 
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of the molecular weight. The viscosity decreased by a factor of 6.25 
and consequently the relaxation frequency increases from 230 to 660 
evcles 

The second “entanglement” relaxation has a stiffness of about 3100 
dynes,/cm’ for a 1 gram per ce solution of the 3.93 xX 10° molecular 
weight solution and about 2650 dynes /em* for the 1.18 X 10° molecular 
weight solution. The variation with temperature, if any, is small. The 
corresponding viscosities n; for the two solutions are nearly equal as 
shown by Fig. 14, and have an activation energy of 4.25 kilocalories 
per mole. The fiual high frequency “short segment” relaxation has a 
high stiffness of 83,000 dynes em’ fora | per cent solution of 3.93 X 10' 
molecular weight. The corresponding viscosities for the two solutions 
shown by Fig. 14 have nearly identical values and an activation energy 
of 4.25 kilocalories per mole, i.e. very closely equal to the “entangle- 
ment” relaxation viscosity. 

These upper two relaxations persist in pure liquid polymers as dis- 
cussed in the next section, although they are spread out over a small 
range of relaxation time values. The highest one can be traced in meas- 
urements of mechanical properties of solid plastics such as polyethylene 
and nylon which indicates that these materials should have rubber like 
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response to mechanical shocks of very short duration. The lowest fre- 
quency configurational relaxation is spread over a wide spectrum of 


relaxation times in pure liquids. 


Measurements of these and other chains in various solvents have 
also been made and the results are discussed, from a chemical point of 
view, in a companion paper by W. O. Baker and J. H. Heiss. It is shown 


that the stiffnesses vary with the polymer chain and the solvent used. 


IV. MEASUREMENTS OF PURE LIQUID POLYMERS 


A. Shear Wave Measurements in Liquid Polymers 


Similar shear wave measurements have been made for pure polyiso- 
butylene liquids of molecular weights from 904 to 10,380 (i.e. from 16 
chain elements to 186 chain elements), by the techniques described in 
Section IL. Some of these results have been discussed in reference (8) 
but the much more comprehensive measurements made since require 
some revisions of the original conclusions. 

The easiest data to interpret are the high-frequency data obtained by 
the shear wave reflectance method. The data of Table I give measure- 
ments of 8 liquids varying in average molecular weight from 900 to 
10,380, at three frequencies and four temperatures. If we plot for ex- 
ample the Maxwell shear stiffness and viscosity for the three frequencies 
and for 25°C as a function of the number of chain elements (here a chain 
element is taken as two adjacent carbon atoms one of which has two 
methyl groups attached and the other two hydrogens) the 4.5-me meas- 
urements are shown by the triangles of Fig. 15. The 14 megacycle 
measurements are shown by the circles and the 24-me measurements by 
the squares, 

An attempt was made to fit these measurements with a two relaxa- 
tion mechanism shown by the figure with two stiffnesses which are taken 
to be independent of the molecular weight and equal respectively to 
1.2 & 10° dynes ‘em’ and 6 & 10° dynes ‘em’*. The best fit is obtained by 
taking the two viscosities ; and m equal and these are adjusted for the 


different molecular weights in such a manner as to best fit the experi- 


mental curve. A fair agreement is obtained except for the range from 60 
to 90 chain elements where the two relaxation model gives too rapid an 
increase of stiffness with increase in the number of chain elements and 
at the high molecular weight viscosity range where the viscosity shows 
a dispersion in values but the model does not. The sum of the two vis- 

These results on the mechanical impedance of long chain molecules in sol 


vents have been presented at the XIIth International Congress of Pure and 
Applied Chemistry by W. O. Baker, W. P. Mason and J. H. Heiss, Sept. 13, 1951 
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cosities m and m2 assumed as a function of molecular weight is shown by 
Fig. 16. The log of the viscosity starts proportional to the molecular 
weight but above a molecular weight of 2,400 the increase is very slow 
and becomes asymptotic to a value of 240 poises. An equation which 
fits the increase in viscosity with molecular weight is 


> 11-8 tanh Z/2370 
7 > Ke 


where Z is the molecular weight. The solid line shows a plot of this curve 
and the circles are the assumed values to obtain a best fit to the meas- 
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ured values. This equation indicates that when Z = 2,370 or 43 chain 
elements the viscosity increases only a small amount more by a chain 
articulation effect and hence in this high frequency range we are dealing 
with a chain length of about 40 elements or 80 carbon atoms. This is 
checked also by a comparison of the static and dynamic Viscosity. The 
total dynamic viscosity due to the two relaxation mechanisms compared 
to the static viscosity does not differ markedly until the number of 
chain elements is more than 40. Above this value other motions than 
that of the shortest chain segment can take place and can add to the 
dynamic viscosity. The static viscosity fits an equation of the same sort, 
but the indicated chain length for the viscous motion is about § times 
that of the shortest segment. 


When a similar process is carried out over the temperature range the 


equations of Fig. 16 are obtained. The static viscosity has an activation 


energy of 16 kilocalories per mole, while the dynamic viscosity has an 
activation energy of about 11.2 kilocalories per mole. 

The relaxation frequencies for the two components are plotted as a 
function of the number of chain segments by the solid lines of Fig. 17. 
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big. 16- Triangles are measured static viscosities and circles are dynamic 
viscosities plotted as a function of molecular weight. Solid lines are a plot of the 
equations given for static and dynamic viscosities 
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For long chain segments the values become asymptotic to 8 X 10° eycles 
and 160,000 cyeles which are not far from the two highest relaxation 
frequencies obtained from the solution measurements of Section III. 
Hence it appears likely that these relaxations are due to the “‘entangle- 
ment’? motion and the twisting motion of the shortest chain segment. 
The increased activation energy is due to the fact that more energy has 
to be applied to the chain segment to break it loose from its equilibrium 
position when it is surrounded by adjacent polyisobutylene molecules 
than when it is surrounded by cyclohexane molecules. The stiffness of 
the chain is due more to the slope of the potential well than to any in 
trinsic chain stiffness as is shown by Fig. 18, which shows the two stiff- 
nesses as a function of temperature. These values are obtained by fitting 


the 15°, 35° and 50° data in a similar manner to that used for the 25 
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6000 





a 
































6 


HAIN EGMENTS 


Fig. 17—Solid lines are mean values of relaxation frequencies for the two 
mechanisms plotted against molecular weight Dotted lines indicate limits of 
regions assumed to obtain a better fit to the measured values 
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data, and the activation energy of 11.2 kilocalories per mole is obtained 
in a similar manner. 

Due to the closeness of the surrounding polyisobutylene molecules, 
one would expect that the relaxation frequencies would not have discrete 
values but would be spread about the center value in some sort of a 
Gaussian distribution. If we approximate this by representing each 
region by two relaxations, one-half, and the other twice the frequency 
of the mean value, as shown by the dotted lines of Fig. 17, the agree- 
ment with the measured values of Fig. 15 is considerably better as 
shown by Fig. 19. A wider distribution yet is indicated. 

For molecular weights greater than 2,000 the shortest chain segment 
viscosity begins to diverge from the static viscosity indicating that 
there are other relaxations for these longer chains. Some data for the 
three longest chain polymers, F, G and H have been obtained by the 
torsional rod method and the results are given in Table Il. These data 
are plotted on Fig. 20 as a ratio of dynamic to statie viscosity plotted 
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against the frequency times the static viscosity. All the viscosity data 
can be represented within the experimental error by a single curve, but 
the stiffness curves appear to require different curves for different tem- 
peratures. On analyzing the data in terms of a distribution of relaxation 
frequencies, a single curve for all temperatures could be obtained if the 
stiffness of each mechanism were independent of the temperature and 
the relaxation frequency were inversely proportional to the static vis- 
cosity, 1.e., had an activation energy variation equal to that* for the 
static viscosity for each relaxation mechanism. This condition holds 
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quite well for frequencies much lower than the relaxation frequencies 
of the smallest chain segment, but as the frequency approaches these 
relaxation frequencies, the stiffness of these polymers increases as the 
temperature decreases. 

A fair approximation to these measured values is obtained by assum 
ing one more “‘configurational” relaxation frequency in addition to the 
two smallest segment relaxations discussed previously. Fig. 21 shows 
calculations of the ratio of dynamic to static viscosity and the shear 
stiffness for 65°C and 25°C. The lowest relaxation frequency is assumed 
to have a stiffness of 6.3 X 10° dynes/em* and a viscosity of 20 poises 
at 65°C. For 25°C the stiffness of 2 X 10’ dynes/em’ is assumed and an 
activation energy of 17.3 kilocalories gives the component a viscosity 
of 607 poises at 25°C, and a relaxation frequency of 5,250 cycles. The 
average value of 16 kilocalories for the static viscosity is a result of the 
sum of the variation due to the two components. Although the agree- 
ment can be improved by assuming distributions of relaxation fre- 


quencies centered around these three primary frequencies, there does 


’ 


not seem to be much doubt of the existence of these primary relaxation 
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Fig. 21--Caleculated values of stiffness and viscosity obtained by adding a 
single ‘“‘configurational”’ relaxation frequeney to the two short chain relaxations 
obtained from high frequency measurements. 
mechanisms which show up as discrete relaxations in long chain mole- 
cules in solution. 

Some measurements have also been made to determine the effect 
of chemical substitutions in the polymer chains. In a previous paper,* 
the high frequency properties of poly-a-methyl styrene were discussed. 
This material has the polyisobutylene chain but with one methyl re- 
placed by a phenyl so that its chain becomes 


H 


For low molecular weights this material is liquid. The shear stiffness of 


this liquid is somewhat higher than for polyisobutylene but has about 


the same change with temperature. The variation of the high-frequency 


viscosity, however, is much larger for poly-a-methyl styrene than for 
polyisobutylene, and corresponds to an activation energy of 23.6 kilo- 
calories. The relaxation region for the shortest chain motion is much 


* See footnote on page 132 
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narrower than in polvisobutylene, which correlates with the smaller 
steric hindrance. 

Measurements have also been made in the 70-ke to 140-ke range for 
two non polar liquids, polybutadiene and polypropylene and one polar 
liquid, polypropylene sebacate. The first two have the formulae 


H 
C 
H H 
polybutadiene polypropylene 
while the polar liquid, polypropylene sebacate, has the formulae 


HH H H CH; H 


Coton kin 


HHH H UH H H 


The measured results are given by Table III. These data are plotted 
on Fig. 22 as a function of the product of frequency times static viscosity. 
For comparison the data for polyisobutylene polymer F is also plotted. 
By extrapolating to low values of the product-frequency times static 
viscosity it is seen that the low frequency ‘‘quasi-configurational”’ stiff- 
nesses of these liquids run from 10° to 1.5 X 10’ dynes per square centi- 
meter. Polyisobutvlene has the greatest stiffness for any value of fre- 
quency times viscosity while polybutadiene has the least. This reflects 
the greatest steric hindrance of polyisobutylene and the smallest for 
polybutadiene which has only hydrogens connected to the carbon chain 
atoms. Another consequence of the larger steric hindrance of the CH, 
groups of polyisobutylene is that the viscosity associated with the short- 
est chain segment motion is largest for polyisobutylene and smallest 
for polybutadiene as can be seen from the ratio of dynamie to static 
viscosities for high values of frequency times static viscosity. 

The activation energy for static viscosity are for polypropylene, poly 
propylene sebacate and polybutadiene respectively 21.2, 12 and 8 kilo- 
calories compared to 16 kilocalories per mole for polyisobutylene. The 
high-frequency activation energies as determined by the dynamic meas 
urements are respectively 11.8, 4.6 and 1.4 kilocalories for polypropylene, 
polypropylene sebacate and polybutadiene. The differences between the 
activation energy for static flow and that for dynamic flow are respec 
tively 9.6, 7.4 and 6.6 kilocalories, which values are all higher than 





s]uUsWIT 


PLO 


ce" 0 


Z4tT'O. 01 
€zOl'O 201 
cr 0 6990°0 | OL 
IS [1t0°O OT 
9 6FZ0'O 
19100 
FOTO O Ol 


SJUIUIIINSPIPY IY 


aywonqgeg sual AdoidAjog 


COFO , Ol 
COCO) Ol 
601 
sOl 
Ol 


Cony) 
Cll 
ENCO 

EO 
6670 
Q610 


CO x ;° 12 OOL 
SCI 009 
ct 00z.‘ 
GES ; 006 

& ’ 4 OOO 


auajAdoidAjog 


19°0 oS 01 
P 


Sb 0] 
tt 1800 Ol 
ar C090 O 

tt 


QY 


suVIpBINg A[Od 





MECHANICAL PROPERTIES OF POLYMERS 











IN DYNES 


CENTIMETER 


STIFFNESS 





PER SQUARE 


SHEAR 














OSITY 


ViSC 


Ps $ 
, 
t¢ POLY PROPYLENE 
q a: SEBACATE 
POLY BUTADIENE +F-+ T 


STATIC 





0.04 + | i + | ee Fee! 
-+7 POLYPROPYLENE 
0.02 


0.01 
0.008 + - + ° + + + + + 
0.006 

10? 2 4 6 10° 2 4 6 10” 2 4 6 108 2 4 6 
PRODUCT OF FREQUENCY TIMES STATIC VISCOSITY 

Fig. 22—Ratio of dynamic to static viscosity and the shear stiffness for four 

polymer liquids plotted against product of frequency and static viscosity 


























2) 
he 
Y 
= 
< 
2 
= 
a 
uw 
Oo 
Oo 
fs 
< 
84 


the 4.8 kilocalories for polyisobutylene. This presumably indicates that 
there is more of a difference between the viscosity flow segment and the 
shortest chain segment in these materials than in polyisobutylene. Since 
no measurements are available over a range of molecular weights, no 


direct evidence has been obtained for the various chain lengths 


B. Longitudinal Wane Measure ments in Liquid Polymers 


Since the increase in shear elasticity for the highest relaxation fre 
quency is so large, it should also appear in longitudinal wave measure- 
ments. Fig. 23 shows a calculation for the 5590 molecular weight liquid 
of the longitudinal velocity assuming that the Lamé elastic constant 


is independent of frequency and that all the variation occurs in the shear 
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lation between measured longitudinal velocity for polyisobutylene 
weight 5590 and that caleulated from shear stiffness measurements 
issuming the Lamé A elastic constant is independent of frequenes 
constant as determined by the shear measurements. The points are 
velocities measured for longitudinal waves and as can be seen, the meas- 
urements agree closely with the calculated values. A’ slightly better 
agreement would be obtained if \ increased by a small amount as the 
frequency increased. As discussed in the next section there is some 
experimental evidence for an increase in AX in nylon 6-6 and in poly- 
ethvlene 
Phe question also arises as to how much of the attenuation is due to 
shear mechanisms and how much due to pure compressional effects. 


From longitudinal velocity and attenuation measurements at 30°C for 


the polymers E, F and G of Table 1, the values of X + 24 and x + 27 


can be determined and are shown by Table IV. The values of w and 7 


can be obtained from Table T by interpolation and are given in columns 


Taste LV 
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$f and 5. Columns 6 and 7 show the values of \ and x, the compressional 
components. A definite longitudinal compressional viscosity is indica- 


ted which however is somewhat smaller than the shear viscosity n 


V. MEASUREMENTS FOR SOLID POLYMERS 


Recently two new methods have been devised for accurately measur 


ing the properties of solid plastics in the ultrasonic frequency range 
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Fig. 24—Normal loss in db per centimeter measured as a function of frequency 
for several rubbers and plastics. 
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Fig. 26—Examples of refraction and focusing effects for sound waves. 


ind these have shown relaxations in such plastics as polyethylene and 
nylon 6-6. The simplest method for measuring one of the properties for 
longitudinal waves, i.e., the attenuation, is to measure the change in 
loss between two transducers in a liquid such as water, caused by insert- 
ing a sheet of the material. This process, used during the last war, 
results in the losses in db per centimeter for several rubbers and plastics, 
shown by Fig. 24. Indications of relaxation mechanisms are given by 
the rubbers and the plastic tenite IIT which is a cellulose acetate butyrate. 


he first fairly accurate method for measuring longitudinal sound 
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velocities’” in plastics was the method of observing the focusing effect 
of a cylindrical lens made of the plastic. Sound waves can be made 
visible by the Debye-Sears technique of using a sound wave as a phase 
diffraction grating. Here light from a slit S; is made parallel by the 
lens Le and passes through the cell parallel to the wave fronts of the 
sound waves as shown by Fig. 25. The compressed parts of the medium 
retard the light waves more than the rarefied parts do and hence the 
medium acts as a phase diffraction grating. If a second slit Ss is used 
which is small enough to pass only the zero order, a light valve action 
is obtained which modulates the light according to the sound wave 
intensity. If now the lens L,; is used which focuses on the median plane 
of the tank, a picture of the sound beam is obtained as shown on Fig. 26. 
The bottom figure shows the focusing effect of a plastic and from the 
focal distance d and the radius of curvature r of the lense, one can c¢al- 
culate the velocity in a plastic compared to the velocity in the water 


r 
Mp = » l — r) ( ] | ) 
t t /( 7 


This method gives velocities good to from 2 to 5 per cent depending on 
£ £ £ 


by the formula 


the attenuation in the lens. 

G. W. Willard’’ has devised recently a more accurate method for 
measuring sound velocities as shown schematically by Fig. 27. Here a 
plastic to be measured is placed half way across the sound beam in the 
liquid and light is sent along the wave front occurring in both the plastic 
and the liquid. If the waves are in phase the retardation in the two light 
gratings, corresponding to sound propagation in both media, add up 
and for a slit selecting the zero order the darkest pattern occurs on the 
photographic plate. If the two waves are just out of phase, the retarda- 
tion is reversed in the two media and the lightest part occurs. With 


this relation it can be shown’’ that the spacing d of light and dark lines 


2W. P. Mason, Piezoelectric Crystals and Their Application to Ultrasonics, 
D. Van Nostrand, 1950, p. 404. It was used in this country by G. W. Willard as 
early as 1940. It was also used in Germany by J. Schaefer ‘‘Eine Neue Method 
Zur Messung der Ultraschallwellen in Festkorpern.”? Diss Strassburg, 1942. By 
making the front surface part of a eylinder, Schaefer also measured the shear 
velocity in a solid 

3G. W. Willard, J. Acous. Soc. Am., 28, Jan. 1951, pp. 83-94. The origin of 
this multiple path interference method goes back to the work of R. Bir (Helvetia 
Physica Physica Acta Bd 13 page 61 (1940)) who attached a piezoelectric ervstal 
to a bar with a 45° end section and set up transverse and longitudinal waves, in 
the bar. These waves produced longitudinal waves in a surrounding liquid and 
by observing the interference pattern between them, the longitudinal and shear 
constants could be determined for an isotropic medium. Willard’s method as 
described above is much more direct and is capable of higher accuracies 
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is related to the wavelength in the liquid A, and the wavelength in the 


solid A by 


(12) 


Vhis corresponds to 4a velocity in the solid compared toa velocity In 


the liquid given by 


(33) 


where fis the frequency. Fig. 28 shows a photograph of a series of lines 
ina transparent plastic and a transparent plastic in the form of a wedge. 
It is seen that bevond the edge of the plastic there is a dark interference 
band for each one in the transparent plastic. This phenomenon is caused 
by the refraction of the sound wave that has traversed the plastic and 
the dark lines are lines of equal phase of the two waves in the liquid 
Phe angle of the dark lines is half the refraction angle. Hence the velocity 


be determined by counting the number of dark bands in the 


Cah also 


liquid beyond the plastic. This makes it possible to measure the veloci 
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Fig. 28 
ind plastics 


Photographs of interference patterns from sound waves in liquids 
ties In opaque plastics. The accuracy of the method is better than 1 per 
cent if the attenuation is low enough to give a number of interference 
lines. For plastics ol high internal loss, the method becomes somewhat 
Inaccurate, 


Typical measurements using this system are shown in Table V. Small 


changes in chemical composition and plasticizer content are shown up 


as can be seen from the table. Of particular interest is the difference 
between nylon 6-6 and polyethylene. Chemically as shown by Fig. 29, 
the two are identical except for the dipoles occurring for every 6 units 
of the ethylene chain. These dipoles have the effect of bonding acd jeu ent 
lavers together and result in a higher shearing modulus 

By attaching shear vibrating ervstals toa right angled prism, as shown 
by the lower part of Fig. 28, with the direction of motion of the erystals 
parallel to the transmitting face, and setting up shear standing waves 
between the two erystals, the shear properties of the plastic can be 
measured. Longitudinal waves are generated in the liquid which inter 


fere with one another and eause dark bands perpendicular to the plastic 
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TABLE \ 


Measured longitudinal velocity and attenuation at 25°C and 2.5 me. 


Longitudinal attenuation in DB/em at 25°C and 2.5 me except as noted. 


Long Shear 

velocity velocity . 
x 107% 10-5 A DB/cm 
cm/sec cm/sec 


Density 
Brass, half hard ; ‘ “ 
Polystyrene 2.36 [8 
Plexigl is 
Fenite II, (cellulose acetate butyrate), 2% 
plasticizer 
Tenite I], 138% plasticizer 
Polyvinal formal 
Polyvinlidene chloride 
Poly N-butyl methacrylate 
Poly I-butyl methacrylate 
Neoprene 
Polyethylene 
Nylon 6-6 (3-30 megacycles) 
Nylon 6-10 (3-30 megacycles) 


Dural, 17 ST 3. 3.1 
ae | 
J 


surface. By determining the spacing of these lines the velocity of the 
shear waves can be determined. 

Another method has also been developed which is more applicable for 
high loss materials. This is a pulsing method and is a modification of the 
method proposed by one of the authors for measuring the properties 
of small crystal specimens.'* Here longitudinal or shear crystals are sol- 
dered to the fused quartz rod as shown by Fig. 30 and a sample to be 
measured is placed between these by means of a liquid such as poly- 
isobutylene which has a high shear elasticity. If the specimen has a 
small attenuation, this can be measured by taking the difference in the 
amplitude of successive reflections. If the specimen has a high loss, this 
does not work and another method has been used which consists in 
sending a pulse from both erystals.’® One crystal is then used to receive 
and it receives the wave sent through the sample and the wave reflected 
from the fused quartz-sample interface. By adjusting the amplitude 
until these two are equal and the frequency or phase of one channel 
until the waves cancel, a ratio of amplitudes and a frequency of half 
wavelength are accurately determined. From these the velocity and 
attenuation can be calculated. 

This method has been applied to measuring the longitudinal and shear 
velocities of polyethylene and 6-6 nylon. The polyethylene was of “‘equi- 
librium’? crystallinity and average molecular weight corresponding to 

“HJ. MeSkimin, ‘Ultrasonic Measurement Techniques Applicable to Small 


Solid Specimens,’’ J. Acoust. Soc. Am., 22, No. 4, July 1950, pp. 413-418. 
*H. J. MeSkimin, J. Acoust. Soc. Am., 28, No. 4, pp. 429-435 
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an intrinsic viscosity in xylene of [n] = 0.89 at 85°C. Fig. 31 shows the 
longitudinal velocity of polyethylene plotted as a function of frequency 
and temperature. The velocity rises with frequency and a dispersion 
is indicated. This is confirmed by the attenuation per wavelength curve 
for two different frequencies plotted as a function of temperature, Fig. 
32. A definite dispersion is seen to occur with an activation energy of 
about 12 + 2 kilocalories per mole. This could occur in either the \ 
constant or the shearing constants yu, but the data of Figs. 33 and 34 
show definitely that it occurs in the shear component. Fig. 33 shows 
the shear velocity for four temperatures plotted as a function of fre- 
quency. This can be fitted for 30°C with a single relaxation mechanism 
having a relaxation frequency of 8 megacycles. To agree with the meas- 
ured attenuation and velocity, there has to be a spreading of the single 
relaxation over a range as also occurs in liquids. The indicated shear 
stiffness below this relaxation frequency is 2.6 < 10° dynes,/em*. Some 
1°] 
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Fig. 29—Spatial structure of nylon 6-6. 
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data on the zero frequency shear modulus is obtained from the Young’s 
modulus for a static pull which is from 30,000 to 50,000 pounds /square 
inch. Since the Young’s modulus is three times the shearing modulus, 
the zero frequency shearing modulus should not exceed 1.1 & 10° dynes 
em’. Hence one may expect that other relaxations will occur at lower 
frequencies, 

Fig. 34 shows the attenuation per wavelength for shear waves. The 
solid line for 30°C represents the calculated attenuation per wavelength 
for the model assumed. If all the lissipation were due to shear mech- 
anisms, the calculated attenuations would occur as shown by the 30°C 
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Fig. 30- Ultrasonic pulse method for measuring the velocities and attenua 
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Fig. 34-—Attenuation per wavelength for shear waves in polyethylene 
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points of Fig. 32. Most of the loss is accounted for by shear mechanisms, 
but it appears that some compressional mechanisms may also be present. 

The mechanism causing the relaxation in the megacvele range for 
polyethylene appears to be the same as for polyisobutylene, namely 
the relaxation of the shortest chain segment that is free to move. The 
chain segment acting appears to be longer than six chain units for similar 
measurements of nylon 6-6 show no relaxations in this frequency range. 
Fig. 35 shows the longitudinal velocity and Fig. 36 the attenuation per 
wavelength for longitudinal waves. Since the attenuation per wave- 
length is still increasing for nylon 6-6 at 25 megacycles a still shorter 
chain segment may be operating for this material. The shear velocity 
and attenuation per wavelength for nylon 6-6 are shown by Figs. 37 
and 38. 


Fig. 39 shows the shear stiffness of polyethylene and nylon 6-6 plotted 
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Fig. 37—-Velocity of shear waves in nylon 6-6 plotted as a function of temper 
ature and frequency. 
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Fig. 38—Attenuation per wavelength for shear waves in nylon 6-6 plotted as a 
function of temperature and frequency 
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against | 7 where 7 is the absolute temperature. Both are plotted for 
S me and 30 me. The dispersion in both materials is evident. Below 
30°C the shear elasticity of polyethylene yaries exponentially with the 
temperature with an activation energy of 2.72 kilocalories per mole 
Above this temperature a deviation occurs due to the approach to the 
melting temperature. Nylon has a smaller variation with temperature. 

Comparing the longitudinal and shear Wave measurements one can 
calculate the Lamé J elastic constant and this is shown plotted on Fig. 40 


for both polvethylene and nylon 6-6 as a function of temperature tor 
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Fig 11 Equivale nt shear and compression il viscosities for polyethylene and 
nvlon 6 6 plotted as t function of frequency for a temperature of 25°C. 


two frequencies. The dispersion of \ for polyethylene is small but is more 
prominent in nylon 6-6. This correlates with the larger compressional 
Viscosity Component present for nylon 6-6 which as shown from Fig. 41 
IS “as large us the shear \ ISCOSILYV According to the structural rearrange 
ment theory of compressional viscosity due to Debye, compressional 
viscosity can enter when some part of the chain can rearrange from one 
stable state to another stable state as a function of pressure. This re- 
arrangement occurs across a potential barrier and hence requires a finite 
amount of time to occur. This lag in the rearrangement results in a 
compressional VISCOSITY and as the frequency is increased, a lrequency 


is found for which the motion can no longer occur in the time of a single 


16 P. Debye, Z. Elektrochem., 45, 1939, p. 174 
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eycle and the A constant increases. It appears from these measurements 
that the dipole binding present in nylon 6-6 allows a greater structural 
rearrangement under pressure than can occur for polyethylene which 


has only linear chains. 


VI. CONCLUSIONS 


Measurements in dilute solutions, in pure polymer liquids and in non 
rigid solid polymers have all shown the presence of a shortest segment 
whose relaxation leads to a crystalline type of elasticity. In dilute 
polymer solutions the presence of a configurational type of relaxation 
and an entanglement relaxation of the shortest chain segment have 
been shown. For pure polymer liquids a quasi-configurational type of 
relaxation has been found for chain lengths greater than 60 segments, 
but for chain lengths less than 40 segments this type of relaxation dis- 
appears. From the difference between the high frequency shear elastici- 
ties measured for polyethylene and nylon 6-6 and the static measure- 
ment of Young’s modulus, it appears that there may be other relaxations 
in these materials for lower frequency ranges. 

For pure polvisobutylene and for nylon 6-6 there appear to be struc- 
tural changes induced by pressure which account for a compressional 
viscosity and a dispersion in the \ elastic constant. This effect is smaller 


for polvethy lene 


APPENDIX BEFECT OF LIQUIDS ON THE PROPAGATION OF SHEAR WAVES 


IN RODS 


For radially symmetric rods, the tangential particle displacement ug 


in the rod is given by 


Ug Ji(kr)e’ rr (1) 


(1A) 


All other displacements are zero. In this equation waves are considered 


to be travelling in the +2 direction with a propagation constant @ = 
A + jB, where A is the attenuation in nepers per em and B the phase 
shift in radians per centimeter. uw is the shear stiffness which may be 
complex to take account of the dissipation within the rod. 
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From the defining relations for the stress strain equation 


= 3 Ou Ug 
or 2 


and the tangential particle velocity dug dt, one can calculate the im- 


pedance Z per square em. of cylindrical surface at a. This relation is 


Z “Pre juk k (ka) 2 


lg w | Ji(ka ka 
J 


Since only the first mode is excited, parameters can be adjusted to 
keep k quite small, te. (ka < .2) and equation (2) can be simplified by 
using power series expansions for the Bessel functions. Neglecting higher 
order terms this results in 
y juak 

dw 

To evaluate the impedance of the liquid surrounding the rod, the tor 
sional wave is first propagated along the length of the rod without the 
liquid, i.e. with Z = 0. Then from equation (3) / 0 and from equa 
tion (1A) 


where Ap and Bo are respectively the attenuation and phase shift: in 
the rod alone. With the small loss in metal and glass rods A 
taken equal to zero and 


» 


b w/V p/p = w/v 


where vp is the velocity of propagation n the rod alone. 


When the liquid surrounds the rod, however, 


For the usual case where (B ; > U - Ag), equation (6) approxi- 


mates 


AB + jAA) = 2B) (—AB + jAA 


where AB is the increase in phase shift per centimeter and AA the in- 
crease in attenuation per ecm, both directly measurable quantities. The 
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then gi 


2BISA + jAB] = *" (aA 


Since AA and AB are the attenuation and phase shift changes per unit 
th. then if fis the length of the rod, covered the total attenuation 
ph ist shift changes will be AA and AB multiplied by 2{ Hence it 


AAy and AB) are the measured attenuation and phase changes, the 


le 
ey 


trie 
mpedance Z becomes 
+ JAB 


Vhis derivation neglects the change of phase occurring at the inter 
section between the rod having no liqnuid and the rod surrounded by 
the liquid, but it can be shown that this is small and moreover, the 
change in the wave on leaving is equal and opposite to that occurring 
on entering and hence this correction cancels out. However, if the 
liquid is viscous enough there is a correction due to the fact that the 
mensured impedance of equation (S) is tora eviindrical surface, whereas 
the desired impedance is the characteristic plane wave impedance. 
Obviously if the radius of curvature is sufficiently large no correction to 
equation (S) need be made. To obtain a suitable criterion one may con 

ider waves propagated into the liquid from the surface of the rod and 
solve for the impedance per square em. of the cylindrical surface. This 
neglects the variation with z, but since the wavelength along the rod is 
quite large, little error results from neglecting variations with 

An outgoing evlindrical wave in the medium may be represented by 


iY, (kr)\e"" ” (kr) (9 


i 
vhere the primes reter to the wave in the liquid and 
pw p w 
re Z 


where Z Ve p' is the pleat e wave impedance of the liquid. 


aring stress 


uv 
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and the tangential velocity may be obtained as before and the complex 


impedance over the cylindrical surface determined to be 


[ ¢ * (Ae (k a / 
yt | Oe ay = | / wH” (kr)’ (11 


or 


Noting that for plane waves Z \ pu’, one may eliminate uw’ from 
(11) and evaluate Z at? a with the result 


Hi? (ka)’ — (ka) 


The above equation can be used to obtain a solution for Z, in terms of 
the measured value of the cylindrical impedance Z for any set of para 
meters that may apply. Except for very heavy loading of the rod, how 
ever, the results of equation (8) may be used directly with little error. 
For example calculations indicate that for | ka | ’ 10, the multiplier 
of Z, of equation (12) is approximately 1.07 Z 7.4° for phase angle of 
(ka)’ of ~25°. For (ka)’ < 10, the correction multiplier rapidly becomes 
important. This same correction is applicable for the torsional crystal, 
but since this is only used for dynamic viscosities less than LO poises, a 


correction is seldom necessary. 








Relay Armature Rebound Analysis 
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Rebound of mechanical structures subse quent to impinging on slops gen- 
erally has deleterious effects on their pe rformance and should, the refore, be 
minimized. A considerable reduction in rebound can often be obtained by 
introducing additional degrees of freedom to the structure. 

{ mathematical treatise of the dynamics of rebound motion of systems 
representing idealized relay armatures is presente d. Normalized differential 
equations of motion and their solutions for the “free” and “impact” inter- 
vals are derived for systems hai ing one, two, and three de grees of fre edom. 
allowing the rebound behavior of a specific system to be calculated. The equa- 
tions of series of rebounds, and possible combinations of such series are con- 
sidered next for systems having one and two degrees of freedom. The field of 
possible rebound maxima is mapped for a practical range of mass distribu- 
tion constants, coe flicients of restitution, and force ratios. A sufficiently broad 
optimum design region 1s indicated 

The results of this analysis have been checked closely on a model and have 


led to appre ciable reduction of armature rebound in re lay de Signs. 


I INTRODUCTION 


In numerous types of mechanisms it is desirable to arrest the motion 
of a member at a particular point and to maintain it in this position. 
One of the simplest means of accomplishing this is to allow the moving 
member to impinge on a fixed member (stop) and to provide forces to 
tension it against this stop. Because the member to be arrested possesses 
kinetic ene rgy and because the stop cannot generally absorb all of this 
energy, the moving member will rebound from the stop. The rebound 
motion generally deteriorates the performance of the mechanism and 
should be minimized. 

Investigation of this phenomenon has been stimulated by the armature 
rebound problem in relay operation, where rebound from the front stop* 
tends to reclose contacts and must therefore be compensated for by 
additional (waste) travel, resulting in deleterious effects on speed and 

* Among relay designers the front stop has been generally referred to as “back- 


this paper the terms front stop and heel stop have been used through- 
sier identification 
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magnetic characteristics. Analysis in this paper will be directed towards 
relay armature systems, but it is also applicable to rebound in similar 
mechanisms. 


Il. STATEMENT OF PROBLEM 


Analysis will be restricted to planar motion of armature systems 
having one, two, and three degrees of freedom as depicted in Figs. 1, 
2. and 3. Generally one stop must be proy ided tor each degree of free 
dom, although in the three degree ot-freedom system of Fig. 3, two of 
the stops have been combined. 

Applied forces Ff, , Ff: , F3 , have been chosen so as to be most easily 
correlated with actual relay designs. 

The initial condition in all cases will be a pure rotation about the 
‘‘heel” just prior to a “zero” impact at the “front”? of the armature. 
The “zero” impact will be followed by rebound motion and impacts at 
the various stops eventually bringing the armature to rest. The object 
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will be to minimize rebound motion at the front, since this is usually 


near the point actuating the relay contacts. 

The basic problem is then to find the response of the armature subject 
to aperiodic but well defined impulses, which are functions of the 
positions and velocities of the system. 


Il ASSUMPTIONS 


In order to facilitate the solution of this problem, the following 
modifying assumptions are made: 

1) As mentioned in the previous section, analysis is restricted to 
planar motion, 


2) The armature is assumed to be a rigid body. 


3) Stops are assumed to be very stiff, massless springs capable of 
energy absorption during impact with the armature. The associated 
coefficient of restitution is assumed constant. Core and stop vibration 
are neglected. 

(4) The tensioning forces FP’; , Fs , 3 are assumed to be constant forces. 

This is fairly closely true for moderate rebound amplitudes of practical 
relay structures. ) 

(5) All displacements are small relative to the dimensions of the 
system and in particular the angular displacement @ is sufficiently small 
so that 


cos @ 


sin 6 


IV. DERIVATION OF EQUATIONS OF MOTION 


The derivation of the equations of motion resolves itself into the 
solution of two different types of intervals: 

|) Free Interval: This is the period during which the armature is 
not in contact with any of its stops and only the tensioning forces are 
acting 

2) Impact Interval: During such intervals the armature Is in contact 
With at least one of the stops. The stiffness of the latter is assumed so 
high that the tensioning forces during this interval may be neglected. 

The three-degree-of-freedom case will be considered first and the 
others subsequently deduced from it by allowing some of the constants 
to approach Zero. 
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A ° Free Inte ri al 


The motion of the armature will be described by the displacement at 
the stop points: 2; , 23, «;.* Let m be the mass and R the radius of 


gvration of the armature about the center of gravity. The latter is 


located by the dimensions (AR, &R, and GR relative to the stop points, 


i.e., the points on the armature which contact the stops in the rest 
position (Fig. 3). 

The equations of motion are derived in \ppendix I and are put into 
dimensionless form: 





where: 


{ + 13 Cn it 
(5 + sr : (1 
{5 ) Cr ff 
Yio, Yo, Yso, are the initial velocities and Ho, Yo, Yso the initial dis 
placements for the free interval in question. 
The equations ol motion for a two-degree-of-freedom system are 


obtained, it F's 0. Then for the two coordinates of interest: 


s [eu + om (F 
[cn + cn (B)] ( 


\ summary of all notations used in this paper is given in Appendix I\ 





176 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1952 


a ; | [Pe 
For a one-degree-of-freedom system jj2 = Cyn + Cx (F = 0, whence 
a 


: as : i 2 + vw} (5 


B. Impact Interval 


The change of velocity at point “7”? due to an impact at “7” is, by 


definition of the coefficient of restitution ‘hk’ 
Ay = — (14+ &k;)a, 


It is assumed here that the action of the stops are true impacts, i.e., 
the changes in velocity take place while there is negligible motion of the 
body. The velocity changes then occur as instantaneous rotation about 
the conjugate axis, leading to the general relation for an impact at 


c 
point er is 


Y jon Y jetn nT K 5M iecn 1) (6) 


The first subscript indicates the coordinate, the second subscript 
indicates the beginning (0) or the end (e) of the free interval described 
by the third subscript. The impact transfer coefficient K,, relating ¢ 
velocity change at point ‘‘7” to an impact at point “7’’: 

( 
Kx = - “(14+ ky) 
Ci 

Iquations (1) through (7) allow any one specific case to be mapped, 

if the mass distribution and force ratio are known. A sample of such 


mapping of rebound motion for a rectangular two-degree-of-freedom 
armature appears in Fig. 4. 


ANALYSIS OF REBOUND PATTERN-—-ONE-DEGREE-OF-FREEDOM SYSTEM 


The rebound pattern for the one-degree-of-freedom system—as de- 


rived in Appendix I[—-consists of an infinite series of parabolic ares of 
diminishing amplitudes. The structure comes to rest after a finite time 
interval. The maximum rebound occurs during the first bounce and 
equals 


(8) 
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where 


The system returns 


VI. ANALYSIS OF REBOUND PATTERN rWO-DEGREE-OF-FREEDOM SYSTEM 


The reason for choosing a two-degree-of-freedom system over a one 


degree-of-freedom system would be, in keeping with the philosophy of 














Front and heel motion of plate type armature. 
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this treatment, to reduce },, the greatest excursion at the front. In 
order to simplity mapping, this maximum excursion will be expressed 
is y,, the ratio of ¥, to) 

1. Thus 2C} 
degree-ol 


ol the 
Wi tl 


as given by Equation (8) for the case of 
is the ratio of the greatest excursion of the two 
{-freedom system under consideration to the greatest excursion 
corresponding pertectly elastic one-degree-of freedom system. 
constants which are functions of the mass 
distribution relative to the stop locations: 


| Introduce two Deasic 


\/ (11) 
Mh constant coupling coefhicient \s 
VW W,,, the two-degree-of-freedom system under consideration here 
has only one such non trivial constant VW, 


represents a mechanical 


Phe second constant represents a force transformation factor from the 
de coordinate to the ** ‘ coordinate: 


pP 


(12 


ln the analysis of the two-degree-of-freedom system only Po is important. 


If there is to be any heel motion, the zero” impact at the front must 
IMmpart a positive velocity to the heel. By Equations (6), (7), and (12), 
this requires that Py» be negative, which in turn implies that hh > 1. 
For the limiting case of (6 Ley M, 0 and no coupling exists 
between the heel and the front. Physically this means that the two 


tops are the centers of percussion of each other and the system will 
act as a simple hinge 


ic 


With the above foundation, it is possible to analyze the patterns of 


motion and maximum rebound amplitudes 


[immediat li Following ee Ue [mpact 
ro” impact at the front, both front and heel will litt: off 
with impact Equation (6) and continue to move in ac- 
vith the tree interval Mquations 


hk ACE ordance 
eordat ('¢ 


OccUrs il 


1). Whether the next impact 
the front or the heel depends on their respective periods, 
ft, al « 4 


Py 
M, 
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where: 


A large value of t/t will result in a series of heel impacts and the 


heel will come to rest while the front is still displaced from the stop 
This will be ealled a complete heel series. A small value of ty ts results in 
a similar complete front series. If t,t is near unity, a limited number of 
impacts on one end are followed by an Impact on the other end, ete 
An analysis of front and heel series follows: 


B. Front Series 


If t/t < 1 a series of front impacts occurs. The impact velocities at 
the front are 


Yi 1, &, &, i 14 
The corresponding time intervals are 


9]. 9] : 


where 
A 


During this time, the heel velocity and displacement are given by 


| 28 


+ 


where 


The velocity and displacement at the heel after a given number of front 
impacts are obtained by a summation of Equations (16). For 
plete front series n —» =~, and 


a Com 
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In addition it is useful to set down energy equations in order to 
simplify evaluation of greatest rebound for the various groups of re- 
bound patterns. The kinetic energy function T is evaluated in Appendix 
I. A potential energy term V—the work done against Ff; and F; from the 
equilibrium position—is introduced. If 7 is the total energy of the 
system prior to the ‘‘zero”’ impact, then 


i - =e Ue = Yye ‘ 
19 Py» (18) 


— 20 (yy 4. fy ) 


2 My ; 2 , 2M 


The energy loss due to n front impacts is 


a(? uM =) (i = ki”) (1 — My)ijie 


For a complete front series n — DQ, and 


a(Zt%) = 0 — atwdite (20) 


If a complete front series follows the “zero” impact, jio = | and 


a(=3 a) (1 — My) (21) 


\fter completion of this ‘initial’ front series, the system maintains 
only one degree of freedom (rotation about the front) until a heel 
impact occurs. By setting yj TF Yo 0 in (21) we obtain the heel 
impact approach velocity 7 Pas. 
Apparently energy loss due to n front impacts is a function of My , 
, and the approach velocity of the first Impact. 


Heel Series 


\n analysis similar to the above can be made for partial and com- 
plete heel series following the “zero” impact. This is demonstrated in 
Appendix IIT, vielding, for hy, k.* 

Dy . 2 


1 


om + j:)? on — k) ) 
APy(1 APy _kKL—k) _ yyy] | 
Ba } B 1+ / 


ky} 2AP, kK 9) | 
. ~ My(1 + k) |} 
‘| B 1 + k) | 


re general form k,; # ky ean be obtained as indicated in Appendix III 
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The energy relationships for heel series are 


] my My» = M2) ii 
Pr. 


For a complete series n — ~*, and 


Pr 7 r) Mit ~ Mi 
1 iP P?, {260 


If a complete heel series follows the ‘‘ zero” Impact, te, 


and 


Finally, for the special case where a complete heel series follows an 
initial complete front series je Py», and 


a(* = i ) Mi(M — 1) (26) 
To 
It is to be noted that the energy loss due to a partial heel series is a 
function of My, Pi, , hk, , and the approach velocity of the first impact, 
but that the equation for a complete heel series does not contain kz . 
Finally, a complete initial heel series is a function of only My and ky, . 


D. Complete Mapping of Problem 


Equations (1) through (26) make it possible to completely map the 
two-degree-of-freedom rebound problem. The relative maximum ampli- 
tude 2CY, and the rebound pattern will be determined. 

Examination of the necessary equations, show that 2CY; is in all 
cases a function of four parameters: hk; , kz, My and Py»f. Of these, ky 
enters only if a partial heel series occurs prior to the time of maximum 
rebound. If it is assumed that for this limited group of cases hk. = k, = k, 
the number of parameters is reduced to three: k, My, Pof. 

In Figs. 5 to 10, 2C Y, is plotted against Py»f for the most useful range 
of 1/8 < My < 1/2.5,0.3 < k < 0.6 and O « Prof < 30. 

As Pyf is increased from zero to infinity (corresponding to an increase 
in the heel tension /’2), the rebound pattern goes through some or all of 


five regions. The criterion for location in any one region is based upon 


the parameter 
Pr 

, f ( - -) 
0 Mi: a oe 
7 l -+ Pywf b k 


1+ 





182 fHE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1992 


Region I—-Complete initial front series for 1 <Q < 1/k. Within this 
l M.)k ; 
region, if My. > the maximum rebound occurs during the 
lL + Prof 
first bounce and 
(1 — Mis)k° : 
~ (28) 


2CY 
1 + Puf 


maximum rebound occurs later, it must occur during a com- 


If the 
plete heel series which follows the initial complete front series. From 


equations 21) and (26 
2CY, = M’ (29) 


1 


By comparing Kquations (28) and (29), the critical requirement 
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Mi, 


while the first rebound maximum, shown in solid lines on Figs. 5 to 10, 


for the latter case is that Pyf > k° —1. It should be noted that 


is always realized, the later rebound given bv (29) is an uppel limit 
shown in dashed lines—and is not always realized. In the dashed re 
gions, phasing is extremely critical and small variations in the param 
eters may cause large variations in maximum rebound. From = an 
engineering standpoint these regions are essentially undesirable 


Region I]1—Partial initial front series for 


h 


This region is one of critical phasing, and attention is limited to 


special cases leading to maximum rebound. These cases occur when a 
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Fig. 6—Relative maximum rebound amplitude for M, 
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heel impact immediately follows the last front impact of the series. 


These cases occur at 


(30) 


and lead to rebound amplitudes 


ICY, = Mi + | My)(k’" — My(1 — k")’] (31) 


In kigs. 5 to 10, these special points are plotted and connected with 
straight dotted lines, which therefore indicate upper limits to rebound. 
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Xegion I11——Partial initial heel series for 


ae l+k 
k © ~kO — k) + Muk(l + k 


This is a region of critical phasing, and values were determined only 
for the maximum cases, where a front impact just precedes the last 
impact of the partial initial heel series. Here: 





Cc 
> 


9 
UO 
N 
. 
Ww 
a) 
a 
is 
ee 
Q 
= 
aq 
‘@) 
Zz 
> 


O 


vy 

1 \ 

; : 
y it 


MAXIMUM REB 


0.4 a 


RELATIVE 
O ° 
oO 


ee | 
*s Lae] \ 
4 \ 
aA \ 
of. + \ + es 
XK, 











Fig. 8—Relative maximum rebound amplitude for M). 
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and: 
2CY; ( VQ 
(1 

k" ere 

Region I\ 
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+ k)(l — k")P) 


Complete initial heel series. A complete initial heel series 


implies that when the heel has come to rest, the front is still away from 


the stop. When the front finally meets its stop, the situation is identical 


with that just prior to the zero impact except that the energy content of 


the system is lower. The pattern must then repeat with diminished 


amplitudes. For this region we recognize two groups of cases. The first 
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group is that in which the front velocity is positive at the completion of 
the heel series. In that case 


—_—__—___— >Q>=- 
Myi(l + k)? é k(1 


and 


(] My)k 
l + Prof 


d 


2c yY; 


For the second group the front velocity is still negative when the heel 
comes to rest from which point on the system acts as a one-degree-of- 


0.30, 





K=0.6 
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Fig. 10-—Relative maximum rebound amplitude for My 
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freedom system until the next front impact. The requirement for this 


group is that 


+ k) 
T Mp1 
and the maximum rebound is given by 
2CY, = My —(1 My»)[((1 — k°) + Mu(1 + k)’ (35) 


It is to be noted that in the upper part of Group 1 the amplitude in- 
creases with successive heel impacts. This can be explored through the 
use of Equation (22). For simplicity of mapping, however, the limit 
given by Equation (35) has been extended back from the lower boundary 
of Group 2 until it intersects the line marking the first rebound ampli- 
tude of Group 1. 

In Figs. 5 to 10 the respective regions have been identified by means 
of the symbols indicated below: 


Region | from Ppyf = 0 to 

Region I] from to A 

Region II] from A to x 

Region IV, Group 1 from to O 

Region IV, Group 2 from to Ppf —- « 


kK. Discussion of Rebound Charts 


Aside from quantitative data contained in Figs. 5 to 10, the following 
general trends are of interest: 

For values of My, > 3, and the values of k under consideration, most 
of the useful range of Pyf involves critical phasing and the rebound 
maxima are relatively high. 

Kor values of } < My < }, consistently controllable rebound ampli- 
tude may be obtained. 

kor values of My < } rebound increases again and the structure 


approaches the one-degree-of-freedom case. 


VII. ANALYSIS OF REBOUND PATTERNS THREE-DEGREES-OF-FREEDOM 


SYSTEM 


Rebound pattern analysis as in Parts V and VI has so far not_been 
performed for the three-degree-of-freedom system, partly because of 
complexity, and partly because for the system of Fig. 3 friction at the 


hinging stop will greatly influence the motion. 
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However, it is felt that the approach and notation of the analysis 
presented here is sufficiently general to allow extension of the rebound 
pattern analysis to the three-degree-of-freedom case. At any rate, with 
the assumption of the magnitudes of frictional forces, the basic equations 


of Part IV may be used to plot any particular case 


VIII. ARMATURE REBOUND MODEL 


In order to verify the formal analysis presented in Parts III, [IV and 
V, a large model of a two-degree-of-freedom system was constructed, 
It consisted essentially of a large bar constrained to move in a plane, 
biased against two stops, and to the ends of which writing pens were 
attached. As rebound conditions were simulated by releasing the bar 
against its stops, chart paper moved at right angles to the bar motion 
and thus produced a record of end displacement versus time. 

By changing spring members and attaching masses to the bar, it was 
possible to vary the mass distribution and the biasing forces. 

The results obtained closely agreed with those suggested by the 
analysis. The maximum rebound amplitudes were generally somewhat 
lower probably due to frictional effects. 


IX. RELAY DESIGN CRITERIA RESULTING FROM ARMATURE REBOUND 
ANALYSIS 


A. Limitation of Analysis 


The assumptions which this analysis is subject to have been described 
under Part IT. As applied to relays and probably the majority of mechani 
cal structures, one assumption is most frequently and severely violated. 
The stops, which have been assumed to be stiff springs associated with a 
definite coefficient of restitution are, in practice, massive bodies which 
dissipate energy through excitation of high frequeney modes of vibration. 
Accordingly, the assumption that the stops are at rest is violated, 
particularly if the mechanism is subject to repetitive (pulsing) im- 
pacts and the stop vibration does not decay greatly in the repetition 
period. 


However, mechanisms designed in accordance with this analysis have 


performed well even under moderate pulsing conditions if the sensitive 


phasing region was avoided. In addition, every effort should be made to 
reduce the amount and duration of stop and mounting structure vibra 
tion by making them stiff, massive, and dissipative, if possible 
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3. Design Crite ria 


1. Type of Armature Structure. 

The selection of the number of degrees-of-freedom for an armature 
structure depends on the expected coefficient of rebound as well as 
practical considerations. 

It can be shown without great difficulty that for very low coefficients 
of rebound the one-degree-of-freedom system is preferable. This is quite 
obvious when one considers the limiting value of k = 0. In this case the 
one-degree-of-freedom system will have no rebound whatsoever, while 
the two-degree-of-freedom system has a heel bounce followed by re- 
bound at the front. The value of k below which the one degree system is 
preferable varies with the mass distribution relative to the stop points, 
being O.18 for a rectangular plate armature with stops located at its 


ed 


xperience indicates that / in most practical relays and similar 


| ‘ 


mechanical structures varies from 0.3 to 0.6. Hence the two-degree-of- 
freedom system is superior in all practical cases to the solidly hinged 
armature, 

\s far as three and higher degree-of-freedom systems are concerned, 
it may be said that generally the greater the number of modes resulting 
in impacts, the quicker the rebound energy can be diverted and dis- 
sipated and the lower theoretical rebound values can be obtained. This 
consideration would favor systems containing many degrees of freedom. 
However, while multi-degree-of-freedom systems can reach very low 
rebound values, their motion (phasing) must be very closely controlled 
or they may prove to be inferior to simpler systems particularly under 
vibratory (pulsing) operation. It is this difficulty which makes it appear 
that the two-degree-of-freedom system offers the best promise with the 
three-degree system also quite promising. By the same reasoning, 
additional spurious rocking modes should be minimized. 

2. Armature Mass. 

The armature mass should be as low as possible. This will minimize 
stop and structure vibration. In addition, in relay applications hght 
armatures tend to increase magnetic drag’ losses of the armature 
during the release motion 

3. Stops and Mounting Structure 

\s discussed before, it is desirable to reduce the amount and duration 
of stop and mounting structure vibration. 

t. Coefficient of Restitution. 

The coefficient of restitution should be kept low. Stops having low 


stiffness should, therefore, be avoided. 
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5. Biasing Forces. 

F should be kept as high as practicable 

For proper energy loss during impacts, the motion between impacts 
must occur outside the region of the compression, 1.e., the armature 
and stop must separate. Therefore, because all practical stops have a 
finite stiffness, the biasing forces (fF; , Fs, ete.) should produce a static 
deflection less than say, arbitrarily, 5 per cent of the maximum expected 
rebound amplitude. 

6. Design Parameters for Two-Degree-of-Freedom Systems. 

As clearly indicated in Figs. 5 to 10 for the practical range of coef 
ficients of restitution, most consistently good results are obtained with a 
coupling factor VM, rs to 4. This factor is most easily adjusted by 
correct placement of the front Stop. 

lor the above range of Jy the force ratio Fs F, should be such as to 
make the product 


l 
6 


(Note: For a rectangular armature structure with the stops placed at 


its edges My. 1, Pe ; and force ratios in the neighborhood of 8 are 
desirable. ) 
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APPENDIX | 


DERIVATION OF BASIC EQUATIONS OF MOTION THREE-DEGREE-OF-FREEDOM 
SYSTEM 


\ 1) Free Interval 


The motion of the armature will be described by the displacement at 
the stop points, X11, %, 23. Let m be the mass and FP the radius of gvra- 
tion of the armature about the center of gravity. The latter is located 
by the dimensions 4,2, &R, and GR relative to the stop points (Fig. 3). 
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The rotation and displacement of the center of gravity is then 


Ty = 


\pplying LaGrange’s Equation to the above, the 
are obtained: 


d (} ” oT os 
dt \0q, aq, 


The quations fs may be solved for Uy, U2, 1&3 


grated, vielding 


(CuF: + CrP, + CisFslt° 
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where 

hl; 

(1 — 4h) 

— bots 

iio , Loo , Xo are the initial velocities, x4 , Xe , Zo the initial displace- 

ments for the free interval in question. Interpretation of the analytic 
results is simplified by the introduction of normalization. Let +, be # 
just before the “zero”’ impact and define 
hoe F, 


— . ee 
Lef rm 


Igm 
F, 

Dividing Equations (d) by war yields the normalized equations 
motion: 


(FP) " y (F \ 
12 (F,) 

(Fy) 
22 (F, 


(F.) 
2 F ) 


0 


‘ai 








(2) Impact Interval 


The change of velocity at point ‘‘7”’ due to an impact 
definition of the coefficient of restitution ‘‘k” 


Ai, = -(l+ kz 


, 

Since this velocity change occurs as rotation about the conjugate 
point as an instant center of rotation, the impact relationships may 
be written, for an impact at point “1” 


Aq =. oe (] + ky) ay 


= —(] + kya &2 
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Cis 
Ch 


(1 + ki)ay 


Similarly it can be shown that impacts at points (2) and (3) follow 
the same pattern. The general impact relations for impact at point ‘‘2”’ 
are then 

Y jon = Y jen 1) t K jAiew 1) (6) 


‘The first subseript indicates the coordinate, and the second subscript 
indicates the beginning (O) or end (e) of the free interval denoted by the 
third subscript. 


The impact transfer coefficient Aj; relating a velocity change at 
point “‘7’”’ to an impact at point ‘‘7’’: 

’ C " 

kum — = + ks) (7) 


i 
APPENDIX II 
\NALYSIS OF REBOUND PATTERNS—ONE-DEGREE-OF-FREEDOM SYSTEM 
‘The equation of motion of this system is 
Yr sce" + Hront) + 


where 


and is measured from the start of the particular interval of free motion 
in question. The impact relationship is 


Yion * Kyirecn 1) 


The motion consists of a series of parabolic ares having periods of 
Y ‘ Y ¢ Y op2 y ‘ —lyyy nm : 
C in general, or 2/C, 2k/C, 2k°/C, ---, 2k” /C. The time elapsed 
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is a convergent series and approaches, for a complete series: 
) 
Lim — {1 4 
( 


nen 
The maximum rebound amplitude in any interval is 
The maximum excursion occurs during the first bounce at 
equals —k°/2C. 
AppeENDIX III 
ANALYSIS OF REBOUND PATTERNS-~TWO-DEGREE-OF-FREEDOM SYSTEM 
The equations of motion of this system are 


y2 . / r 
yn = sAt a Yronl tT Yon 


wy = Bt" UO Yoont’ t Y20n 


’ 


where A = Cy + 


B 


F, 
t measured from the start of the par- 
7 ticular free interval in question. 

A. Complete Front Series 


At the beginning of a front series 


In a manner analogous to that for the one-degree-of-freedom system 


each front impact reduces 7, to —kyj, . Therefore, after the n*® impact, 


Yion 
and the time elapsed in the n*® interval is 
2k? 


re . 1 ( 
A 
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At the heel, from (g), the heel velocity preceding the n‘* impact is 
Y2e(n—1 


Y20(n 1 = a Bt (k) 

The velocity change during the (n 1) interval is then equal to 

BT, . From Equations (6), (7) and (12), the change in velocity during 
the n*® impact is —Pype(l + ki)kf jie 

The total change of y between impacts is then 

ryy ‘ rl. 
Y20n 1 yee Py(1 + ky)ky Yieod 
Similarly in preceding intervals 


Y20(n—1) — Y20(n—2 BT, Py2(1 t ky )ky 


Yied 


Y2 - Yor = BT, 


Y2a0 


Pll + ki) ico 
By addition of the above 


‘ Py r ky) kuti 


Y201 


n—l n—1 

Y20n — Yr20 ans B>- aw 2 Pil + ky) > kT ja 
m=} m=) 
2f 


3 n 1 
ale Yieo 


kt — Pr 
m= 1 


n-1l 
(1 + ki)jio 2, k” 


m=) 
The summations may be evaluated, yielding 


2B k, — kr 1— ky |. 
Yr20n Yoo = | i ; = — Py(l + hy) ly (1) 
d v1 —T 1 


Y20n Y 2 


To evaluate the displacements at the heel, Equation (g) yields 


n--1 Y 201 


nT na + 4BT1 
Adding these expressions for intervals 0 to n; the total change in y is 
1 


n nl 
3 — . ve m2 
pam on = - Yoom mT 3 B , ie 
m=] m=1 


2h (1 -- kz") 
oo 1 7 
Atl — ki) Yidl/20 
4 [2B = ake + Ha”) 
A*(1 — k,)? 
; 2Pph — kt 


— kp? + ki*”) 2 
A(1 — ky)? et 
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Expressions for an initial series may be obtained by setting 7,4 I 


Ya = Yro = O, and, finally, for an initial complete series m — = 


hence k” — 0, and Equations (1) and (m) become 


ss 2hy Bk, > 
~ ACL + k)? | , 2| 
eo = _ > 10 + k ) 


B. Complete Heel Series 


and 


For heel series, Equations (1) and (m) may be used by interchanging 
the initial velocities, accelerations, and impact transfer coefficients for 
those relating to heel motion: 


| | 2A ky — ke 
Yn —~ he * | T_-k 


Qke(l — kt 
Y- & 


Yon - /a= - 


2A(k2 — 2kp*? + hk") — 2Myke(l — ke 


Br :2) - BP» 
An initial heel series occurs when the heel strikes first after the “zero” 
impact. The first heel impact then occurs 7; 2P»/BA 


+ fy) after 
the zero impact and the initial conditions are 


Py(l + hy) 


Substitution of the above into (n) yields 


+ ky | 2A Py , 
Yin = —hy + is | AP a (y v2 -My2(1 + k, 


l = ky B 


The corresponding expression for jo, is quite involved. For the special 
case of k = ky = ke 
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AP,.(1 + k)’ asia AY 
B _B l—k 


If the initial 


(’. Partial Front Series 


The worst rebound occurs when heel and front impacts occur nearly 
simultaneously, with the front hitting first. From Equation (m) for an 
initial front series, this requires that 

B l k’ 
Q= 30) 
AP;, I: fen 

\fter » front impacts conditions are given by Equations (14) and 

(19) with y, y 0, and 


Myjz 2M pk" ye 
»2 Po 


12 


My) (1 . — 


This may be solved for #2 = Pwd k"). The maximum front excur- 
sion now possible-is that for a complete series of heel impacts. The 
above value of a in Equation (24) vields 


2CY Vi, + (J Myo)(k°’ ’ M,(1 —_ k")*) (31) 


D. Partial Heel Series 


The worst rebound occurs again when heel and front impacts occur 
nearly simultaneously, with the front hitting first. From Equation (9) 
for an initial heel series, this requires that 


B 
AP 





RELAY ARMATURE REBOUND ANALYSIS 
After n heel impacts 7% Py(l + k)k" and from Equations (19 
and (23) 
T+V 
a hl 
Ty 


This may be solved for ‘ A’) and after the 


front impact immediately following: 
ie = Pol + DA" — Poll + lk — Mul + (0 — | 


The maximum front excursion now possible is that for a complete 
series of heel impacts. The above value for gj in Equation (24) yields 


2CY; = 1 tl + I] is a 
My(1 


APPENDIX LV 
SUMMARY OF NOTATION 
A= Cn + Cr - 
b 


front tensioning force 
heel tensioning force 
coefficient of restitution at vertical front stop 


coefficient of restitution at vertical heel stop 
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vertical front stop location relative to Cw. 
vertical heel stop location relative to e.g. 
horizontal heel stop location relative to ¢.g. 


mass Ol armature 


Pu 
Myf ly l + Ay 
1+ Pi bo hy 


radius of gyration of armature about center of gravity 


Lam 


F, 


time 


basic period of front after ‘zero’? impact 
basic period of heel after “zero” impact 
duration of nt” free interval 

vertical front displacement 

vertical heel displacement 

horizontal heel displacement 


front velocity just prior to “zero” impact 


yreatest excursion (rebound) of front 
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A New Telephone Carrier System for Medium-Haul Circuits. R. 8. 
CaruTHers ', H. R. Huntiey’?, W. E. Kant, and L. Prperson'. 
Elec. Eng., 70, pp. 692-693, Aug. 1951. 

Low terminal costs and single-cable operation make this the first economically 
practical carrier system for medium-haul telephone circuits. Performance is not 


sacrificed for economy. 


A .15-Kw 500-M eqac ycle Grounded-Grid Tridode.* C. Ei. Fay', D. A. 
A. Haze’, and R. J. Kirener'. Proc. J.R.E., 39, pp. 800-803, July, 
L951. 


Short-Cut Method Aids Figuring Exhaust and Collecting Systems. W. 
H. FoGuie’. Heating, Piping and Air Cond., 23, pp. 75-78, July, 1951. 

A simple, workable method of determining static losses in industrial exhaust 
and collecting systems is explained here by means of a sample problem and its 
step-by-step solution. The method might be described as an “averaging out” 
process, whereby all duct sizes and lengths are averaged out with the efm and 
velocity to give a total linear footage of the average duct size and an average 


velocity. 


Areing at Electrical Contacts on Closure. Part I. Dependence Upon 
Surface Conditions and Circuit Parameters.* L. H. Germer'. Jl. Appl. 
Phys. 22, pp. 955-964, July, 1951 

In a low-voltage circuit the occurrence of an are between approaching elec- 
trodes is dependent upon the nature of the surfaces and upon the circuit induct- 
ance. For carbon surfaces, or noble metal surfaces which have been “activated” 
by operation in various organic vapors resulting in a carbonaceous layer, the 
limiting circuit inductance is somewhat above 10°*h, which is much higher 
than the limiting inductance for clean noble metal surfaces. This activation by 
organic vapors occurs for noble metals only and for certain vapors; for example, 
benzene derivatives. In the case of silver and benzene \ apor, it has been shown 
that the activation is due to adsorption of benzene onto a greasy surface layer 
and its decomposition there by the heat of subsequent closures. A metal surtace, 
which has been activated by organic Vapor, remains active indefinitely if there 

*A reprint of this article may be obtained on request. 

' Bell Tel. Labs 

A. 2. a 5 GO 

3W. E. Co. 
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ng at the surfaces; 


but with continued operation and accompanying 
arcing, the activating material is burned away, and the surface returns to the 
inactive condition if no activating vapor is supplied. 

Are voltages, which are independent of current and of ambient gas, as far as 
tested, have been measured for a number of metals and for carbon; the are 
voltage for carbon is quite erratic in the range between 20 and 30 volts, but for 
each of a number of metals the are voltage is steady. 

Arcing at noble metal surfaces, similar to that induced by carbonaceous 
material from organic vapors, can be produced also by insulating particles or 
insulating films. The active condition gradually disappears with continued are- 
ing, unless there is a steady supply of insulating material to the surface. 

The minimum are current has been measured to be 0.6 amp for active silver 
ind for carbon, and 0.03 amp for inactive silver. These are the currents at 


which an established are is extinguished. 


Tron-Silicon Alloys Heat Treated in a Magnetic Field.* M. Gorrrz'. 
Py lppl Phys... 2. pp QO4- 905, July, LOL, 

Heat treatment in a magnetic field has been found effective for iron-silicon 
illoys between two per cent and ten per cent silicon, the highest maximum per- 
meability being obtained at about 6.5 per cent silicon. In a single crystal of this 
composition, magnetized parallel to a (100) direction, the hysteresis loop is 
squared by the magnetic anneal and the maximum permeability is increased 
from 50,000 to 3,800,000, the highest value yet reported. 


Domain Boundary Motion in Ferroelectric Crystals and the Di-electric 
Constant at [igh Frequen y. C. Wrrrer Letter to the Editor. Phys. 
Rev.. 83, p. fos, July 15. 1951. 


1 Wethod for Determining the Propagation Constants of Plastics at 
Ultrason Frequencu s.* H. J. McSximin!. J. Acoust. Soc. Am., 23, pp. 
120-434, July, 1951. 

A pulse technique particularly suited to dissipative materials is described for 
measuring attenuation and phase-shift constants of plastics, using either trans- 
verse or longitudinal waves in the Irequency range of 5-50 me. 

A thin wafer of the material under test is placed between two identical fused 
silica buffers; and waves generated by quartz crystals at the ends of the assem- 
bly are transmitted simultaneously through the specimen in both directions. 
Comparison of transmitted and reflected components by means of a special 


balancing circuit provides information from which the complex propagation 


constant can be calculated, and hence dynamic rigidities and viscosities. 


Illustrative data for polvethylene and Nylon are given. 


of this article may be obtained on request. 
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Some Antecedents of Quality Control. E. C. Moninat. Ind Quality 
Control, 8, pp. 1O-L1, July, 1951. 


Traveling-Wave Amplifier Measurements. F. E. Rapewirre'. Electronics, 
24, pp. 110-111, Aug., 1951. 

tapid sweep-frequency technique used at 4,000 me can be applied to all 
broad-band amplifier measurements. Oscilloscope display shows transmission 
| 


accurate to 0.1 db and return-loss values up to 40 db 


Kirchhoff’ s Formula. its | eclor Analogue, and Other Field Kquivale nee 
Theorems.* S. A. SCHELKUNOFF!. Communications on Pure and A pplied 


Math., 4, pp. 438-59, June, 1951. 


Remarks Concerning Wave Propagation in Stratified Media.* S. A. 
SCHELKUNOFF'. Communications on Pure and Lpplied Vath., 4, pp. 
117-128, June, 1951 


An Achromatic Doublet of Silicon and Germanium.* R. G. Treutine!, 
Se Opt. Soc. Am., 41, pp. 454-456, July, LOS]. 

The semi-metals germanium and silicon have high transparency and high 
refractive indices over a wide range of infrared wavelengths and are stable to 
normal atmospheres, Their relative indices and dispersions make achromatic 
combinations possible ; and designs are given for axially corrected doublets of 
relative apertures f:2 and f:1. The optical homogeneity of the materials is dis- 
cussed: compositional variations are not considered an optical hazard, but there 
is evidence of structural imperfections in some specimens whose effect on optical 


properties remains to be evaluated. 


On the Motion of Gaseous Tons in a Strong Klectric Field. I.* G. H., 
WaNNIER!. Phys. Rev., 83, pp. 281-289, July 15, 1951 

This paper applies the Boltzmann method of gaseous kinetics to the prob- 
lem of positive ions moving through a gas under the influence of a static, uni- 
form electric field. The ion density is assumed to be vanishingly low, but the 
field is taken to be strong; that is, the energy which it imparts to the ions is 
not assumed negligible in comparison to thermal energy. Attention is focused 
upon the computation of velocity averages, and the drift velocity in particular, 
rather than a complete knowledge of the entire velocity distribution. It is shown 
in Sections C and E that the problem so formulated is completely soluble if 
the mean free time between collisions of ions and molecules is a constant: this 
is the case for the so-called polarization force between ions and molecules which 


predominates over other forces at low temperature. A method for obtaining 


averages to any desired accuracy in the general case is developed in Section D. 


The method is a plied to the hard sphere n ‘| for the high field range and 
DP} | 


~ 


cA reprint of this article may be obtaine d on request. 
1 Bell Tel. Labs. 
4 Bell Tel. Labs., Retired. 
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mass ratio 1. An application of the resulting formula (43) to experimental mate- 


rial has been published earlier. 


Evidence for the Noneubie High Te mperature Phase of BaTiO; a 
Woop!. Letter to the Editor. J. Chem. Phys., 19, p. 976, July, LOSI. 


Seven-Leaque Oscillator.* FB. ANpDERSON!'. Proc. T.R.E., 39, pp. SS1— 
890, Aug., 1951. 

A bridge-type RC oscillator is described which is continuously adjustable 
over a frequency range of 20 cps to 8 me in one sweep of a two-gang linear 
potentiometer control. Tracking requirements of the two-gang control are not 
severe. The output is available in four phases, and the frequency is an approxi- 
mately logarithmic function of the linear potentiometer setting. Practical limits 
of the frequency range are tentatively 0.01 cps and 10 me. Accuracy of setting 
of the order of one per cent is attainable with ordinary components. Frequency 
stability is of the order of 2 per cent per db of tube gain variation. 


Semi-Conductor Surface Phenomena.* W. H. Brarratn!. Semi-Con- 
ducting Materials, WH. 1K. Hentscn, ed., pp. 37-46. Proceedings of a con- 
ference held at the University of Reading (July 10-15, 1950) London, 
sutterworths, 1951. 

Developments in the understanding and interpretation of phenomena occur- 
ring at the surface of a semi-conductor are reviewed. The development starts 
with the Mott-Schottky theory of the space charge layer. Bardeen’s concept of 
a space charge layer due to ‘surface states’ explained the independence of recti- 
fieation on contact potential and Meyerhoff’s small values for the contact poten- 
tial between n- and p-type silicon, Shockley and Brattain observed that this 


contact potential difference increased with impurity concentration. Tllumination 


of a silicon surface produced hole and electron pairs in the space charge layer. 
The potential of the surface changed until the photocurrent was balanced by a 
conduction current. The relation between photo-current and potential change 
was of the same form as a forward characteristic for a rectifying contact. In an 
experiment similar to Becquerel’s, using water as an electrolyte, the surface 
may be biased in the reverse direction. When so biased the response to modu- 
lated light gives the differential resistance of the space charge layer. This re- 
sistance increases rapidly with reverse bias and the time constant of the layer 
increases, both agreeing qualitatively with theory. This experimental method 
of measuring changes in surface potential caused by illumination permits deter- 
mination of the properties of the space charge layer at the free surface of a semi- 
conductor, 


The Calculation oO} Traveling Wave-Tube Gain.* C. C. Cuter. Proc. 
TRE., 39, pp. 914-917, Aug., 1951. 


A reprint of this article may be obtained on request 
t Bell Tel. Labs 
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Essential information for claculation of traveling-wave-tube gain is sum- 
marized and condensed in this paper. The important relations are documented, 


presented in a concise form for simplified computation, and developed as a 


nomograph. The conclusions have been found to be in agreement with measure- 
ments on six different tube designs. 


Thermal Variation of Young’s Vodulus tm =SOME Fe Vi Mo Alloys.* 
M. E. Fine! and W. C. Exuuis'. References. Jl. Metals, 3, pp. 761-764, 


Sept., 1951. 


A Broad-Band Transcontinental Radio Relay System.* ‘TV. J. Grueser' 
and A. C. Perrerson'. Elec. Eng., 70, pp. SLO-S815, Sept., L9SL. 

Spanning the continent from coast to coast, this microwave relay system 
provides six channels, each of which can carry one television circuit or hundreds 
of telephone circuits. Some features of this vast network are descirbed. 


An Improved Telephone Set.* A. H. INGuts' and W. L. Turrnenw'. 
Elec. Eng., 70, pp. 770-775, Sept., 1951. 

The familiar telephone set has undergone numerous changes which will pro- 
vide better service at lower cost than do present models. Increased transmitting 
and receiving gain, better sidetone control, broader frequency response, faster 
dialing, simple ringing control, and a trim appearance are some of the features 
of the new design. 


The Institutes for Basic Research—Their Contribution to National 
Strength. M. J. Keviy', pp. 11-23. Applied Research is Not Enough, 
(booklet). Addresses at the Dedication of the Institutes for Basic Re- 
search, The University of Chicago, May 16, 1951. 


The Crystal Clock, W. A. Marrison'. Science Marches On, JAMES 
STok.ey, ed., N. Y., Ives Washburn, Inc., 1951, pp. 303-309. 


Observations of Zener Current in Germanium p-n Junetions. K. B. 
McAree!, E. J. Ryper'’, W. SHockiey', and M. Sparks'. Letter to the 
Editor. Phys. Rev., 83, pp. 650-651, Aug. 1, 1951. 


Experimental Radio-Telephone Service for Train Passengers.* N. 
Monk’. Proc. 1.R.E., 39, pp. 873-881, Aug., 1951. 

Experimental public radio-telephone service for train passengers was in- 
augurated by the Bell Telephone System several years ago. Initial installations 
operated in conjunction with a series of urban mobile base stations. More re- 
cently, highway mobile systems have been used for this service, and this paper 


describes a typical train installation operating through a highway channel. All 


*A reprint of this article may be obtained on request. 
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of these early systems utilized an attendant on the train. The cost of providing 
an attendant has, in some Cuses, been found excessive. Consequently, experl- 
ments have been initiated in which a coin box is used on the train. The arrange- 


ments for this purpose are also described. 


The Magneto-Resistance hi ffeet in Oriented Single Crystals of Ger- 
manium.* G. L. Pearson! and H. Sunw!. Phys. Rev., 83, pp. 768-776, 
Aug. 15, 1951. 


This paper describes an extensive study of the magneto-resistance effect in 


germanium as a function of crystal orientation. Experimental measurements 
establish the constants involved in the dependence of the effect on orientation 
of magnetic field and electric relative to the crystal axes. The measurements 
ire internally consistent with existing phenomenological theory based on cubic 
crystal symmetry, in which terms involving the magnetic field to higher than 
the second order are neglected. It is shown that such deviations as do occur 
wise from higher terms in the field, since an extension of the phenomenological 
theory to the fourth order predicts their symmetry. Relations are established 
between the experimentally observed phenomenological constants and those 
constants appearing in existing magneto-resistance electronic theories. It is con- 
cluded that no electronic theory yet worked out is entirely consistent with 
experiment. The present electronic theories are special cases of a very general 
theory recently proposed by Shockley, and it is possible that agreement can 
be obtained as soon as the computational difficulties of the latter theory are 


overcome 


New Phenomena of Electronic Conduction in Semi-Conductors. W. 


SHock.ey!. Semi-Conducting Materials, H. IX. HeNiscu, ed., pp. 26-36. 
Proceedings of a conference held at the University of Reading (July 
10-15, 1950) London, Butterworths, 1951. 

The semi-conductors silicon and germanium may be discussed as insulators 
the electronic structure of which is disturbed. Excess electrons, which act as 
negatively charged current carriers, may be present as may be ‘holes’ or places 
where electrons are missing from the valence-bond structure. Holes act as 
positively charged current carriers. In ordinary electronic conduction the flow 
of current carriers is substantially incompressible so that the density of carriers 
remains constant. When a new transistor phenomenon known as ‘carrier injec- 
tion’ occurs, however, the total density of holes and electrons may be greatly 
increased and this modulation of the electronic structure may be used both for 
scientific investigation and for practical amplification. In particular, carriers 

injected ata predetermined time and place into a known uniform elec- 
and their transit time to another place accurately timed by detecting 
val with a ‘collector’ point. Drift velocities and mobilities may be 


ly in this way with a directness unattainable by pre-transistor 
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methods involving conductivity and Hall effect. These new experiments and 
their related theories may furnish a foundation for a new engineering science of 
transistor electronics. 


Domain Patterns on Nickel.* H. J. Wituiams' and J. G. WaLKer'. 
Phys. Rev., 83, pp. 634-636, Aug. 1, 1951. 

Domain patterns have been observed on two single crystals of nickel cut in 
the form of hollow parallelograms. The length of the sides were parallel to the 
(111) directions in one specimen and to the (110) directions in the other. The 


crystals show domain structures with the three types of domain boundaries 


which are to be exper ted from a material hav ing the directions of « asy magneti- 


zation along the (111) directions. Domain boundary mc vement under the influ- 
ence of an applied magnetic field was observed. 


Polymorphism in Potassium Niobate, Sodium Niobate, and Other ABO; 
Compounds.* FE. A. Woop'. Bibliography. Acta Cryst., 4, pp. 353-362, 
July, 1951. 

The first part of this paper presents the results of optical and X-ray studies 
of the perovskite-type crystals, potassium niobate and sodium niobate. Potas- 
sium niobate is orthorhombic at room temperature, changing to tetragonal at 
about Zao. and cubie near 435°C. Sodium niobate Is orthorhombi at room 
temperature, changing to tetragonal at about 370°C. and to cubie at about 
640°C. 

The second part of the paper discusses relations among the structures of the 


ABO; compounds. 


Subjective Sharpness of Additive Color Pictures.* MM W. BALDWIN', 
Proc. I.R.E., 39, pp. 1173-1176, Oct., 1951. 
This is a report on the first numerical results to come from 


a laboratory 


experiment on the subjective sharpness of additive three-color pictures. The 
sharpness factor is isolated by using out-of-focus projection (of slides) instead 
of actual television transmission. 

An observer’s acuity for defocus is greatest for the green component and least 
for the blue component, in an additive three-color picture. When the same 
picture is reproduced in monochrome (white, red, green, or blue) at the same 
brightness, the observer’s acuity for defocus is equal to that found for the 
green component. 


Frequency-Modulation Terminal Equipment for the Transcontinental 
Relay System.* J. G. Cuarrer! and J. B. Maaeto'. Elec. Eng., 70, pp. 
880-883, Oct., 1951. 

To meet the exacting requirements of the new transcontinental microwave 
relay system, specially designed frequency-modulation terminal equipment was 


*A reprint of this article may be obtained on request. 
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constructed. The terminal transmitter converts either message or television 
signals to a frequency-modulated signal centered on 70 me and the terminal 
frequency-modulation receiver recovers these signals, thus providing a link 
between the relay system and other telephone facilities. 


Observer Reaction to Low-Frequency Interference in Television Pictures.* 
A. D. Fowuer'. Proc. [.R.E., 39, pp. 1332-1336, Oct., 1951. 

This paper presents results of tests to determine how much low-frequency 
interference can be tolerated in black-and-white television pictures. Various 
levels of single low-frequency interference were superimposed on a locally trans- 
mitted television picture. Observers viewed the picture and rate the disturbing 
effect of each level of the interference. Ratings were made in terms of preworded 
comments ranging from ‘not perceptible” to “unusable.” Interfering frequencies 
from 48 to 90 cycles per second were employed. 

Just visible interference appears as a flicker. The rate of flicker is the differ- 
ence between interfering and 60-cycle field frequencies. The most disturbing 
interference produced a flicker rate of 5 or 6 cycles per second. To be tolerated, 
peak-to-peak amy litude of this interference had to be 54 db weaker than the 
peak-to-peak amplitude of the television signal (including synchronizing pulse). 
For flicker rates of 0.5 and 12 cycles per second, the amount of interference 
which could be tolerated was larger by 14 and 3 db, respectively. 


Arcing at Electrical Contacts on Closure. Part IT. The Initiation of an 
Arc.* L. H. Germer'. Jl. Appl. Phys., 22, pp. 1183-1139, Sept., 1951. 

The capacity of the plates of an oscilloscope charged to 35 or 40 volts is dis- 
charged repeatedly by approaching electrodes of carbon, active silver, and inac- 
tive silver. Facts about the discharges, which are ares of very short duration, 
are inferred from resulting open circuit potentials and calculated electrode 
separations. 

The separation at the first arc varies in different experiments but corresponds 

i the average to 2 nominal electric field of 0.6 * 106 volts/em for carbon or 

ilver and to 2 * 10% volts/em for inactive silver. Each are is initiated 

very small number of field emission electrons. The hypothesis that a single 

electron may perhaps be sufficient is consistent with observations at later stages 
of each closure when the electrodes are closer and the field much higher. 

The earlier observation, that the potential across a short are is constant and 
independent of current, is not true if the are time is sufficiently short. For 
active silver a time comparable with 2 X 10> see is required to establish the 
steady are voltage characteristic of later stages of ares which last longer than 
this. The initial time during which the potential is decreasing toward its final 
steady value is 100 times the transit time of a silver ion across the gap. 
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Computation of Control Limits for p-Charts When the Samples Vary in 
Size. H. L. Jones®. Ind Quality Control, 8, pp 26-27, Sept _ 195] 


The Design of Switching Cirewts. W. Keisrer', A. E. Rirenie', and 
S. H. Wasusurn!. NV. Y., Van Nostrand, 1951. 556 pp. (Bell Telephone 
Laboratories Series). 

This is the first published textbook in its field. It presents, first, the funda- 
mental design principles of switching circuits composed of discrete-valued 
switching elements. Most of the discussion concerns two-valued elements, with 
greatest emphasis placed on electromagnetic relays. Chapters cover basic circuit 
paths, the logical interpretation of requirements, and the techniques of con- 
structing networks to fulfill these requirements. The symbolic methods of 
Boolean algebra and its application to the design of combinational and sequen- 
tial circuits is covered. Later chapters cover various unifunctional circuits such 
as selecting, connecting, translating, counting, and lockout. Final chapters dis- 
cuss methods of synthesising unifunctional circuit building blocks into larger 
circuits and systems 


Measurement of the Elastic Constants of Silicon Single Crystals and 
Their Thermal Coefficients. H. J. MceSkimin', W. L. Bonp', E. Burner’, 
and G. K. Tea’. Letter to the Editor. Phys. Rev., 83, p. 1080, Sept. 1, 
LOL. 

Interest in the properties of silicon single crystals arising from their use as 
semiconductors has led us to make measurements of the elastic constants of 
two single crystals. Measurements of velocities of propagation for both shear 
and longitudinal waves were made in the crystals as described in a recent paper 
by McSkimin. Frequencies in the range 8-12 me/see were used. 

The three independent elastic constants were evaluated, a density of 2.331 
(measured by pycrometer) being used. Data and formulas used are summarized 
in Table I. Two crystals were measured-—as indicated--with data obtained 
from the larger one being used to determine the elastic constants. Check meas- 
urements were made for the smaller crystal; and despite the less accurate ‘‘pulse 
overlap” technique used for two of the measurements, velocity agreement to 
within 0.15 per cent was obtained. 

Both crystals were of a high degree of crystalline perfection as shown by 
etching and X-ray tests 

Domain Wall Relaxation in Nickel. W. P. Mason’. Letter to the 
Kditor. Phys. Rev., 83, pp. 683-084, Aug. 1, 1951. 

Phase Transitions in Ferroelectrics. B. Marrutas'. National Research 
Council, Comm. on Solids. Phase Transformations in Solids. Ed. by 
R. Smoluchowski, J. E. Mayer, W.A. Weyl N: Y.. Wiley, 1951. 660 pp 
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e name ferroelectrics are classified those materials which exhibit 
LnoOmaves phenomenologicalls similar to the magnetic hehavior of 
magnetics. Perhaps it would have been more logical to use the term 


electrics, thus emphasizing the similarity in the dielectric behavior to 


known ferroelectrics are listed first, and then there follows 


n of the various theories which have been created to explain them. 


Data on Random-Noise Re quirements for Theater Television.* pP 

Mena’. Jl. SMP.TE., 67, pp. 89-107, Aug., 1951. 
Provisional evaluation of permissible random noise for theater television is 
onsidered from several sources of information. These cover broadcast television 
experience and the graininess in motion picture film; the requirements deduced 
from the various sources generally agree. For broadcast television, a frequency 
weighting and limit on weighted noise power have been used. The finer picture 
detail of theater television presumes a lower permissible random noise. Changes 
in weighting curve are discussed. A limit figure of noise is suggested, which is 
ble to graininess effects in motion pictures, though slightly more severe 


than present published performance on camera tubes. 


1 Spatial Harmonic Traveling-Wave Amplifier for Six Muillimeters 
Wavelength.* S. MituMaAn!. Proc. I.R.E., 39, pp. 1035-1043, Sept., 1951. 
This paper describes a traveling-wave amplifier in which the electren beam 
interacts with a spatial harmonic of an electromagnetic wave propagating along 
array of resonator slots. The result is a considerable reduction in operating 

m voltage for a given physical separation of the circuit elements. This type 
mplifier operating at about 1,200 volts has J ielded net power gains of about 

IS db in the 6-mm wavelength region. A magnetic field of about 1,600 gauss is 
sufficient for proper beam focusing. Aside from small variations of gain with 


frequency that is caused by internal reflections, the bandwidth is of the order 


of 3 per cent 


Form of Transient Currents in Townsend Discharges with Metastables.* 
J.P. Mounar'. Phys. Rev., 83, pp. 933-940, Sept. 1, 1951. 

The form of the current is caleulated for a Townsend discharge stimulated 

pulsed light beam, with particular reference to the current component 

by metastable effects. The caleulation is directed particularly to the 

ment of methods for quantitative interpretation of current patterns 


1 EX] rimentally. 


es of ¥ Proce SSES vo} Klectron I-mission Employing Pulsed Town- 
ond D charges on a VWillisecond Time Seale.* J P. MoLNari. Phys. 
Rev., 83, pp. 940-952, Sept. 1, 1951. 
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The relative amounts of electron emission from the cathode in a Townsend 


discharge caused by ions, photons, and metastables have been studied exper- 


mentally for several cathodes in argon, using pulsed-light stimulation of the 


discharge. The current initiated by metastables exhibits a slow build-up and 
, thus permitting easy separation from the faster rising effects of gas 


ionization and electron emission by photons and ions. Time constant stu 


decay 


lies of 


the slow component y ielded a diffusion constant for metastable argon atoms of 


45 cm? sec™ at one millimeter pressure, The efficiencies of electron emission by 


metastables and ions was found to be closely the same, while the quantum 


vield for photon emission was found to be generally smaller 


klectrical Prope rties of aFe.O3 and akFe.O0 Containing Titanium 
Morin'. Phys. Rev., 83, pp. 1005-1010, Sept. 1, 1951. 
Electrical conductivity, Hall effect, and Seebeck effect | 


1 nave been measured 
on two sets of polycrystalline samples of aFe:O; and aFe,O; containing from 
0.05 to 1.0 atomic per cent titanium (n-type impurity). One set of samples con- 
tained 0.6 atomic per cent excess of iron (n-type impurity), the 


» second set con- 
tained 0.6 atomic per cent deficiency of iron (p-type impurity 


The conductivity of pure aFe.O; is independent of this amount of stoichio- 


metric deviation. The slope of the log conductivity vs reciprocal temperature 
plot is 1.17 ev and the intercept at 1/T 0 is 2.1 & 10¢ ohm 
temperature conductivity varies from 

pure ake.OQ; to 0.3 ohm em! for ake,O 
titanium. 


ecm Room 
extrapolated) lo! 
containing 1.0 atomic per cent 


The measured Hall voltages seem to result entirely from magnetization of 


the samples, which are weakly ferromagnetic, and disappear above the ferro- 
magnetic Curie temperature. 

The temperature variations of the Fermi level are determined from Seebeck 
data. The temperature variations of carrier concentration are determined fi 
Fermi level and of mobility from carrier concentration and conductivity for 


] 
| 
some samples. Carrier concentration results indicate that each added titanium 
ion donates approximately one electron to the conduction process. Mobilities 


are found to be less than 2.0 em?/volt sec, suggesting that conduction lves 


> 


electrons in the d level of iron. 


Acceptance Inspection of Purchased Material.* J. E. Paumer? : 
G. D. Patrerson'. [nd Quality Control, 8, pp. 15-19, Sept , 195] 


This paper describes some of the principles and procedures employ 
in the form of components or finished pr 


inspection of purchased material 
gvely with procedures used in the Be ! Sys- 


The authors’ « xperience has been lar Li 
tem, and the illustrations have therefore been « 


rawn from this source. It is felt 


however, that considerations leading to the choice necifie inspection tech- 


. 
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niques will be generally applicable even though the number and volume of items 
purchased and the number of suppliers involved may in some cases differ widels 
In this presentation, stress has been placed on a discussion of the broader gauge 
factors underlying the engineering planning of inspection procedures rather than 


on specific sampling and control techniques. 


Analysis of Audio-Frequency Atmospheries.* R. WK. Porrer'!. Proc 
I.RA., 39, pp. 1067-1069, Sept., 1951. 

Sound portrayal techniques used in studies of speech and noise reveal the 
structure of atmospheric disturbances well known to long-wave radio and ocean- 
cable engineers as ‘‘whistlers,”’ ‘‘swishes,” and tweeks.” It is suggested that 
renewed investigation of these effects, using modern analyzing tools, might yield 


information of considerable scientific interest. 


Reflection of Llectromagnetic Waves from Slightly Rough Surface s.* 
8S. O. Rice!. Communications on Pure and Applied Math., 4, pp. 351-378, 
Aug., 195] 

Color Television and Colorimetry *W.T. WINTRINGHAM!. Proc I.RE., 
39, pp. 1135-1172, Oet., 1951. 

The high lights of the history of color measurement and of color photography 
ire reviewed. Following this introduction, the principles of modern three-color 
colorimetry are developed from a hypothetical experiment in color matching. 


The conventional theory of ‘“‘perfect color reproduction” by color television is 


built up from colorimetric background. Some of the difficulties to be expected 


in applying colorimetry to color television are brought out 
Finally, there is some discussion which tends to show that colorimetry may 
not be a sufficiently powerful tool to provide answers to all of the questions which 
will arise in the reproduction of scenes in color by television. The advantage of 
colorimetry as a background is indicated, however. 
\ reprint of this article may be obtained upon request. 
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